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This  contract  covered  the  period  1  January  1976  through  12  December  1980. 

The  period  1  January  1980  to  12  December  1980  represented  a  no  cost  extension. 

During  1980,  we  completed  the  proposed  studies  of  collislonal  redistri¬ 
bution  in  the  thallium-argon  system.  The  work  Included  experimental  and 
theoretical  studies  of  both  spontaneous  and  stimulated  light  scattering  by 
atoms  in  a  collislonal  environment.  The  principal  work  included  studying  the 
effects  of  radiative  trapping  on  the  shape  of  the  spontaneous  redistribution 
line  profile,  and  developing  and  applying  a  theoretical  model®  to  the  data 
previously  obtained  on  stimulated  Raman  scattering  (SRS)  and  stimulated 

collision-induced  fluorescence  (SCF).  The  major  results  on  Tl~Ar  summarized 

A 

here  are  discussed  in  detail  in  the  Ph.D.  thesis  of  H.  Raymer.* 

Previously  we  had  published^  the  first  quantitative  comparison  between  a 
collislonal  redistribution  profile  and  an  emission  profile  (conventional  line 
broadening  study  of  B.  Cheron,  R.  Scheps,  and  A.  Gallagher,  Phys.  Rev.  A  15, 
651  (1977)).  For  redistribution  in  thallium  we  found  generally  good  agreement 
between  the  shapes  of  the  profiles  obtained  by  the  two  experimental  methods, 
but  an  uncertainty  remained  as  to  the  absolute  scale  (or  calibration)  of  our 
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method.  [As  Indicated  at  the  end  of  this  report,  good  agreement  is  to  a  large 
extent  unique  for  thallium  which  has  an  Sj/2  upper  state.  We  now  know  that  in 
general  scattered  profiles  and  emission  profiles  are  only  the  same  under  quite 
restrictive  circumstances  (like  angular  and  polarization  averaging).]  The 
uncertainty  in  absolute  scale  resulted  from  the  presence  of  radiative 
trapping.  We  observed  as  our  signal  photons  emitted  at  535.0  nm  to  the  P3/2 
level  of  thallium.  Since  many  ground-state  thallium  atoms  were  in  the  line- 
of-sight  of  our  detection  optics  (along  the  laser  beam,  in  the  backward 
direction),  photons  which  were  initially  emitted  at  377.6  nm  could  be 
reabsorbed  (trapped)  and  subsequently  emitted  at  535.0  nm,  thus  causing  our 
detected  signal  at  535.0  nm  to  be  anomalously  large.  This  effect  could  give 
rise  to  at  most  an  Increase  by  a  factor  of  two,  since  the  two  resonance  lines 
have  nearly  equal  oscillator  strengths.  After  further  experimentation  we 
discovered  that  optical  pumping  by  the  laser  caused  a  significant  population 
in  the  ^3/2  *evel>  which  tended  to  reduce  the  detected  signal  at  535.0  nm  by 
radiative  trapping.  The  competition  between  these  two  trapping  effects  was 
studied  by  observing  emission  in  both  the  backward  and  side  directions,  and  it 
was  concluded  that  the  two  effects  largely  cancelled  each  other,  and  that  the 
good  (30Z)  reported^  agreement  between  the  redistribution  and  emission 
profiles  was  somewhat  fortuitous,  though  understandable. 

In  our  research  on  stimulated  scattering,  we  made  progress  in  the 


analysis  of  existing  data  on  SRS  and  SCF.  We  further  developed  the  theory  of 

stimulated  scattering  by  pulsed  laser  sources  to  account  for  these  data.  The 

3  WC 

theoretical  work  was  initiated  to  explain  a  discrepancy  reported  earlier  in  r, 
the  SCF  turn-on  behavior,  associated  with  the  transient  nature  of  the 


excitation.  The  theoretical  model  which  was  developed  proved  to  be  capable  of 
explaining  several  other  important  aspects  of  the  stimulated  scattering. 
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including  argon  buffer  gas  pressure,  laser  detuning,  thallium  number  density, 
and  laser  bandwidth.  [Again,  details  of  these  studies  are  described  in  M. 
Raymer's  thesis.1] 

During  1980,  when  the  final  stages  of  this  contract  were  being  completed, 
work  by  the  co-principal  investigator  (K.  Burnett,  J.  Cooper,  R.  Ballagh  and 
E.  W.  Smith,  Phys.  Rev.  A  22_,  2005  (1980);  K.  Burnett  and  J.  Cooper,  Phys. 

Rev.  A  22,  2027  (1980),  and  Phys.  Rev.  A  22,  2044  (1980))  indicated  that  the 
analysis  of  redistribution  experiments  in  terms  of  absorption  (or  emission) 
profiles  as  performed  earlier  in  this  contract  (refs.  2  and  5,  in  particular) 
can  often  be  erroneous  and  is  at  best  misleading.  Effects  of  Zeeman- 
degeneracy  mean  that  the  polarization  and  angular  distribution  of  the 
scattered  radiation  is  very  complicated  and  the  concept  of  profiles  has  to  be 
generalized.  (This  has  recently  been  confirmed  experimentally  for  the  Sr-Ar 
system  —  P.  Thomann,  K.  Burnett  and  J.  Cooper,  Phys.  Rev.  Lett.  45 ,  1325 
(1980)).  Fortunately,  for  the  thallium-argon  system  the  upper  state  is  an 
Sjy2  state  ^  isotropic  scattering  and  a  unique  comparison  with  the 
absorption  profile  is  in  this  case  possible  (but  not  for  non-S-states ! ) . 

Thus,  although  the  major  work  on  Tl-Ar  was  completed,  the  final  period  of 

the  contract  was  devoted  to  a  further  study  of  Zeeman  degeneracy  effects 

(especially  for  high  power  lasers  following  J.  Cooper,  R.  J.  Ballagh  and  K. 

Burnett,  Phys.  Rev.  A  22,  535  (1980)).  The  results  of  this  w  >rk  (see  ref.  8) 

showed  that  the  scattering  was  a  triplet  for  polarizations  ,  el  to  the 

laser  but  a  doublet  for  polarizations  perpendicular  to  the  laser  —  not  just  a 

3 

triplet  as  initially  supposed. 

The  ultimate  conclusion  is  that  in  the  future  in  order  to  understand 
radiative  collisions,  in  particular,  and  the  effects  of  laser  interactions 
with  atomic  systems,  in  general,  we  will  need  to  know  much  more  detail 
concerning  Zeeman  degeneracy  and  specific  interatomic  interactions. 
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This  paper  reports  on  our  studies  of  near-resonant  scattering  of  laser  light  in  a  coliisional  environment.  A 
pulsed  dye  laser  with  a  peak  power  of  S3  MW/cm1  was  tuned  near  the  460.73-nm  resonance  line  of  strontium 
and  the  side  emission  was  observed  from  an  oven  containing  both  strontium  vapor  and  argon  buffer  gas.  The 
emission  was  con- posed  of  three  spectral  components:  Rayleigh  scattering  at  the  frequency  of  the  laser  o>L, 
fluorescence  at  the  resonance  frequency  oig  of  strontium,  and  a  third  component  at  2tnt  —  These  three 
components  have  been  studied  as  a  function  of  the  frequency  and  intensity  of  the  laser  and  also  as  a  function 
of  the  argon  buffer  gas  pressure.  While  the  Rayleigh  emission  was  found  to  vary  as  A~5  (A  =  to0— o>L),  the 
fluorescence  component,  which  was  produced  by  Sr-Ar  collisions,  was  found  to  be  asymmetric  with  the  sign 
of  A  as  predicted  by  line-broadening  theory.  By  measuring  the  ratio  of  the  intensities  of  the  fluorescence  and 
Rayleigh  components,  we  were  able  to  measure  directly  the  coliisional  redistribution  function,  important  in  the 
study  of  radiative  transfer  in  stellar  and  planetary  atmospheres.  At  high  laser  intensities  all  three  components 
were  found  to  saturate.  The  results  were  compared  with  the  theoretical  predictions  of  Mollow's  steady-state 
theory.  Theoretical  fits  for  the  high-intensity  results  were  obtained  when  the  coitisional  cross  sections  were 
taken  to  be  considerably  smaller  than  in  our  low-intensity  measurements.  We  believe  the  discrepancy  lies  in 
the  use  of  a  steady-state  theory  for  a  transient  experiment.  Effects  of  radiative  trapping  and  spatial  averaging 
are  also  discussed. 


I.  INTRODUCTION 

In  two  previous  Letters1"’  we  reported  on  the 
first  observation  of  the  spectral  resolution  of  light 
scattered  by  an  atomic  vapor  near  a  resonance 
line,  and  on  the  first  experimental  study  of  the 
spectral  redistribution  of  light  by  collisions  as  a 
function  of  the  detuning  of  the  incident  light.  This 
paper  is  a  summary  and  extension  of  these  studies, 
giving  details  of  the  previous  measurements  and 
reporting  later  studies  of  saturation  behavior. 
Basically,  in  all  of  our  measurements  we  illumin¬ 
ate  a  cell  containing  Sr  vapor  (10*4-10*1  Torr)  and 
Ar  buffer  gas  (10-500  Torr)  with  a  pulsed,  tune¬ 
able  dye  laser  whose  wavelength  is  near  the 
460.73-nm  (lP°  -  lS0)  resonance  line  of  Sr,  and  we 
analyze  the  spectrum  of  scattered  light.  In  this 
Introduction  we  would  like  to  put  our  measurements  . 
into  perspective  by  pointing  out  their  relation  to 
other  studies. 

When  a  monochromatic  light  wave  impinges  on 
an  atom  near  one  of  its  resonance  (absorption) 
lines,  part  of  the  light  is  scattered.  In  the  ab¬ 
sence  of  collisions  and  af  low  intensities,  the  spec¬ 
trum  of  the  scattered  light  is  a  single  line  centered 
on  the  incident  wavelength,  known  as  Rayleigh 
scattering.  If  the  lower  state  is  the  ground  state, 
the  scattered  light  is  monochromatic  for  a  station¬ 
ary  atom1  and  Doppler-shifted  for  a  moving  atom/ 

If  the  lower  state  decays,  the  spectrum  of  the  scat¬ 
tered  light  is  broadened. **"  If  the  atom  that  scat¬ 
ters  the  light  collides  with  other  atoms,  the  scat¬ 
tered  light  has  a  second  component  emitted  near 


the  resonance  frequency  that  is  known  as  fluore¬ 
scence.  In  a  simple-minded  description  this  com¬ 
ponent  is  due  to  the  absorption  of  the  incident  light 
in  the  collision-broadened  wing  of  the  absorption 
line  and  its  subsequent  emission  near  line  center. 
The  theory  of  coliisional  redistribution  has  been 
developed  from  basic  principles  (the  Kramers- 
Heisenberg  dispersion  formula,  and  collision 
dynamics)  by  Omont,  Smith,  and  Cooper8  and 
Huber9  and  is  applicable  in  the  central,  impact 
region  of  the  spectrum.  No  explicit  theory  for  re¬ 
distribution  exists  outside  this  region.  The  know¬ 
ledge  of  the  spectral  and  spatial  redistribution  of 
light  scattered  by  atoms  in  the  presence  of  col¬ 
lisions  is  of  fundamental  importance  in  the  study 
of  radiation  transfer  in  stellar  and  planetary  at¬ 
mospheres/’  10  In  fact  we  found  it  quite  surprising 
that  no  direct  measurements  of  the  redistribution 
function  had  been  made  earlier. 

Low -intensity  absorption  and  emission  measure¬ 
ments  have  been  done  for  many  years  in  atomic 
vapors  broadened  by  a  foreign  gas.  In  fact,  both 
the  impact  and  the  quasistatic  theories  of  line 
broadening"'11  were  developed  to  deal  with  these 
absorption"  and  emission19  spectra.  It  is  intuitive¬ 
ly  clear  that  the  "excitation  functions”  (i.e.,  the 
dependence  of  the  intensities  of  the  various  spec¬ 
tral  components  of  scattered  light  on  the  incident 
frequency)  are  related  to  the  absorption  and  emis¬ 
sion  spectra.  In  the  impact  region  this  relation 
can  be  shown  explicitly  from  the  general  theory,8’9 
and  outside  the  impact  region  our  experimental 
results  suggest  a  similar  connection. 
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When  the  incident  light  is  strong,  even  {or  a 
noncolliding  stationary  atom  and  monochromatic 
incident  field,  three  components  appear.  This 
was  predicted  first  by  Mollow1'  and  studied  since 
theoretically  by  many4’7’17*22  who  obtained  essen¬ 
tially  the  same  results.  These  components  were 
observed  in  atomic  beam  experiments  using  very 
narrow  cw  dye  lasers,21*”  and  recent  results 
agree24  with  the  theory.  In  addition  to  the  Rayleigh 
component  mentioned  above,  there  is  a “three  - 
photon’’  component  that  corresponds  to  absorption 
of  two  incident  photons  and  the  emission  of  another 
photon  (a  resonant  Raman  effect).  Conservation 
of  energy  then  requires  the  emission  of  one 
“fluorescence”  photon  on  resonance.  As  the  in¬ 
tensity  is  raised,  all  three  components  eventually 
saturate. 

Scattering  at  high  light  intensity  including  the  ef¬ 
fects  of  collisions  has  been  studied  using  several  the¬ 
oretical  approaches.  Mollow24’ 27  introduced  phe  - 
nomenological  damping  terms  into  his  treatment 
of  the  collisionless  atom  scattering  in  a  strong 
steady -state  light  field.  This  procedure  is  ex¬ 
pected  to  give  correct  results  if  the  “optical”  cross 
section,  which  is  related  to  the  transition  prob¬ 
ability  per  incident  photon,  is  independent  of  the 
(strong)  incident  flux.  Another  approach,  to  some 
extent  complementary,  was  taken  by  Yakovlenko 
and  co-workers,24*30  who  studied  the  dynamics  of 
atomic  collisions  in  the  presence  of  a  classical 
(strong)  incident  field.  They  show  that  for  weak 
fields  the  absorption  reduces  to  the  impact  or 
quasistatic  limits,  implying  a  constant  (in  intensity) 
but  frequency -dependent  optical  cross  section. 

For  high  intensities  they  predict  a  change  of  the 
optical  cross  section.  A  similar  study  was  done 
by  Kroll  and  Watson31  and  Lau.” 

There  is  yet  another  point  of  view  of  the  scatter¬ 
ing  process  that  unifies  and  simplifies  many  cal¬ 
culations:  that  of  the  dressed  atom  and  dressed 
molecule.22’33  The  dressed  atom  state  is  the  (adi¬ 
abatic)  eigenstate  of  the  atom  in  a  strong,  single¬ 
mode  external  field— in  the  rotating -wave  approxi¬ 
mation.  Light  scattering  can  be  viewed  as  spon¬ 
taneous  emission  causing  transitions  among  the 
various  states.22,34  The  absorption  process  during 
collisions  can  be  viewed  as  a  simple  transition 
among  these  states:  the  ground  state  of  the  ab¬ 
sorber  before  the  collision  and  its  excited  state 
with  one  less  photon  in  the  field  after  the  collision. 
A  similar  view,  put  forward  originally  by  Yakov¬ 
lenko  el  al .24*30  provides  a  unified  picture  of  ab¬ 
sorption,  collision-induced  absorption,  and  radia¬ 
tive  (“optically  induced”,  "switching”)  col- 
lisions.1**30’3*’34 

A  number  of  recent  experiments  indicate  the 
importance  of  collisional  redistribution  of  light. 


Rousseau  el  al.  have  seen  two  different  time  de¬ 
pendences  in  their  pulse -excited  Raman  spectra 
in  iodine.37*30  Several  authors  have  analyzed  this 
experiment.34*  1°*,: 2  Recently  Li  ran 13  found  that  in¬ 
deed  the  shifted  (Raman)  component  follows  the 
laser  pulse  in  time  and  the  fluorescence  component 
decays  with  the  lifetime  of  the  upper  state,  as 
predicted.31’ 12  Collisional  redistribution  has  been 
observed  in  resonance  scattering  from  sodium 
in  high-pressure  helium,41  from  sodium  in  a  dense 
high-temperature  plasma, 15  and  also  from  stron¬ 
tium  vapor. 14  Collisional  redistribution  has  been 
studied  as  an  excitation  mechanism  in  cesium- 
rare-gas  mixtures.17  In  two-photon  absorption, 
a  “fast”  off -resonance  component  and  a  “slow" 
on-resonance  component  were  studied  by  Grisch- 
kowsky. 14  The  formation  of  molecules  that  occurs 
during  collisional  absorption  of  light  has  been 
utilized  for  spectroscopy  by  Gallagher13  and 
others19  for  many  years. 

In  Sec.  II  we  review  the  relevant  theory  as  it 
exists.  The  experiment  is  described  in  Sec.  in, 
and  the  results  are  analyzed  in  Sec.  IV.  Basically 
we  studied  the  three  emission  components  as  a 
function  of  laser  frequency,  laser  intensity  and 
buffer  gas  pressure.  There  are  also  separate 
sections  on  our  "auxiliary"  experiments:  studies 
of  radiation  trapping,  measurements  of  absorption 
and  emission  spectra,  and  the  measurement  of  the 
laser  spatial  profile.  Section  V  presents  a  brief 
summary  and  discussion. 

II.  THEORY 

A.  Low  laser  intensities:  collisional  regime 

Consider  an  atom  and  an  incident  monochromatic 
laser  field,  which,  at  the  position  of  the  atom,  is 

BU)  =  ttE  cos(u>£f) ,  (1) 

where  iL  is  the  polarization  vector,  o>x  is  the  In¬ 
cident  frequency,  and  £  is  the  electric  field. am¬ 
plitude.  The  incident  intensity  of  the  laser  in  units 
of  photons  cm*2  sec*1  is 

IL  =  (2) 

The  incident  light  frequency  u't  is  taken  to  be  near 
the  atomic  resonance  frequency  u0l  and  the  de¬ 
tuning  A  is  defined  by 

A  *  w0  -  .  (3) 

We  now  define  the  Rabi  frequency  on  resonance: 

n  «rf£/R,  (4) 

where  d  =  c(l \tL  •  r  |2)  Is  the  electric  dipole  matrix 
element  between  the  two  relevant  atomic  state*. 
Thus  the  laser  intensity  is  proportional  to  the 
square  of  the  Rabi  frequency,  IL  a  O2.  We  make 
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the  unrealistic  assumption  that  both  states  are 
nondegenerate.  When  a  lSn  state  is  coupled  by 
linear  polarized  radiation  to  a  */*  state  (as  in  our 
studies)  this  is  a  good  approximation  in  the  absence 
of  collisions.  In  presence  of  collisions  the  picture 
has  to  be  modified.90  We  also  assume  that  the 
lower  of  these  states  1 1)  is  the  ground  state,  and 
the  upper  state  1 2)  has  a  total  damping  rate  y 
(half  width  at  half  maximum)  given  by 


the  duration  of  the  collision  is  small  compared  to 
the  time  between  collisions.  In  addition,  it  is  ap¬ 
plicable  only  for  the  frequency  region  close  to  line 
center.  Mihalas'  has  shown  that  the  boundary 
frequency  A,,  beyond  which  the  impact  approxima¬ 
tion  fails,  is  given  (for  a  Van  der  Waals  interac¬ 
tion)  by 

Af*1.65»’s/Jw;/ye',\  (6) 


y  =  YM  +  >e  +  Vf  (5) 

The  unit  of  y  is  rad/sec.  The  rate  y„  is  due  to  the 
natural  radiative  decay  of  the  excited  state,  while 
yt  and  y,  result  from  elastic  and  inelastic  col¬ 
lisions,  respectively. 

These  damping  processes  are  shown  schematic¬ 
ally  in  Fig.  1.  The  driving  laser  field  is  repre¬ 
sented  by  the  solid  lines,  spontaneous  decay  by 
the  wavy  line,  and  collisions  by  double  lines. 

Elastic  collisions  are  usually  viewed  as  collisions 
which  interrupt  the  phase  of  the  emitted  wave 
train.  However,  in  the  context  of  our  work  an  elas¬ 
tic  event  can  also  be  considered  as  a  collision 
which  transfers  populations  from  the  laser-excited 
virtual  level  to  the  real  level.  In  this  sense  these 
events  are  actually  “quasielastic"  since  the  col¬ 
lisions  must  provide  or  absorb  the  energy  dif¬ 
ference  ®A.  In  contrast  to  this,  inelastic  col¬ 
lisions  (quenching  collisions)  are  those  events 
which  cause  a  net  population  transfer  from  the 
excited  state  to  the  ground  state  in  a  two -level 
system  or  to  nearby  levels  in  a  multilevel  sys¬ 
tem. 

The  theory  has  been  worked  out  in  detail  for  the 
impact  region  of  the  spectrum.®,'>,*,<*  The  impact 
theory  is  valid  only  for  pressures  low  enough  that 


where  Nt  is  the  density  of  the  perturbers  (in  our 
case  the  argon  atoms),  yt  is  the  half  width  at  half 
maximum  (HWUM)  of  the  broadening  due  to  elastic 
collisions,  and  r  is  the  mean  relative  velocity  of 
the  emitting  and  perturbing  atoms.  The  collisional 
linewidth  is  usually  determined  by  experiment,  or 
by  specific  theoretical  calculations.  For  example, 
in  Van  der  Waals  broadening* 

yt-NtoMv,  (7) 

where 

oi  =  4.04(C,/r)2/*.  (8) 

Here  C,  is  the  Van  der  Waals  interaction  constant. 

Inside  the  impact  region,  for  stationary  atoms 
excited  from  their  ground  state  by  weak  mono¬ 
chromatic  light,  the  intensity  distribution  of  the 
scattered  light  integrated  over  all  angles  and 
summed  over  polarizations,  in  steady  state,  in 
units  of  photons  cm-1  sec"1  is  given  by  Huber’s 
Eq.  (4.8),® 


1  (y„/ 2Wdw 


y«-T  T/[,  jl. 

L  L  y  f+Yt  (w0-«r+rJ  ’ 


(9) 


damping  rates.  The  straight  arrows  represent  the  laser 
field  at  frequency  while  the  wavy  arrow  represents 
radiation  damping  of  the  upper  level  at  the  spontaneous 
rate  y*.  The  double  arrows  represent  collisional  damp¬ 
ing:  elastic  collisions  take  the  atom  from  the  virtual 
level  (dashed  line)  to  the  real  level  at  a  rate  y*.  while 
inelastic  (quenching)  collisions  take  the  atom  from  level 
1 2)  to  level  1 1)  at  a  rate  y, , 


where  N  is  the  density  of  scatterers  (in  our  case 
the  strontium  atoms),  and  the  symbols  are  defined 
above. 

The  first  term  in  this  equation,  which  occurs  at 
the  incident  frequency  wt,  is  known  as  Rayleigh 
scattering  and  corresponds  to  the  process  shown 
In  Fig.  2(b).  The  incident  light  is  elastically  scat¬ 
tered  off  a  virtual  state  produced  by  the  laser 
(shown  as  a  broken  line).  Since  the  ground  state 
has  no  natural  width,  the  Rayleigh  scattering  has 
the  same  spectral  profile  as  the  excitation  source. 

The  second  term  in  Eq.  (9)  is  the  collision-in¬ 
duced  fluorescence  and  corresponds  to  the  process 
shown  in  Fig.  2(d).  The  laser-induced  virtual  state 
is  collisionally  transferred  to  the  real  excited 
state  which  then  fluoresces  with  a  Lorentzian  pro¬ 
file  centered  at  w0  and  has  a  width  equal  to  the  sum 
of  the  natural  and  collisional  widths.  In  the  ab¬ 
sence  of  collisions,  the  fluorescence  vanishes 
and  the  scattered  light  maintains  the  profile  of 
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ric.  2.  (a)  Two  atomic  levels  separated  by  the  unper¬ 
turbed  resonant  frequency  u>,  plus  the  laser-induced 
Stark  sh.'t  A;  ihl  Rayleigh  scattering  from  the  ground 
•late  at  the  laser  frequency  u»t;  (c)  Rayleigh  scattering 
from  the  excited  state;  (di  collision- induced  fluorescence 
at  •  O;  tel  “three-photon"  scattering  at  o>j=  2uia- 
t-,  *  At  and  consequent  fluorescence  at  + 


the  incident  light,  even  U  the  excitation  source 
is  narrower  than  and  tuned  inside  the  natural 
width  of  the  resonance  line.2**1 

In  an  atom  whose  upper  and  lower  levels  have 
spatial  degeneracy,  the  angular  distribution  of  the 
radiation  and  its  polarization  properties  are  com¬ 
plicated  (cf.  Omont  et  al.*).  In  our  experiments 
the  lower  level  is  a  ‘S0  state,  and  the  upper  one  is 
‘P,;  therefore,  the  Rayleigh  scattering  has  a  clas¬ 
sical  dipole  distribution,  but  the  resonance 
fluorescence  distribution  depends  on  the  details 
of  the  depolarization  due  to  collisions.90  In  our 
experiment  this  is  further  complicated  by  the  fact 
that  our  sample  is  optically  thick  on  resonance; 
therefore  the  light  that  reaches  our  detector  is 
multiply  scattered.  We  have  not  studied  these 
depolarization  effects. 

When  the  laser  is  tuned  far  off  resonance,  |  A  j 
»  y,  the  two  components  can  be  resolved  and  the 
spectrally  integrated  intensities  of  the  Rayleigh 
scattering  I„  and  collisionally  induced  fluorescence 

1,  are 


Thus  in  the  impact  region,  the  integrated  inten¬ 
sities  of  both  the  Rayleigh  scattering  and  fluor¬ 
escence  behave  as  A*2. 

Note  that  Eq.  (11)  predicts  that  the  irtegrated 
fluorescence  signal  Ir  should  be  linear  in  laser 
intensity  and  also  linear  in  the  density  of  pertur- 
bers  Np  (via:N,).  However,  the  linearity  holds 
only  at  low  laser  intensity  and  low  pressure.  In 


Sec.  IV  we  will  see  that  at  higher  laser  inten¬ 
sity  or  at  high  pressure,  saturation  occurs. 

As  |  A  |  becomes  larger  than  A,,  the  impact  ap¬ 
proximation  fails.  Theories  for  the  redistribu¬ 
tion  function  outside  the  impact  region  have  not 
been  conveniently  developed.  (Huber9  has  written 
down  the  integrals  in  terms  of  the  time  develop¬ 
ment  operator,  but  they  have  not  been  solved.) 
However,  since  for  j  Aj  collisions  of  the  scat¬ 
tering  atom  do  not  change  the  Rayleigh  intensity, 
we  expect  Eq.  (10)  to  hold  outside  the  impact  re¬ 
gion.  Indeed  far  from  line  center,  it  goes  over  to 
the  well-known  Rayleigh  scattering  formula.52 

While  collisions  do  not  affect  the  Rayleigh  in¬ 
tensity  far  from  resonance,  the- fluorescence  com¬ 
ponent  depends  heavily  on  the  collision  dynamics. 
Basically  it  depends  on  the  ability  of  the  collisions 
to  supply  (or  absorb)  the  energy  difference  between 
the  laser  photon  and  the  scattered  photon.  Indeed, 
far  from  resonance,  the  collision  can  no  longer 
provide  this  energy  difference  efficiently  and  the 
cross  section  is  expected  to  drop  off.12  Even  be¬ 
fore  the  fall -off,  non-Lorentzian  shapes  similar 
to  those  seen  in  emission  and  absorption  experi¬ 
ments  are  expected.  Accordingly,  we  will  replace 
yt  in  Eq.  (11)  by  the  "empirical"  elastic  collision 
rate  ye(A),  which  is  a  function  of  the  detuning  and 
will  be  indicative  of  a  non-Lorentzian  line  shape. 
The  cross  section  in  Eq.  (7)  also  becomes  a  func¬ 
tion  of  detuning,  denoted  by  <JE( A). 

Lisitsa  and  Yakovlenko50  have  shown  the  connec¬ 
tion  between  oE(A)  and  the  static  theory  of  line 
broadening.  It  should  be  emphasized  that  the 
optical  cross  section  oe(a)  is  not  a  cross  section 
in  the  usual  sense  (although  it  does  have  dimen¬ 
sions  of  cm2).  The  actual  rate  of  production  of 
excited  atoms  is  described  by  the  quantity  ot( A)02/ 
A2.  Therefore  it  is  proportional  to  the  incident 
laser  flux  in  dimensionless  units. 

B.  High  User  intensity;  colliaonlew  regime 

We  now  discuss  the  distribution  of  scattered 
light  when  the  laser  intensity  is  allowed  to  become 
large,  but  when  no  collisions  are  present  and 
Doppler  broadening  has  been  neglected.  Recently 
there  has  been  considerable  theoretical  interest 
in  the  spectrum  of  scattered  light  at  high  incident 
intensities.*’ ,,,e*22  When  fl2  +  A2 »  y2 ,  the  inten¬ 
sity  (photons  cm*3  sec*1)  in  steady  state  for  mono¬ 
chromatic  incident  light  is  given  by  Mollow**  as 


N(y*/2)fl3du>r  .,  v  0«/4  /  */t  \ 
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where  the  widths  of  the  Lorentzian  emission  pro¬ 
files  are  given  by 


\  £  (\ 

ln2+A2^ 

■J. 

(FTa5/ 

and  the  other  terms  are  defined  in  Sec.  IIA  with 
the  exception  of  6  which  is  the  ac  Stark  shift  given 

by 

«*A{[l+(fl/A)3]l/2_l}.  (14) 

Note  that  at  low  intensities  the  shift  becomes 

6  =  (0«|A|),  (15) 

which  is  the  well-known  result  derived  from  sec¬ 
ond-order  perturbation  theory.53 

There  are  four  components  in  the  emission  spec¬ 
trum  given  by  Eq.  (12).  The  first  two  terms  are 
centered  at  the  frequency  of  the  incident  laser. 
These  terms  are  identified  as  the  coherent  and  in¬ 
coherent  parts  of  the  Rayleigh  scattering,  respec¬ 
tively,  and  are  shown  at  uiL  in  Fig.  3.  The  inco¬ 
herent  Rayleigh  component  is  unimportant  at  low 
intensities  but  becomes  dominant  at  high  laser 
intensities.  The  third  and  fourth  terms  are  cen¬ 
tered  at  oi  =  -  (A  +  6)  and  =  u>t  +  (A  +  6).  These 

“side  bands,”  which  have  recently  been  observed 
experimentally  by  a  number  of  researchers, i,23‘25 
are  due  to  higher-order  photon  interactions.  The 
process  leading  to  these  side  bands  is  depicted  in 
Fig.  2(e).  With  the  laser  tuned  to  a>L  *  u>0,  two 
laser  photons  are  absorbed,  a  photon  is  emitted 
and  an  excited  state  produced.  This  “three-pho- 


FIG.  3.  Sketch  of  the  emission  spectrum  of  a  two-level 
•tom  strongly  driven  near  resonance  by  a  monochroma¬ 
tic  excitation  source  at  frequency  wt.  In  the  presence 
of  collisions,  the  spectrum  consists  of  the  four  compo¬ 
nents  shown  by  solid  lines:  coherent  (narrow)  and  incohe¬ 
rent  (broad)  Rayleigh  scattering  at  u>L,  "three- photon” 
emission  at  u,  and  resonance  fluorescence  at  w,  produced 
by  both  elastic  collisions  and  the  ,<thrce-photon”  pro¬ 
cess.  In  the  absence  of  collisions  the  fluorescence 
(broken  line)  has  the  same  magnitude  as  the  component 
•t  <•>). 


ton”  component  is  given  by  the  third  term  and  is 
shown  as  the  profile  at  u>3  in  Fig.  3.  Since  an  ex¬ 
cited  state  is  produced  for  every  photon  generated 
at  w3,  there  will  also  be  emission  at  u>0*6  from 
the  excited  atom.  This  is  the  origin  of  the  other 
side  band  which  is  shown  as  the  broken  profile  at 
uj0  in  Fig.  3.  Thus  the  two  side  bands  have  equal 
intensities.  When  the  excitation  laser  is  tuned  to 
resonance,  these  two  peaks  occur  as  side  bands 
located  at  wt  +  fl  and  Therefore  when  the 

laser  is  tuned  to  resonance,  one  sees  the  Rayleigh 
emission  at  wt  with  two  side  bands  each  shifted  by 
the  Rabi  frequency  from  the  resonance.  At  this 
point  it  is  useful  to  think  of  the  side  bands  as  re¬ 
sulting  from  the  modulation  of  the  resonance  fre¬ 
quency  uj0  by  the  Rabi  frequency  fl.22 

Having  discussed  the  higher-order  photon  inter¬ 
actions,  we  can  now  understand  the  dependence  of 
the  two  Rayleigh  components  on  laser  intensity. 

The  coherent  Rayleigh  scattering,  linear  in  laser 
intensity,  is  scattering  off  the  ground  state  by  a 
single  laser  photon.  The  incoherent  Rayleigh 
scattering,  cubic  in  laser  intensity,  indicates  the 
following  process:  Two  laser  photons  are  absorbed 
to  give  an  excited  state  as  in  Fig.  2(e)  which  Ray¬ 
leigh  scatters  a  third  laser  photon  as  in  Fig.  2(c). 
As  the  excited  state  has  a  width,  this  gives  rise 
to  the  width  s  of  the  incoherent  Rayleigh  scatter- 
ing. 22,54 

When  R2+  A2  »  y],  the  side  bands  can  be  resolved 
from  the  central  Rayleigh  scattering  and  Eq.  (12) 
can  be  integrated  over  &>,  giving  the  integrated 
Rayleigh  intensity  IR  (the  sum  of  the  coherent  and 
incoherent  components)  and  the  integrated  side 
bands  /,  and  IF,  due  to  the  “three -photon”  scatter¬ 
ing  and  the  subsequent  fluorescence,  where 

/  /i«\ 

A2  +  fi4  (16) 


/  _/  -  v  (y<r/8)Q4 

h  -  b  -  « ( “ {/)  (Ai + nV2)  •  (17) 

As  before,  the  Rayleigh  scattering  can  be  viewed 
as  emission  from  the  laser-excited  virtual  level. 
At  low  intensities  when  the  excited  state  popula¬ 
tion  is  small,  we  can  associate  a  population 
(N/2)[fl2/(A2  +  R2)]  with  this  virtual  level.  If  we 
let  this  “state”  radiate  at  a  rate  yN,  we  obtain 
Eq.  (16).  An  analogous  argument  can  be  made  at 
higher  laser  intensities  by  including  scattering  off 
the  excited  state. 

Equation  (16)  predicts  that  at  high  intensities  the 
Rayleigh  intensity  IR  will  saturate.  This  process 
can  be  viewed  as  an  equalization  of  the  population 
in  the  virtual  and  real  states.  It  is  a  manifesta¬ 
tion  of  the  limitation  on  the  spontaneous  scatter- 
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ing  intensity  by  the  coupling  constant  to  the  vac¬ 
uum,  y„. 

Equation  (17)  shows  that  at  low  intensities,  the 
“three -photon”  scattering  and  resulting  fluores¬ 
cence  are  proportional  to  the  square  of  the  laser 
intensity.  At  higher  intensities,  these  components 
also  saturate  due  to  population  equalization  of  the 
ground  and  excited  states.  Note  that  since  the 
three-photon  process  corresponds  to  the  genera¬ 
tion  of  a  photon  and  an  excited  atom,  the  process 
is  a  type  of  Raman  scattering  and  therefore  will 
have  positive  gain  for  a  noninverted  system. 
Stimulated  emission  of  the  component  at  u>3  has 
been  predicted55*5"  and  observed.50"53  Also,  the 
emission  of  /3  and  IF,  the  absorption  of  two  laser 
photons  and  the  return  of  the  atoms  to  the  ground 
state  is  a  four-photon  parametric  process  that  has 
been  discussed  in  the  literature.55 

C.  High  laser  intensity:  collisional  regime 

In  this  section,  we  consider  the  effects  of  high 
laser  intensities  when  collisions  are  also  present. 
A  number  of  papers  have  treated  collisions  through 
phenomenological  relaxation  terms,  extending  the 
treatment  of  the  collisionless  atom  interacting  with 
radiation.25*27*35*52*55  In  all  these  treatments  it  was 
assumed  that  a  strong  incident  electromagnetic 
field  has  no  direct  influence  on  the  rate  of  colli¬ 


sions,  yB\  the  main  effect  is  to  saturate  the  frac¬ 
tion  of  virtual  levels  and  the  population  difference 
in  the  two -level  system.  Recently,  collisions  of 
atoms  in  the  presence  of  the  incident  field  were 
treated  more  accurately,  and  it  has  been  pointed 
out  that  when  the  fields  get  strong  they  can  directly 
influence  the  collision  rates.25*32  However  the 
emission  spectrum  was  not  explicitly  computed. 
Therefore  we  have  chosen  to  treat  the  relaxation 
constant  y£  in  Mollow’s  treatment25* 27  as  a  free 
parameter. 

Since  we  are  primarily  concerned  with  scatter¬ 
ing  at  optical  frequencies,  the  major  collisional 
effect  is  due  to  elastic  collisions.  Inelastic  col¬ 
lisions  also  have  to  be  considered  but  are  gen¬ 
erally  less  frequent  as  they  require  large  momen¬ 
tum  transfers.  At  optical  frequencies  inelastic 
collisions  that  transfer  the  atom  from  the  ground 
state  to  the  excited  state  are  usually  negligible, 
since  KuJkT  is  large. 

We  quote  the  results  of  Mollow’s  calculation,25*27 
in  which  he  added  phenomenological  relaxation 
terms  to  the  equations  of  a  strongly  driven  two- 
level  system.  The  use  of  phenomenological  re¬ 
laxation  terms  gives  rise  to  Lorentzian  line  shapes 
and  is  therefore  only  applicable  in  the  impact  re¬ 
gion  of  the  spectrum  [cf.  Eq.  (6)].  When  Sl2  + A2 
»  y2,  the  scattered  spectrum  (in  photons  cm*3  sec*1) 
in  steady  state  is 
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where  the  parameter  i\  is  defined  as 

nm  (r*  «-y/+y*)/(2y#+ 2y,) .  (19) 

The  detuning  parameter  is  defined  as 
B'=A  +  6,  (20) 

and  the  widths  of  the  profiles  are 
,  %/2»in2  +  2A2\ 

(21) 

,  ,  ,/(l  +  »j)02+27jAM 

S,  =  (r#+y/)(- — — ) . 

The  four  terms  in  Eq.  (18)  are  analogous  to  the 
terms  given  in  Eq.  (12)  for  the  collisionless  case. 
Therefore  we  will  refer  to  these  terms  as  coher¬ 
ent  Rayleigh  scattering,  Incoherent  Rayleigh 
scattering,  “three -photon”  scattering,  and  fluo¬ 
rescence,  respectively.  One  effect  of  the  colli¬ 
sions  is  to  increase  the  widths  of  the  Lorentziar 
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profiles.  However,  the  intensities  of  the  various 
components  are  also  affected.  The  largest  change 
occurs  in  the  fluorescence  which  can  now  occur 
both  by  the  three -photon  process  and  by  direct 
collisional  transfer  from  the  laser  excited  virtual 
level.  Thus  it  is  always  larger  than  the  three- 
photon  component,  as  shown  in  Fig.  3. 

We  can  now  see  from  Eq.  (18)  that  the  condition 
of  “weak  excitation”  can  be  precisely  stated  as 

»jn2«A*.  (22) 

When  the  laser  is  "weak"  and  the  collision  rates 
are  low,  the  spectrum  given  by  Eq.  (18)  reduces 
to  the  low  intensity  limit  given  by  Eq.  (9).  When 
02+  A2 »  y2,  the  components  can  be  resolved  and 
the  integrated  Rayleigh  intensity  lK  is  again  given 
by  Eq.  (16).  While  the  collisions  affect  the  coher¬ 
ent  and  incoherent  Rayleigh  intensities  separately, 
the  sum  f*  is  unaffected  because  Rayleigh  scat- 
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tering  can  occur  off  ground  or  excited  states;  thus 
collisions  do  not  affect  the  number  of  scattering 
centers.  The  integrated  side  bands  f3  and  /,  are 
given  by 

.  N(y,/4)n*(n'-A)[»,(a'-A)+Ai 

- WTaWT7W) - ’  l23) 

,  .  Af(y,/4)fl=(n'  +  A)[»,(fi'+A)-A] 

If - j&rmtFTWi - •  () 

When  rj  depends  on  the  detuning,  as  will  be  the 
case  outside  the  impact  regime,  the  theory  has 
not  been  worked  out.  Wc-  expect  that  the  integrated 
intensities  will  be  given  by  Eqs.  (23)  and  (24),  but 
that  the  line  shape  will  not  be  Lorentzian.  In  fact, 
we  expect  that  as  long  as  the  strong  laser  inten¬ 
sity  has  no  direct  influence  on  the  collision  rate, 
tj(A)  can  be  obtained  from  low -intensity  measure¬ 
ments  of  o£(A), 

It  should  be  noted  that  a  simple  theory,3*  using 
dressed  states  gives  all  the  results  for  integrated 
intensities  correctly  and  also  is  valid  in  the  tran¬ 
sient  regime. 

111.  EXPERIMENTAL  APPARATUS 

The  strontium  resonance  line  (5s5/>  ‘Pf-5s2  *S0) 
is  located  at  460.73  nm  (21698.482  cm-1).  The 
oscillator  strength  of  this  transition  is  /=1.94 
corresponding  to  a  radiative  lifetime  of  4.68 
nsec.86  Since  strontium  is  predominantly  an  even- 
even  nucleus  (93%),  there  is  no  hyperfine  struc¬ 
ture.  Therefore,  in  the  absence  of  collisions,  this 
system  will  respond  to  a  linearly  polarized  laser 
beam  as  a  two-level  system. 

A  nitrogen -laser -pumped  dye  laser  was  tuned 
near  this  460.73-nm  resonance  line  as  shown  in 
.•ig.  4  and  focused  into  an  oven  containing  stron- 
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FIG.  4.  Simplified  energy  levels  for  Sr.  The  laser, 
with  frequency  uL,  is  tuned  near  the  5s21Sn-  5s5p'Pl 
resonance  line  located  at  21698.482  cm'1  (460.73  nm). 


FIG.  5.  Schematic  diagram  of  experimental  apparatus. 

tium  vapor  and  argon  buffer  gas.  A  schematic  of 
the  experimental  apparatus  is  shown  in  Fig.  5. 

The  oven  was  generally  run  with  a  buffer  gas 
pressure  of  10  -500  Torr  and  in  a  temperature 
range  of  400-500  °C.  At  500  °C  the  vapor  pres¬ 
sure  of  strontium  is  ~8  x  10'3  Torr  and  the  Dop¬ 
pler  width  (half  width  at  e"  height)  of  the  460.73- 
nm  resonance  line  is  2.8xl0'=  cm-1  (5.9  xlO"1  nm). 
The  side  emission  from  the  excited  region  in  the 
oven  was  spectrally  analyzed.  Details  of  the  ap¬ 
paratus  are  discussed  below. 


A.  Dye  laser 

The  dye  laser  was  pumped  by  a  modified  version 
of  an  N2  laser  described  by  Woodward  etal 81  in 
the  standard  transverse  pumping  arrangement  with 
a  xlO  beam  expander  and  an  Echelle  grating  (316 
lines/mm).  The  N.  laser  output  was  focused  into  the 
dye  cell  with  a  simple  quartz  lens  (18  cm  focal  length). 
The  dye  cell  was  a  variation  of  a  basic  design  by 
Drullinger6"  and  consisted  of  a  stainless -steel 
block  fitted  with  a  quartz  entrance  window  and  two 
30'  wedged  side  windows  which  were  canted  at  5° 
to  prevent  multiple  reflections.  All  three  windows 
were  sealed  to  the  stainless  cell  with  Teflon  O 
rings.  The  dye  was  flowed  through  the  cell  with  a 
system  containing  parts  made  of  stainless  steel. 
Teflon,  and  polyethylene  and  including  a  Millipore 
filter.  The  output  of  the  dye  laser  was  almost 
totally  linearly  polarized  with  a  linewidth  of  1.4 
cm*1  (0.03  nm)  and  a  half -angle  beam  divergence 
of  1  mrad.  The  laser  had  a  5  nsec  pulse  length 
and  an  energy  of  ~0.2  mj  at  the  oven  entrance.  The 
dye  laser  beam  was  focused  into  the  oven  with  a 
x5  beam  contractor. 
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B.  Own 

The  vapor  cell  oven  was  made  from  2-in.-diam 
stainless -steel  tubing  in  the  shape  of  a  cross  with 
20-cm  legs  and  with  asbestos -insulated  nichrome 
heater  wire  wound  around  the  legs.  Each  leg  of  the 
cross  was  fitted  with  a  window.  Thermocouples 
and  a  temperature  controller  were  used  to  main¬ 
tain  the  temperature  to  within  ±2  °C.  The  stron¬ 
tium  was  held  in  a  stainless -steel  boat  with  5  mm 
holes  in  each  end  to  allow  passage  of  the  laser 
and  3  x  10  mm  slits  in  the  sides  to  allow  the  side 
emission  to  be  studied.  Since  we  were  looking  at 
Rayleigh  emission  with  the  same  frequency  as  the 
laser,  spurious  scatter  had  to  be  carefully 
avoided.  To  facilitate  this,  the  windows  through 
which  the  laser  passed  were  set  at  an  angle  to 
avoid  back  reflections.  To  protect  the  windows 
in  the  oven,  water  cooling  coils  were  attached  to 
the  stainless  cross  near  the  windows.  In  addition, 
a  buffer  gas  of  Ar  of  at  least  10  Terr  kept  the 
windows  from  being  coated  with  Sr. 

C.  Detection  system 

The  side  emission  was  studied  with  an  // 10  de¬ 
tection  system.  By  using  a  periscope  the  excita¬ 
tion  region,  which  was  approximately  a  cylinder 
of  110  pm  diameter,  was  imaged  parallel  onto  the 
slit  of  a  0.3-m  McPherson  model  218  monochro¬ 
mator  using  two  camera  lenses.  By  scanning  one 
of  the  lenses  in  the  detection  optics,  we  mea¬ 
sured  the  spatial  profile  of  the  Rayleigh  scatter¬ 
ing  at  low  intensities.  Using  this  profile,  we 
estimated  the  peak  power  of  the  laser  at  the  cen¬ 
ter  of  the  excitation  region  to  be  55  MW/cm2. 

The  resolution  of  the  monochromator  was  0.03  or 
0.06  nm  depending  on  the  grating  installed.  The 
output  of  the  RCA  7326,  S-20  photomultiplier, 
which  was  placed  at  the  exit  of  the  monochroma¬ 
tor,  was  fed  into  a  gated  (10  Msec)  analog  inte¬ 
grate -and -hold  circuit  and  recorded  on  a  chart 
recorder  with  a  1.5-sec  time  constant.  The  mono¬ 
chromator  measured  the  wavelength  to  a  reset¬ 
table  accuracy  of  0.01  nm. 

IV.  EXPERIMENTAL  STUDIES 

When  the  dye  laser  was  tuned  near  the  1Pf-lS0 
resonance  transition  at  460.73  nm,  three  compo¬ 
nents  were  observed  in  the  side  emission  as  pre¬ 
dicted  by  the  theoretical  discussion  of  Sec.  II. 
Figure  6  shows  the  spectral  distribution  of  the 
side  emission  for  two  different  laser  intensities 
with  the  laser  tuned  to  460.56  nm  (A=  -8.0  cm*1), 
the  blue  side  of  the  resonance.  The  central  com¬ 
ponent,  which  occurs  at  the  excitation  frequency 
is  the  Rayleigh  scattering.  The  width  of  this 
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FIG.  6.  Emission  spectra  for  the  laser  tuned  to  the  blue 
side  of  resonance  U,=-8.0  cm'1),  (a)  Full  laser  inten¬ 
sity,  (b)  laser  attenuated  by  a  factor  of  50.  The  factors 
give  the  detection  sensitivity.  Observed  are  the  Rayleigh 
scattering  at  lol,  the  resonance  fluorescence  at  u>0  and 
the  “three-photon”  scattering  at  w3. 

component  is  0.07  nm  which  is  determined  pri¬ 
marily  by  the  resolution  of  the  monochromator. 

The  component  near  u>}=  u>L-  A  is  due  to  three- 
photon  scattering.  Note  that  at  high  laser  intensity 
it  is  considerably  shifted  due  to  the  ac  Stark  ef¬ 
fect.  Since  our  laser  is  pulsed,  this  effect  shows 
up  both  as  a  shift  and  a  broadening  as  seen  pre¬ 
viously  by  others.'9’70  The  fluorescence  component 
is  located  near  the  resonance  frequency  u>0.  This 
component  is  a  result  of  both  three -photon  scat¬ 
tering  and  collisional  redistribution.  At  high  laser 
intensity  this  component  is  also  shifted  and  broad¬ 
ened.  Since  the  lifetime  of  the  excited  state  is 
comparable  to  the  laser  pulse  length,  only  part  of 
the  fluorescence  is  emitted  during  the  time  the 
laser  is  on.  Thus  the  shift  and  broadening  of  the 
fluorescence  due  to  the  ac  Stark  effect  are  not  as 
large  as  those  of  the  three -photon  component, 
which  is  emitted  only  during  the  laser  pulse. 

We  studied  these  emission  components  as  a  func¬ 
tion  of  argon  pressure,  laser  detuning,  and  laser 
intensity.  All  the  data  to  be  presented  will  be 
given  as  the  integrated  intensities  of  the  resolved 
spectral  components:  /„,  and  /,.  These  values 
were  obtained  by  graphically  integrating  scans 
such  as  those  in  Fig.  6  and  are  thus  limited  in  pre¬ 
cision  by  noise  in  the  apparatus.  We  estimate  the 
precision  of  these  measurements  to  be  i20‘«. 
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A.  Kmksion  vrnus  detuning  at  low  laser  intensities 

We  studied  the  effect  of  detuning  on  the  integrated 
Rayleigh  emission  lR  and  integrated  fluorescence 
/,  at  low  laser  intensities  over  a  large  range  of  A. 
The  oven  was  run  at  a  temperature  of  530  °C  and 
an  argon  buffer  gas  pressure  of  50  Torr.  At  this 
temperature  the  number  density  of  strontium  was 
estimated  to  be  9.7  xlO13  cm"  from  vapor  pres¬ 
sure  data.71  In  making  the  measurements,  the 
linearity  with  laser  intensity  was  checked  for  each 
point.  This  check  eliminated  the  possibility  of 
multiphoton  and  saturation  effects  and  guaranteed 
that  we  were  in  the  low  intensity  regime.  In  par¬ 
ticular  it  guaranteed  that  the  fluorescence  was  due 
only  to  collisional  transfer  from  the  virtual  level 
and  not  to  the  higher-order  three-photon  process. 

1.  Rayleigh  emission  versus  detuning 

Since  we  were  always  in  the  region  where  |  A| 

»y,  we  expected  from  Eq.  (10)  that  the  Rayleigh 
intensity  should  fall  off  as  A*2.  The  results  are 
shown  in  Fig.  7  for  detunings  on  both  sides  of 
resonance.  The  dashed  curve  shows  the  predicted 
A*2  dependence  adjusted  for  the  overall  gain  of 
the  detection  system.  We  see  that  for  two  decades 
of  detuning,  the  integrated  Rayleigh  intensity  goes 
as  A*2.  Note  that  A  was  never  large  enough  for 
the  quadrupole  allowed  'S0  -  lD2  transition  (cf. 

Fig.  4)  to  affect  the  Rayleigh  scattering.  In  Fig. 

7  scattering  from  the  surface  of  the  oven  (of  the 
order  of  IR  at  A=  200  cm")  was  subtracted.  We 
were  not  able  to  measure  the  Rayleigh  scattering 
at  A  »  +  5  cm*1  since  the  red  absorption  wing  of 
the  atoms  surrounding  the  excitation  region  at¬ 
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FIG.  7.  Rayleigh  scattering  (spectrally  integrated) 
versus  detuning  of  the  laser  from  resonance  at  low  laser 
intensity.  Detunings  to  the  blue  and  red  sides  of  reson¬ 
ance  are  shown  by  triangles  and  circles,  respectively. 
The  broken  line  shows  the  expected.-*2  Lorentzian  fall- 
off. 


tenuated  the  Rayleigh  intensity.  The  effect  of 
ground-state  atoms  surrounding  the  excitation  re¬ 
gion  and  the  resuiting  radiative  trapping  will  be 
more  fully  discussed  in  the  next  section  when  we 
consider  tnc  fluorescence  component. 

2.  Fluorescence  versus  detuning 

In  our  experiment,  the  impact  approximation  is 
expected  to  be  valid  only  for  |a|  <4.6  cm".  This 
was  determined  using  the  measurements  of  Penkin 
and  Shabanova73  and  our  Eq.  (6).  Thus  we  expect 
that  for  |a|  >4.6  cm'1,  the  impact  approximation 
will  fail  and  the  collision-induced  fluorescence 
will  no  longer  be  Lorentzian.  Experimentally  we 
find  that  indeed  the  fluorescence  is  strongly  fre¬ 
quency  dependent  and  asymmetric.  The  composite 
spectral  scans  of  Fig.  8  show  that  the  fluorescence 
component,  near  u>0  (slightly  displaced  by  the  ac 
Stark  effect)  is  an  order  of  magnitude  larger  when 
the  laser  is  tuned  to  the  red  side  of  the  resonance 
than  when  tuned  to  the  blue. 

Figure  9  shows  experimental  curves  for  the  in¬ 
tegrated  intensity  of  collision-induced  fluores¬ 
cence  as  a  function  of  A.  The  data  at  5 -cm*1  bor¬ 
der  the  impact  region.  The  units  of  / F  are  rela¬ 
tive  to  fR  of  Fig.  7.  The  vertical  error  bars  were 
determined  by  the  scatter  of  several  experimental 
readings  for  each  point,  while  the  horizontal  bars 
simply  show  the  1.4  cm"1  uncertainty  in  frequency 
due  to  the  laser  linewidth.  Also  drawn  is  a  broken 
curve  for  a  Lorentzian  fall -off  corresponding  to  a 
constant  yE  in  Eq.  (11),  where  yE  is  taken  from  the 
pressure-broadening  measurements.72  The  broken 
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FIG.  8.  Spectra  showing  asymmetry  of  fluorescence, 
at  u>p.  as  a  function  of  the  laser  frequency  (al  White 
light  absorption  profile  showing  location  of  resonance. 

(b)  and  (c)  show  the  laser  detuned  from  resonance  by 
plua  and  minus  15  cm",  respectively.  While  the  Ray¬ 
leigh  emission  is  symmetric  with  detuning,  the  collision- 
induced  fluorescence  is  highly  asymmetric. 
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FIG.  9.  Collision-induced  fluorescence  (spectrally  inte¬ 
grated)  as  a  function  of  detuning  of  the  laser  from  res¬ 
onance.  Dotunings  to  the  blue  and  red  sides  of  reso¬ 
nance  are  shown  by  triangles  and  circles,  respectively. 
The  broken  line  shows  thezi*1  dependence,  valid  in  the 
impact  region,  extended  to  Urge  detunings.  Impact 
approximation  is  only  valid  for  |A|  <4.6  cm"1. 

curve  was  corrected  (or  self -absorption  of  the 
sample  (see  Sec.  IV  A  3  on  radiative  trapping). 

It  is  seen  from  Fig.  9  that  yB  is  strongly  fre¬ 
quency  dependent,  as  discussed  in  Sec.  Q  A.  In 
particular,  when  the  laser  is  tuned  to  the  red  side 
of  resonance,  the  fluorescence  is  initially  much 
larger  than  when  the  laser  is  tuned  to  the  blue. 

This  “red  wing"  is  very  similar  to  the  results 
found  in  absorption  and  emission  experiments. 

This  similarity  is  due  to  the  fact  that  all  these 
experiments  effectively  measure  the  collisional 
line  shape  of  the  transition.  In  particular,  our 
experiment  can  be  viewed  as  an  absorption  event 
followed  by  an  emission  event 

The  origin  of  the  red  wing  can  be  understood 
from  very  basic  ideas  involving  potential  curves. 
Figure  10  shows  a  schematic  representation  of 
the  potential  curves  for  the  strontium -argon  col¬ 
lision.  Since  the  Sr(‘P)  level  is  more  polarizable 
than  the  Sr(*S),  the  potential  curve  for  the  Sr('P)- 
Ar('S)  collision  will  be  more  attractive  than  the 
potential  curve  for  the  Sr(*S)-Ar(‘S)  collision. 

When  the  laser  is  tuned  to  the  red  side  of  the  res¬ 
onance,  the  frequency  of  the  laser  will  be  equal  to 
the  difference  potential  at  some  intcrnuclear  sep¬ 
aration  and  there  will  be  a  large  overlap  between 
the  upper  and  lower  state  wave  functions.  Thus, 
according  to  the  Franck-Condon  principle,  a 
transition  will  be  induced.  This  side  of  resonance 
is  referred  to  as  the  quasistatic  wing.13  On  the 
other  hand,  when  the  laser  is  tuned  to  the  blue  side 
of  resonance,  the  overlap  is  less  so  that  the  col¬ 
lision-induced  fluorescence  is  smaller.  This  side 
of  resonance  is  referred  to  as  the  antistatic  wing. 11 
When  the  User  is  tuned  far  off  resonance  the  re¬ 
pulsive  parts  of  the  curves  will  become  important. 


FIG.  10.  Schematic  drawing  showing  the  general  be¬ 
havior  of  the  potentials  for  a  collision  between  the  stron¬ 
tium  i'P  or  3S)  atom  and  the  ground-state  argon  atom. 

For  instance  as  seen  from  the  schematic  of  Fig. 

10,  when  the  laser  is  tuned  far  off  resonance  to 
the  blue,  the  *2  upper  curve  and  the  ‘2  lower 
curve  run  parallel  so  that  the  difference  potential 
is  constant  and  the  transition  probability  increases. 
In  our  experiment  this  is  probably  the  reason  why 
the  blue  wing  becomes  more  significant  for  large 
detunings. 

In  order  to  bring  out  the  dependence  on  A  more 
clearly,  the  quantity  oE( A)  is  plotted  versus  A  in 
Fig.  11.  This  was  determined  by  the  ratio  of  Ir/ 


FIG.  11.  Cross  section  oh  tor  collision-induced 
fluorescence  as  a  function  of  detuning  at  930 'C.  De tun¬ 
ings  to  the  blue  and  red  sides  of  resonances  are  shown 
by  triangles  and  circles,  respectively.  The  dashed 
curves  show  the  prediction  of  the  theory  of  Szudy  and 
Ilaylis  (Ref.  13)  using  s  Van  dcr  Waats  interaction  con¬ 
stant  Ct=9.»  *10*51  cm'sec**:  ta)  the  impact  region, 
tbl  the  quasistalic  region,  and  tci  the  antistatic  region. 
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/  k  from  Figs.  7  and  9.  Account  was  also  taken  of 
the  attenuation  due  to  radiative  trapping  as  dis¬ 
cussed  in  the  next  section.  The  rate  of  collisional 
quenching  was  determined  self-consistently  from 
the  measurements  of  Sec.  IV  B  to  be  yt  =  0.29y„. 

We  have  compared  our  results  with  the  predic¬ 
tions  of  the  line -shape  theory  of  Szudy  and  Baylis.13 
Using  their  predictions  for  the  line  shape  from  a 
R~*  Van  der  Waals  potential,  we  obtained  the 
dashed  curves  for  the  (ai  impact  region,  (b)  quasi¬ 
static  region,  and  (c)  antistatic  region  when  we 
used  a  value  of  C„  =  9.8  xlO'31  cm6  see'1.  As  can 
be  seen  in  Fig.  11  the  value  of  aE  in  the  impact 
region  (|  A|  <  4.6  cm*1)  is  then  predicted  to  be 
4.4  xl0*u  cm2  for  530  “C.  This  is  higher  than  that 
obtained  at  530  °C  from  the  data  of  Penkin  and 
Shabanova72  (2.3  xlO'H  cm2)  and  also  Farr  and 
Hindmarsh71  (3.6  xlO*1'1  cm2).  The  differences  are 
within  the  expected  uncertainty. 

Under  our  experimental  conditions  we  found  that 
the  broadening  was  entirely  due  to  Sr-Ar  collisions 
and  not  Sr-Sr  collisions.  This  was  expected  in  the 
impact  region  from  measurements  of  collisional 
widths  for  Sr-Sr  collisions  and  Sr-Ar  collisions. 
The  values72'77  of  these  widths  at  530  °C  and  50 
Torr  of  argon  are  y£(Sr-Sr)=  5.6  xlO7  rad  sec'1 
and  ye(Sr-Ar)=  1.2  xlo,J  radsec*1.  This  expecta¬ 
tion  was  confirmed  by  repeating  the  experiment 
with  the  strontium  density  lowered  by  a  factor  of 
8  and  observing  no  change  in  the  shape  of  the  pro¬ 
files  in  Fig.  9.  This  is  consistent  also  with  our 
data  on  fluorescence  versus  argon  pressure  at  a 
detuning  of  9.4  cm'1  (see  Sec.  IV  B). 

For  large  detun1  tgs,  a  question  arose  as  to  the 
mechanism  responsible  for  populating  the  excited 
state.  Since  our  laser  had  a  spectral  width  of 
0.03  nm,  the  possibility  of  excitation  by  the  spec¬ 
tral  wings75'7’1  needed  to  be  considered,  in  addi¬ 
tion,  the  spectral  wings  of  the  laser  can  produce 
a  continuum  of  excited  states  between  tlwL  and 
K<j)0  which  could  be  collisionally  transferred  to 
the  excited  state.  When  the  laser  is  tuned  very 
far  from  resonance,  this  “cumulative  wing  effect” 
could  dominate  the  “normal”  collision -induced 
fluorescence  arising  from  the  virtual  state  at 

Direct  excitation  bv  spectral  components  at 
the  transition  frequency  was  expected  to  be  negli¬ 
gible  because  the  oven  had  an  optical  depth  of  -800, 
thus  absorbing  all  components  within  the  absorp¬ 
tion  line  before  reaching  the  excitation  region.  To 
verify  this  expectation  and  to  show  also  that  the 
“cumulative  wing  effect”  was  negligible  under  our 
conditions,  we  placed  a  0.03-nm  bandpass  mono¬ 
chromator  in  the  laser  beam  to  filter  out  any 
spectral  wings.  When  the  fluorescence  versus 
detuning  data  of  Fig.  9  were  taken  again  with  the 
filter  in,  no  differences  were  seen.  Thus  the  col¬ 


lision-induced  fluorescence  could  be  considered  to 
arise  entirely  from  collisional  transfer  from  the 
virtual  level  at  tiwL  to  the  excited  state,  where 
the  location  of  the  virtual  level  is  known  only  to 
within  the  line  width  of  the  laser. 

1.  Radiative  trapping 

The  small  cylindrical  region  of  Sr  atoms  that  is 
excited  by  the  laser  is  surrounded  by  unexcited  Sr 
atoms.  Most  of  our  experiments  were  done  at 
densities  where  this  unexcited  region  is  many  ab¬ 
sorption  lengths  deep  (-800).  The  simplest  way  to 
estimate  the  number  density  of  strontium  atoms 
and  therefore  the  optical  depth  is  by  measuring 
the  equivalent  width  of  the  absorption  line. M  This 
measurement  can  also  be  checked  with  vapor- 
pressure  curves.  As  a  result  of  the  radiative 
trapping,  the  fluorescence  component  of  the  scat¬ 
tered  intensity  was  heavily  attenuated.  The  loss 
of  intensity  comes  from  two  separate  sources: 
first,  radiative  trapping  lengthens  the  lifetime  of 
the  excitation  so  that  any  quenching  mechanism 
that  is  present  becomes  more  effective;  second, 
because  of  the  diffusion  of  the  radiation,  the  size 
of  the  radiative  region  increases,  therefore  that 
part  of  the  emission  that  is  collected  on  the  slit  of 
the  monochromator  decreases. 

In  order  to  verify  that  the  trapping  did  not  change 
the  shape  of  the  curves  in  Fig.  9,  but  only  the  over¬ 
all  gain,  we  checked  the  linearity  of  the  fluores  - 
cence  signal  in  laser  intensity  at  every  data  point 
by  inserting  a  neutral  density  filter  into  the  laser 
beam  and  observing  the  appropriate  linear  drop 
in  fluorescence  signal.  At  high  laser  intensities 
and  small  detunings,  the  ac  Stark  shift  by  the  laser 
can  be  large  enough  to  shift  the  emission  frequency 
outside  the  absorption  equivalent  width,  reducing 
the  amount  of  trapping  and  leading  to  an  increase 
in  fluorescence  signals.  The  linearity  check 
justified  our  using  a  higher  laser  intensity  to  pro¬ 
duce  signals  at  large  detunings  and  then  normaliz¬ 
ing  all  the  points  in  Fig.  9  to  some  fixed  laser  in¬ 
tensity.  That  the  fluorescence  was  linear  in  laser 
intensity  indicated  that  no  significant  Stark  effects 
were  occurring  for  these  measurements. 

It  was  mentioned  in  the  previous  section  that  the 
broken  line  in  Fig.  9  was  corrected  for  self -ab¬ 
sorption  of  the  sample.  This  was  done  by  multiply¬ 
ing  Eq.  (11)  by  a  constant  attenuation  factor  deter¬ 
mined  as  follows.  By  lowering  the  oven  tempera¬ 
ture  to  400  °C  (strontium  density  *5x10“  cm*3), 
we  substantially  reduced  the  amount  of  trapping. 

At  this  low  density,  however,  fluorescence  could 
be  observed  only  at  small  detunings.  When  the 
emission  was  observed  at  both  plus  and  minus 
5  cm*1  detuning,  which  borders  the  impact  region, 
the  fluorescence  relative  to  the  Rayleigh  scatter- 
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ing  was  within  a  factor  of  3  of  that  predicted  by 
Eqs.  (10)  and  (11),  indicating  an  absence  of  trap¬ 
ping.  Thus  we  used  the  ratio  of  the  relative  fluo¬ 
rescence  signals  with  and  without  trapping  as  our 
attenuation  factor. 

4.  Comparison  of  redlstrihntion  line  shape  with  absorption  and 
emission  experiments 

Absorption  experiments  measure  the  extinction 
coefficient  which  is  the  attenuation  of  a  well -col¬ 
limated  monochromatic  light  beam  at  frequency 
w£.  On  the  other  hand,  in  emission  experiments, 
the  atoms  are  excited  near  line  center  by  mono¬ 
chromatic  light,  and  the  fluorescence  spectrum 
is  measured  off  line  center.  Equation  (9)  shows 
that  in  the  impact  region  of  the  spectrum  the  line 
shapes  measured  in  absorption  and  emission 
should  be  the  same  as  the  profile  determined  by 
the  scattering  experiment  described  in  this  paper. 
Outside  the  impact  region,  the  theory  has  not  been 
explicitly  developed.  We  have  done  two  crude  ex¬ 
periments  to  determine  if  any  major  differences 
exist  between  these  various  experiments  in  the 
quasistatic  and  antistatic  region. 

Figure  12(a)  shows  a  scan  of  the  absorption  pro¬ 
file  of  the  460.73-nm  resonance  line  of  strontium. 
In  order  to  see  the  spectral  wings  we  operated  at 
a  temperature  of  580  °C  and  595  Torr  of  argon. 

The  scan  was  obtained  with  a  dc  broadband  light 
source  passing  through  the  oven  and  focused  onto 
the  monochromator  slit  As  in  our  redistribution 
experiments,  we  see  a  large  red  wing.  Figure 
12(b)  shows  a  scan  of  the  emission  profile  for  the 
same  temperature  and  pressure  conditions  when 


FIG.  12.  (a|  Scan  of  the  absorption  profile  of  the 
460.73-nm  line  in  strontium  at  SfiO’C  and  695- Torr 
argon,  (h)  Emission  profile  for  the  same  experimental 
conditions  with  Ihe  laser  tuned  lo  459.9  nm.  Doth  pro¬ 
files  show  a  large  red  wing  as  do  the  redistribution  re¬ 
sults  ( Fig.  91.  The  factor  gives  the  increase  in  detector 
sensitivity. 


the  laser  was  tuned  to  459.9  nm.  The  Rayleigh 
emission  is  seen  at  the  laser  frequency  as  ex¬ 
pected,  and  is  as  narrow  as  our  resolution.  As 
in  both  the  redistribution  experiment  (Fig.  9)  and 
the  absorption  experiment  [Fig.  12(a)],  the  fluo¬ 
rescence  line  shape  again  has  a  strong  red  wing. 
Note  that  for  this  higher  temperature  and  pres¬ 
sure  the  radiative  trapping  is  evident  from  the  dip 
in  the  fluorescence  near  oj0.  We  have  determined 
the  absorption  profile  [  -Iog(///0)]  from  Fig.  12(a) 
and  also  the  emission  profile  from  Fig.  12(b)  and 
find  that  the  shape  of  the  profiles  is  the  same  as 
the  redistribution  profile  shown  in  Fig.  9  to  within 
a  factor  of  2. 

B.  Emission  versus  argon  pressure 

We  have  studied  the  fluorescence  component  as 
a  function  of  argon  buffer  gas  pressure  with  the 
laser  tuned  to  460.93  nm  (A  =  9.4  cm'1).  The  study 
was  done  at  a  temperature  of  400  °C  so  that  effects 
of  radiative  trapping  were  minimal.  Agreement 
with  theory  was  obtained  when  the  inelastic  rate 
y,  was  taken  to  be  linear  in  argon  pressure.  Ex¬ 
perimental  and  theoretical  curves  showing  the 
fluorescence  versus  Ar  pressure  at  three  differ¬ 
ent  laser  intensities  are  shown  in  Fig.  13. 

The  theoretical  curves  were  computed  from  Eq. 
(24)  with  two  adjustable  parameters:  y,  and  O' 
(laser  intensity  for  curve  a).  The  laser  intensities 
for  the  lower  curves  were  determined  by  the 
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FIG.  13.  Collision- Induced  fluorescence  Ispcctrally 
integrated!  as  a  function  of  argon  buffer  gas  pressure  tor 
three  different  laser  intensities,  (a)  Full  laser  intensity, 
tbl  laser  attenuated  by  4.5.  ic)  laser  attenuated  by  20. 

The  solid  curves  were  computed  from  Eq.  (24)  with  two 
common  parameters:  the  quenching  rate  if  and  the 
laser  intensity  for  curve  (al. 
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•mown  attenuations  of  the  neutral  density  filters 
used  in  the  experiment.  The  value  of  yt  deter¬ 
mined  by  this  fitting  technique  is  yt/yt  -  0.007 
1 0.002.  The  value  of  fls  was  determined  to  be  flV 
A’ =  0.1  ±0.03  (0.5  MW/cm-).  The  radiative  decay 
rate  y„  was  taken"'1  to  be  1.07  xlO8  rad  sec*1  and 
the  elastic  collision  race  yE  was  determined  self  - 
consistently  from  the  data  in  Fig.  11. 

The  general  form  of  the  curves  for  low  laser  in¬ 
tensity  can  be  understood  timili/ntirely  by  consider¬ 
ing  the  rate  scheme  shown  in  Fig.  1.  Following 
Eq.  (17)  it  was  pointed  out  that  the  “population"  in 
the  upper  virtual  level  can  be  taken  as  (N/ 2)(02/ 
A1).  This  “population”  is  then  transferred  to  the 
upper  state  at  a  rate  yE  due  to  elastic  collisions. 
Thus  the  transfer  rate  from  the  lower  state  to  the 
upper  state  is  yr  =  (y£/2)(fl2/A2).  There  is  also  a 
route  for  laser-induced  collisional  depopulation  of 
the  excited  state.  Thus  we  expect  that  the  curves 
should  be  linear  at  low  Ar  pressure  with  the  slopes 
increasing  with  increasing  laser  intensity  and  that 
each  curve  should  saturate  to  some  maximum 
fluorescence  level,  due  tc  equalization  of  the  ex¬ 
cited  and  ground  state  populations.  In  the  absence 
of  quenching,  all  curves  would  saturate  to  the 
same  level.  This  is  because  the  transfer  rate  yT 
will  eventually  dominate  the  radiative  decay  rate 
y„.  However,  in  the  presence  of  quenching,  there 
is  a  competition  between  the  quenching  rate  y ,  and 
the  transfer  rate  yT  between  levels  so  that  each 
curve  saturates  to  a  different  intensity  as  seen  in 
Fig.  13.  In  fact,  if  the  fluorescence  is  taken  to  be 
given  by  the  decay  rate  y„  times  a  steady -state 
population,  we  obtain  Ir  =  NyK  [  y r/(y/  +  y*)].  It  is 
interesting  to  note  that  Eq.  (24),  which  was  used 
to  compute  the  curves  in  Fig.  11,  reduces  to  this 
expression  at  low  laser  intensities. 


C.  Emission  versus  laser  intensity 

We  studied  the  saturation  of  the  three  emission 
components  at  high  laser  intensities.  The  oven 
was  run  at  490 °C  with  an  argon  buffer  gas  pres¬ 
sure  of  10  Torr.  Thus,  the  number  density  of 
argon  was  Af(Ar)  =  1.3  xl0lT  cm*3  and  the  number 
density  of  strontium  was  N(St)=  2.5  xlO13  cm*3. 

In  our  Initial  comparisons  with  the  theoretical 
predictions  given  by  Eqs.  (16),  (23),  and  (24),  we 
assumed  the  laser  intensity  was  constant  spatially. 
However,  we  found  that  at  high  degrees  of  satura¬ 
tion,  there  were  deviations  from  the  theory  due  to 
the  nonlinear  effects  of  the  spatial  wings  of  the 
laser.  Therefore,  we  decided  to  average  the  theo¬ 
retical  predictions  over  the  spatial  profile  of  the 
laser.  The  method  used  for  the  spatial  averaging 
is  discussed  below. 


/.  Spatial  arerage 

In  the  coordinate  system  shown  in  Fig.  14,  the 
laser  beam  travels  along  the  z  direction  (out  of 
the  page)  and  the  detection  system  is  along  the  x 
direction.  For  an  emission  profile  t(r),  where 
r  is  the  radial  distance,  the  observed  intensity 
profile  is 

/(v)  =  2  j\(r)  dx  =  2  j~  ■  (25) 

By  using  the  Abel  inversion,"'  «(»•)  can  be  obtained: 

.(,1-1  £  (2.) 

To  determine  the  laser  profile  lL(r)  in  the  interac¬ 
tion  region,  we  scanned  the  collection  lens  in  the 
y  direction  while  the  monochromator  observed 
the  Rayleigh  emission.  The  laser  was  run  at  low 
intensities  so  the  Rayleigh  emission  was  linear 
in  laser  intensity.  Then  by  using  Eq.  (26)  we  found 
that  the  laser  profile  lL(r )  (photons  cm*2 sec*1) 
could  be  approximated  by 

/  r2  \ 212 

(27) 

over  the  range  of  rs  700  pm.  Outside  this  range 
our  signal  was  in  the  noise.  Here  lL  =IL(r  =  0), 
and  r0  was  found  to  be  110  pm.  The  spatial  pro¬ 
file  of  the  laser  beam  was  then  folded  into  Eqs. 
(16),  (23),  and  (24)  using  Eq.  (25).  In  this  way, 
theoretical  predictions  for  the  three  emission 
components  as  a  function  of  laser  intensity  lL 
were  obtained,  as  presented  in  the  following  sec¬ 
tions. 


FIG.  14.  Coordinate  system  for  spatial  averaging. 

The  laser  beam  travels  along  the  s  axis  (out  of  the  page) 
and  the  detection  system  is  in  the  x  direction. 
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2.  Ray  kith  component 

The  dependence  of  the  spectral  integrated  Ray¬ 
leigh  intensity  on  the  incident  laser  intensity  is 
given  by  Eq.  (16).  This  equation  shows  that  at  low 
laser  intensities,  the  dependence  is  linear  in  laser 
intensity.  However,  when  0:*  A2,  saturation  be¬ 
gins  to  occur  and  the  Rayleigh  intensity  falls  be¬ 
low  the  linear  dependence,  eventually  becoming 
constant. 

We  studied  the  saturation  effect  of  the  Rayleigh 
scattering  at  several  detunings  of  the  laser  on 
both  sides  of  resonance.  The  results  are  sh  •«-,n 
in  Fig.  15  where  it  is  seen  that  the  saturation  is 
less  rapid  when  the  laser  is  tuned  far  from  reson¬ 
ance  than  when  it  is  tuned  near  to  resonance.  The 
theoretical  curves  shown  by  the  broken  lines  were 
obtained  assuming  a  rectangular  spatial  profile 
for  the  laser  and  are  therefore  plots  of  Eq.  (16) 
with  two  free  parameters— the  incident  laser  in¬ 
tensity  and  the  overall  detection  sensitivity.  All 
three  curves  use  the  same  values  for  these  param¬ 
eters;  only  A  was  changed  as  measured  experi¬ 
mentally.  The  fit  was  obtained  with  an  incident 
laser  intensity  of  35  MW/cm2.  When  a  spatial 
average  was  taken  as  described  in  Sec.  IV  C 1,  the 
theore  tical  curves  shown  by  the  solid  lines  were 
obtained.  The  fit  was  obtained  with  a  laser  inten¬ 
sity  of  105  MW/cm2  at  the  center  of  the  excitation 
region.  The  deviation  of  the  fit  with  and  without 
spatial  averaging  is  comparable  to  our  experimen- 


0  05  1.0 

NORMALIZED  LASER  INTENSITY.  IL/l£ 


FIG.  15.  Rayleigh  scattering  (spectrally  Integrated) 
versus  laser  Intensity  for  several  detunings.  Detunings 
to  the  blue  and  red  sides  of  resonance  are  shown  by  tri¬ 
angles  and  circles,  respectively,  (a)  A  =  ±0.17  ntn,  (b) 

A  =  i0.36nm,  (c)  A  =  ±0.75  nm.  Near  resonance,  satura¬ 
tion  is  seen.  The  theoretical  curves  were  obtained  from 
Eq.  (16)  assuming  a  rectangular  spatial  profile  (broken 
curves),  and  using  the  measured  spatial  profile  given  by 
Eq.  (27)  (solid  curves).  See  Sec.  IV C2  for  discussion. 


FIG.  16.  Three-photon  scattering  (spectrally  inte¬ 
grated)  versus  laser  intensity  for  a  detuning  of  0.17  nm 
on  the  red  side  (a)  and  the  blue  side  (b)  of  resonance. 
Solid  curves  are  theoretical  predictions  using  Eq.  (23) 
with  no  free  parameters.  The  insert  shows  the  region 
at  low  intensity  expanded  by  five.  The  intensity  of  the 
three-photon  component  is  predicted  to  be  quadratic  in 
the  region  below  a  normalized  laser  intensity  of  0.01. 

See  Sec.  IVC3  for  further  discussion. 

tal  error.  Note  that  for  |a|  =0.75  nm,  spatial 
averaging  does  not  change  the  theoretical  curve. 

3.  Three-photon  component 

Theory  predicts  that  at  low  intensities  the  three- 
photon  component  is  quadratic  in  laser  intensity 
but  eventually  saturates  and  levels  off  as  does  the 
Rayleigh  intensity.  As  can  be  seen  from  Eq.  (23), 
the  integrated  intensity  J3  of  the  three -photon  com¬ 
ponent  not  only  depends  on  the  detuning  and  inten¬ 
sity  of  the  laser,  but  also  to  a  small  degree  on  the 
collisional  rates. 

We  have  observed  this  three -photon  component 
at  the  predicted  position  for  both  positive  and 
negative  values  of  the  detuning  A  and  for  various 
laser  intensities.  In  particular  at  A  =  ±0.17  nm 
we  have  measured  the  intensity  dependence  of  this 
component.  The  results  are  shown  in  Fig.  16.  . 
The  vertical  scale  is  the  same  relative  scale 
used  in  plotting  the  Rayleigh  intensity  in  Fig.  15. 
The  curves  shown  in  Fig.  16  are  zero  parameter 
predictions  for  the  integrated  intensity  /,  of  the 
three-photon  component  using  Eq.  (23)  and  the 
values  for  the  fixed  parameters  determined  as 
follows.  The  overall  sensitivity  of  the  detection 
system,  spatial  average  and  incident  laser  inten¬ 
sity  are  the  same  as  used  *o  fit  the  Rayleigh  sat¬ 
uration  data  in  Fig.  15.  Also  y,  was  determined 
from  our  data  presented  in  Sec.  IV  B.  Finally, 
yt  was  determined  from  our  low  intensity  data  on 
the  fluorescence  (Fig.  11).  We  see  that  this  zero 
parameter  fit  is  within  our  20'ib  experimental  ac¬ 
curacy. 
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4.  I  lmirrui'iu  v  comiiuuent 

We  studied  the  laser  intensity  dependence  of  the 
spectrally  integrated  fluorescence  at  various  de¬ 
tunings  on  both  sides  of  resonance.  Data  at  A  = 
±0.17  nm  are  shown  in  Fig.  17.  The  relative  units 
of  the  vertical  scale  of  Fig.  17  arc  the  same  as  in 
Fig.  15  for  the  Rayleigh  intensity. 

The  solid  curves  in  Fig.  17  are  theoretical  pre¬ 
dictions  using  Eq.  (24)  with  the  spatial  averaging 
procedure  discussed  in  Sec.  rVCl.  The  fits  were 
obtained  using  the  elastic  collision  rate  oe  as  a 
free  variable.  All  other  parameters  were  de¬ 
termined  by  previous  data  as  for  I3.  The  best 
fit  to  the  data  on  the  red  side  of  resonance  was 
obtained  with  cr£  -  320  x  10'16  cm2.  This  is  about 
a  factor  of  3  smaller  than  expected  from  our  low- 
intensity  measurements  of  the  fluorescence  shown 
in  Fig.  11.  On  the  blue  side  of  resonance  the  best 
fit  was  for  o£=  18  x  10*16  cm2  which  is  a  factor  of 
8  smaller  than  our  low-intensity  measurement. 

The  reason  for  this  discrepancy  probably  lies 
in  the  fact  that  we  have  applied  a  steady -state  the¬ 
ory  to  a  pulsed  experiment.  Unlike  the  Rayleigh 
and  three -photon  components,  which  are  only 
emitted  during  the  laser  pulse,  the  fluorescence 


FIG.  17.  Resonance  fluorescence  (spectrally  integrat¬ 
ed)  versus  laser  intensity  for  a  detuning  of  0.17  nm  on 
the  red  side  (a)  and  the  blue  side  (b)  of  resonance.  Satu¬ 
ration  on  the  red  side  is  very  pronounced.  The  theore¬ 
tical  curves  were  obtained  from  Eq.  <24>  assuming  a 
rectangular  profile  (broken  curve)  and  using  the  mea¬ 
sured  spatial  profile  given  by  Eq.  (27)  (solid  curve). 

See  8ec.  1VC  4  for  discussion. 


decays  exponentially  after  the  laser  has  turned  off. 
As  discussed  in  Sec.  IV  A  3  there  is  considerable 
radiation  trapping  due  to  the  surrounding  stron¬ 
tium  atoms.  At  high  laser  intensities  and  small 
detunings,  the  ac  Stark  effect  is  large  enough  that 
the  fluorescence  emitted  during  the  laser  pulse 
is  essentially  untrapped.  However,  the  fluo¬ 
rescence  emitted  after  the  laser  pulse  is  heavily 
trapped,  so  that  the  total  fluorescence  is  less  than 
expected,  leading  to  a  smaller  value  of  o£  than 
expected.  We  are  presently  attempting  to  extend 
our  analysis  to  a  transient  formulation  using  the 
theory  of  Courtens  and  SzBke.-" 

We  also  have  calculated  the  theoretical  pre¬ 
dictions  for  the  fluorescence  component  without 
spatial  averaging  and  using  an  incident  laser 
power  of  35  MW/cm2  determined  by  the  fit  to  the 
Rayleigh  saturation  data  presented  in  Fig.  15.  The 
resultant  theoretical  curves  are  shown  in  Fig.  17 
by  the  broken  lines.  All  other  variables  are  the 
same  as  for  the  theory  shown  by  the  solid  lines. 
Without  spatial  averaging,  the  theory  predicts 
that  the  fluorescence  at  A  =  +  0. 17  nm  will  saturate 
heavily.  However,  when  we  take  the  spatial  aver¬ 
age,  the  effect  of  extreme  saturation  is  washed 
out  by  the  spatial  wings  as  seen  by  the  solid  curve. 
The  reason  that  the  effect  of  the  wings  is  so  much 
more  important  in  the  fluorescence  than  in  the 
Rayleigh  emission  is  that  the  slope  of  the  fluor¬ 
escence  at  low  intensities  is  larger  so  that  the 
wings  are  heavily  weighted  in  the  average. 

We  have  also  studied  the  fluorescence  as  a  func¬ 
tion  of  the  intensity  of  the  laser  at  larger  de- 
tunings,  where  there  is  less  saturation  and  the 
dependence  is  more  linear.  Good  theoretical  fits 
could  be  obtained  using  Eq.  (24),  but  again  all 
values  of  ofi(A)  determined  from  these  fits  were 
less  than  the  values  expected  from  the  low  inten¬ 
sity  data  shown  in  Fig.  11.  As  discussed  above, 
we  feel  these  differences  are  due  to  the  steady  - 
state  assumption  coupled  with  the  effects  of  radia¬ 
tive  trapping. 

V.  SUMMARY  AND  DISCUSSION 

Our  studies1*2  were  the  first  ones  that  actually 
observed  the  collisional  redistribution  of  scattered 
radiation.  The  knowledge  of  this  function  is  im¬ 
portant  to  the  theory  of  radiative  transport  and 
line  formation  in  stellar  atmospheres.  It  has  been 
calculated  from  first  principles  in  the  impact 
region,  but  no  firm  theory  exists  in  the  static  re¬ 
gion.  Our  experiments  for  low  incident  inten¬ 
sities  can  be  summarized  in  terms  of  an  em¬ 
pirical  redistribution  function  that  is  similar  to 
Eq.  (9): 
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x  (o(a)  -  (o£)  •  (28) 

where  yfi(A)  can  be  obtained  from  Fig.  11.  Here 
/(io)  is  the  emission  line  shape,  normalized  to 
/ /(o))db)  =  1.  We  actually  measured  only  pieces 
of  this  formula,  i.e.,  /fay£(A)/4!,  ac  l/A2, 
and  in  a  separate  experiment  /(oj),  where  we  also 
found  /(A)  ccy£(A)/A2  within  our  experimental  range 
and  error.  We  have  also  shown  that  experiments 
with  variable  perturber  density  can  establish  the 
inelastic  collision  rate  y„  though  we  have  done  no 
experiments  to  determine  the  actual  cause  of 
quenching. 

Radiation  trapping,  one  of  the  major  obstacles 
in  our  experiments,  is  an  interesting  subject  by 
itself.  We  conducted  enough  experiments  to  con¬ 
vince  ourselves  that  we  can  account  for  the  in¬ 
fluence  of  trapping  consistently.  It  is  interesting 
to  note  that  experiments  similar  to  ours  can  be 
conducted  to  study  details  of  the  space -time  be¬ 
havior  of  radiative  trapping.82  By  measuring  the 
space-,  time-,  and  spectrally -resolved  scattered 
light  after  a  small  cylindrical  region  is  excited  by 
a  pulsed  laser,  the  most  fundamental  quantity  in 
radiative  transfer,  the  Green’s  function,  can  be 
measured  experimentally. 

The  high-field  behavior  of  the  Rayleigh  and 
three -photon  scattering  has  been  studied  and 
agreement  was  obtained  with  Mollow’s  steady -state 
theory.  While  the  integrated  Rayleigh  component 


was  found  to  be  “collisionless,  ”  the  third  compon¬ 
ent  was  found  to  have  a  weak  collisional  depen¬ 
dence.  The  fluorescence  component  was  also 
found  to  fit  the  steady -state  theory  quite  well,  but 
the  cross  section  for  the  elastic  collision  was 
smaller  than  that  predicted  by  our  low  intensity 
measurements.  We  feel  that  this  discrepancy  is 
due  in  part  to  the  transient  nature  of  the  laser  ex¬ 
citation  because  (a)  there  is  not  enough  time  for 
collisional  equilibration,  and  (b)  during  the  strong 
laser  pulse  the  fluorescence  is  ac  Stark  shifted 
and  therefore  escapes,  but  after  the  pulse  is  over 
the  rest  of  the  fluorescence  is  trapped  and  partly 
lost  to  the  detection  system.  It  is  interesting  to 
note  that  according  to  recent  studies30*32  it  is  ex¬ 
pected  that  the  elastic  collision  rate  should  have 
some  intensity  dependence.  We  continue  our 
studies  in  this  direction. 

Note  added  in  proof.  Recently  Nienhuis  and 
Schuller  (private  communication)  have  extended 
the  theory  of  collisional  redistribution  into  the 
far  wing  where  the  impact  limit  is  invalid.  They 
find  that  their  predictions  are  in  qualitative  agree¬ 
ment  with  our  observations  of  the  asymmetry  in 
the  collisional  redistribution. 
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We  have  observed  collisional  redistribution  of  laser  light  scattered  by  a  three-level 
atom,  thallium,  in  the  presence  of  argon.  The  two  observed  spectral  components  corre¬ 
spond  to  Raman  scattering  and  collision- Induced  fluorescence.  We  have  seen  the  growth 
of  both  components  from  spontaneous  to  stimulated  scattering  as  the  laser  intensity  is  in¬ 
creased.  While  the  gain  for  the  Raman  component  agrees  with  theory,  the  observed  gain 
of  the  collision-induced  fluorescence  Is  over  an  order  of  magnitude  lower  than  predicted. 


Since  the  experimental  study  by  Rousseau,  Pat¬ 
terson,  and  Williams,1  there  has  been  consider¬ 
able  theoretical  and  experimental  interest  in  the 
collisional  redistribution  of  near-resonant  scat¬ 
tered  light.2*6  For  a  two-level  atom  in  the  pres¬ 
ence  of  collisions,  the  spectrum  of  the  scattered 
light  at  low  intensities  consists  of  two  compon¬ 
ents;  elastic  scattering  at  the  frequency  of  the 
incident  light,  referred  to  as  Rayleigh  scatter¬ 
ing,  and  inelastic  scattering  at  the  frequency  sep¬ 
aration  of  the  two  levels,  which  is  called  colli¬ 
sion-induced  fluorescence.  For  a  folded  three- 
level  system,  as  shown  in  Fig.  1,  there  will  also 
be  two  emission  components  near  the  2-3  fre¬ 
quency:  Raman  Stokes  scattering  at  w,  as  well 
as  collision-induced  fluorescence  at  u»„.  Mol- 
low7  has  predicted  that  both  of  these  latter  com- 


FIG.  1.  (a)  First  three  energy  levels  of  thallium. 

With  the  laser  tuned  near  the  377.6-nm  resonance 
line,  we  observed  scattered  light  near  the  535.0-nm 
line,  (b)  Schematic  representation  of  electron  Raman 
scattering  at  the  Stokes  frequency  w,  from  an  fncldeot 
laser  at  frequency  (c)  Schematic  representation 
of  collision- Induced  fluorescence.  Collisions  (which 
In  our  case  were  Tl-Ar  collisions)  transferred  T1 
atoms  from  the  laser-induced  virtual  level  to  level  2, 
resulting  in  fluorescence  at  When  level  3  Is  Initial¬ 
ly  unpopulated,  both  of  these  components  can  become 
stimulated. 


ponents  can  have  gain.  Therefore  the  possibility 
exists  for  simultaneous  generation  of  stimulated 
emission  at  both  frequencies,  w,  and  u23.  In  an 
experiment  by  Wynne  and  Sorokin,*  there  are  in¬ 
dications  that  these  two  processes  were  occur¬ 
ring  but  the  two  components  could  not  oe -well  re¬ 
solved  and  the  population  mechanism  for  the  flu¬ 
orescence  was  not  determined. 

This  Letter  reports  on  the  spectrally  resolved 
observation  of  the  simultaneous  generation  of 
stimulated  Raman  scattering  (SRS)  and  stimulated 
collision-induced  fluorescence  (SCF).  Using  a 
dye  laser  tuned  near  the  377.6-nm  l61Pl/g-7iSl/g) 
resonance  line  of  Tl,  we  observed  the  growth  of 
both  of  these  components  (near  the  535.  O-nm 
emission  line)  from  the  linear  (with  pump  laser) 
regime,  where  the  scattering  is  spontaneous,  to 
the  exponential  regime,  where  the  scattering  be¬ 
comes  stimulated.  In  addition,  we  have  studied 
the  collisional  dependence  of  the  SCF  and  SRS. 
Recently  SRS  in  vapors  and  gases  has  been  used 
by  a  number  of  researchers  as  an  efficient  means 
of  down-conversion.®  This  study  is  an  effort  to 
understand  the  effects  of  collisions  on  such  stim¬ 
ulated  scattering. 

For  a  low-intensity  monochromatic  laser  at 
frequency  w£,  detuned  far  from  resonance,  the 
spontaneous  Raman  scattering  at  u>s  has  a  steady- 
state  intensity  (in  photons  cm'5  sec*1)  given  by10 

U*iNr”<ia/ A)*  (1) 

where  N  is  the  number  density  of  scattering 
atoms  (in  our  case  Tl  atoms),  r,23  is  the  spon¬ 
taneous  decay  rate  from  level  2  to  level  3,  A. 

=  «*,  -  <dl  is  the  detuning,  Q  *  txE/fl  is  the  Rabi 
frequency  associated  with  the  incident  laser  field 
E,  and  m  is  the  dipole  matrix  element  between 
states  1  and  2.  Similarly,  the  collision-induced 
fluorescence  at  has  a  steady-state  intensity 
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given  by10 


/,  =  |Ny,(A)(n/A)*yJ,**/yjw  (2) 


where  >a(A)  is  the  rate  of  quasielastic  collisions 
(in  our  case  Tl-Ar  collisions)  which  make  up  the 
energy  difference  A  needed  to  produce  an  atom  in 
level  2  and  Ym^Ym^+Ym*3-  We  have  written  the 
collision  rate  as  a  function  of  the  detuning  to  ac¬ 
count  for  the  non-Lorentzian  dependence  outside 
the  impact  regime.11  Effects  of  spatial  degener¬ 
acy  have  not  been  included  in  Eqs.  (1)  and  (2). 

Mollow7  has  predicted  that  both  fR  and  I,  will 
have  gain  when  level  3  is  initially  unpopulated. 
For  laser  intensities  where  ft2  <  A2,  we  expect 
the  gain  of  the  Raman  scattering  at  w,  to  be 


c2  >^22n2 

>  R  A1’ 


(3) 


where  y„  is  a  measure  of  the  Raman  linewidth 
(to  be  discussed  later).  For  the  collision-in¬ 
duced  fluorescence  at  the  gain  is  expected 
to  be 


ksCF  ■ 


£  w  C*  r.(A)  O 
2  v„2[y[/nl/\A 


t 

f  • 


(4) 


Here  y0  is  the  Doppler  linewidth  and  U  is  the 
peak  height  of  the  normalized  Voight  profile,12 
both  for  the  2-3  transition.  One  can  think  of 
Yd/vU  as  an  effective  width  for  the  Doppler-plus- 
colli sion -broadened  transition,  2-3. 

Expressions  (1)  and  (2)  for  the  spontaneous 
emission  and  Eqs.  (3)  and  (4)  for  the  gain  can  be 
used  with  simple  photon  propagation  equations12 
to  solve  for  the  single-pass  stimulated  outputs 
(in  photons  sec'1)  of  the  SRSand  SCF  collected 
by  a  solid  angle  a,: 

f  srs  =  (fRAo,/4jr^SRS)[exp(gSRSL)-  1] 

+/RA(a,-a,)L/4*  (5) 

and 

f  SCF  =  (/f  AO,/4j7^'sl;,)[eXp(|l,'sf;pL)  —  1] 

+  IFA(aa-  ar,)L/4i,  (6) 

where  A  and  L  are  the  area  and  length  of  the  ex¬ 
citation  region,  a,  is  the  solid  angle  formed  by 
the  excitation  region,  and  it  is  assumed  that  a, 

•  In  both  equations,  the  second  term  ac¬ 
counts  for  the  spontaneous  emission  which  exits 
from  the  sample  outside  of  the  gain  region,  but 
which  still  enters  our  collection  angle  or2.  We 
note  that  Eqs.  (5)  and  (6)  predict  that  both  the 
Sits  and  SCF  will  initially  have  a  linear  depen¬ 


dence  on  laser  intensity  when  spontaneous  scat¬ 
tering  is  dominant,  but  will  eventually  grow  ex¬ 
ponentially  when  the  scattering  becomes  stimu¬ 
lated.  Effects  of  saturation  or  population  de¬ 
pletion  have  not  been  included. 

The  apparatus  used  in  this  experiment  is  simi¬ 
lar  to  that  described  in  detail  by  Carlsten,  Sztike, 
and  Raymer.2  A  tunable  dye  laser,  pumped  by 
an  Nj  laser,  was  tuned  near  the  62Ju1/1-72S1/l 
resonance  line  of  T1  at  377.6  nm.  The  dye  laser 
had  a  pulse  duration  of  10  nsec,  a  spectral  width 
of  0.03  nm,  and  an  energy  of  50  pj  inside  the  ex¬ 
citation  region.  The  beam  was  focused  to  330  pm 
diam  giving  a  power  density  of  2  MW  /cm2  inside 
an  oven  containing  0.1  Torr  of  T1  vapor  and  5  to  80 
Torr  of  Ar  buffer  gas.  The  oven  ihput  and  output 
windows  were  put  at  an  angle  to  avoid  back- re¬ 
flections,  which  would  affect  the  growth  of  *he 
stimulated  emission.  The  length  of  the  T1  vapor 
region  was  -2.5  cm.  The  emission  region  was 
then  imaged  onto  the  slit  of  a  monochromator  of 
0.06  nm  resolution  with  an  / /10  optical  system 
that  was  capable  of  viewing  the  emission  at  right 
angles  to  the  laser  beam  or  along  the  laser  beam 
direction. 

When  the  incident  laser  is  tuned  on  resonance, 
only  one  spectral  component,  centered  at  u;2„  is 
observed  in  the  emission.  We  have  studied  this 
emission  at  535  nm  in  both  the  side  and  forward 
.  directions  when  the  laser  was  tuned  to  the  377.6- 
nm  resonance  line  of  Tl.  At  low  laser  intensities 
we  observed  that  both  the  side  and  forward  emis¬ 
sions  were  linear  in  laser  intensity.  Above  - 10 
kW/cm2,  the  forward  emission  became  stimulat¬ 
ed  and  eventually  saturated,  allowing  a  maximum 
photon  conversion  efficiency  of  60%. 

In  order  to  study  the  collisional  effects  upon 
this  stimulated  scattering,  we  tuned  the  laser 
0.14  nm  to  the  red  side  of  resonance.  We  were 
then  able  to  resolve  spectrally  the  Raman  emis¬ 
sion  at  o>s  and  the  collision-induced  fluorescence 
at  u.2s.  The  dependence  of  these  two  resolved 
components  on  laser  intensity  is  shown  in  Fig. 

2.  The  Ar  pressure  for  these  data  was  20  Torr. 

In  the  side  direction  [Fig.  2(b)J  both  the  Raman 
scattering  / R  and  collision-induced  fluorescence 
Ir  were  linear  in  laser  intensity.  By  measuring 
the  ratio  lF/l*t  we  obtained  an  absolute  measure 
of  the  collisional  redistribution  function.  From 
Eqs.  (1)  and  (2),  this  ratio  is  expected  to  be 

(7) 

Using14 )  „  =  6.7  a  10T  rad  sec*1  and  our  measured 
value  for  /,//„  we  obtain  >,(0.14  nm  to  red) *7.4 
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normalized  laser  intensity  . 


FIG.  2.  Dependence  of  Raman  scattering  and  colli¬ 
sion-induced  fluorescence  on  laser  Intensity  when  the 
laser  was  tuned  off  resonance  (0.14  nm  to  the  red  side 
of  the  377.6- nm  resonance  line),  (a)  Growth  of  the 


diam  of  the  excitation  region.  It  is  important  to 
note  that  in  calculating  the  SRS  gain  we  have  used 
the  result  of  Akhmanov,  D'yakov,  and  Pavlov1* 
that,  for  forward  SRS  in  a  dispcrsionless  medium, 
a  broad-band  laser  has  the  same  gain  as  a  mono¬ 
chromatic  laser.  This  is  because  the  intehsity 
variations  of  the  Stokes  emission  follow  those  of 
the  pump  laser  as  the  two  pulses  travel  with  the 
Bame  velocity  through  the  medium.  Thus  for  >  R 
we  have  used  the  Raman  linewidth  (which  is  pre¬ 
dominantly  the  Doppler  width  for  the  1-3  separa¬ 
tion).  This  results  in  an  SRS  gain  which  is  240 
times  larger  than  that  calculated  by  taking  yR 
equal  to  the  laser  linewidth.  We  see  that,  while 
the  predicted  initial  exponential  growth  of  the 
Raman  scattering  is  quite  close  to  our  experi¬ 
mental  results,  the  gain  for  the  collision-induced 
fluorescence  is  more  than  an  order  of  magnitude 
less  than  predicted  by  Eq.  (6).  We  do  not  know 


stimulated  Raman  scattering  (SRS)  and  stimulated  colli¬ 
sion-induced  fluorescence  (SCF)  from  spontaneous, 
linear  scattering  in  the  forward  direction.  The  solid 
curves  are  theoretical  curves  obtained  from  Eqs.  <5) 
and  (6)  for  the  SRS  and  SCF,  respectively,  (b)  Spon¬ 
taneous,  linear  scattering  in  the  side  direction  for 
comparison.  The  laser  intensity  4°  was  2  MW/cm2. 


x  10*  rad  sec'1.  This  Is  a  factor  of  3  larger  than 
the  value  determined  by  the  emission  line-shape 
measurements  of  Cheron,  Scheps,  and  Gallagher.1* 
According  to  recent  results  of  Nienhuis  and  Schul¬ 
ler,"  we  expect  the  two  experiments  (collisional 
redistribution  and  collisional  line  broadening)  to 
give  the  same  value  for  ys(A)  over  the  experi¬ 
mental  range  that  we  studied.  The  existing  dis¬ 
crepancy  may  be  due  to  the  transient  nature  of 
our  experiment  as  well  as  neglect  of  degeneracy 
in  Eq.  (7). 

In  the  forward  direction,  as  in  the  side  direc¬ 
tion,  the  emission  is  initially  linear  at  low  inten¬ 
sities.  However,  eventually  both  the  Raman 
emission  and  the  collision-induced  fluorescence 
grow  exponentially  when  the  scattering  becomes 
stimulated,  and  at  still  higher  intensity  satura¬ 
tion  occurs.  Using  Eqs.  (5)  and  (6)  multiplied 
by  the  pulse  duration  we  calculated  the  expected 
exponential  growth  for  the  Raman  scattering  and 
collision-induced  fluorescence.  The  results  are 
shown  as  the  solid  curves  in  Fig.  2(a).  The  verti¬ 
cal  scale  was  considered  a  free  parameter,  but 
agreed  to  within  a  factor  of  .4  with  an  absolute  - 
calibration  estimate.  The  laser  profile  was 
taken  to  be  uniform  in  intensity  over  the  300- pm 


the  reason  for  the  discrepancy,  but  problems  may 
arise  from  use  of  the  steady -«tate  theory  as  well 
as  from  the  assumption  of  a  constant  spatial  pro¬ 
file. 

In  order  to  study  further  the  predictions  of 
Eqs.  (5)  and  (6),  we  measured  the  dependence  of 
both  the  SRS  and  SCF  as  a  function  of  Ar  pres¬ 
sure.  Profiles  similar  to  those  in  Fig.  2(a)  were 
taken  at  5,  20,  and  80  Torr  of  Ar  buffer  gas,  at 
both  * 0.14  nm  detuning.  While  the  SRS  showed 
no  pressure  dependence,  the  SCF  was  highly  pres¬ 
sure  dependent,  as  expected  from  the  theoretical 
discussion.  In  particular,  the  collision-induced 
fluorescence  was  linear  in  Ar  pressure  at  low 
intensities,  while  the  high-intensity  behavior  is 
described  by  noting  that  the  SCF  gain  was  roughly 
linear  in  Ar  pressure.  Thus,  curves  calculated 
from  Eq.  (6),  but  with  the  gain  decreased  by  an 
overall  factor  of  30,  were  in  good  agreement 
with  these  data.  These  results  show  that  the  ex¬ 
perimental  dependence  of  the  SRS  and  SCF  on  the 
collision  rate  was  in  accord  with  theory. 
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Abstract.  We  report  preliminary  studies  on  the  self-broadening  of  the  377-6  nm  thallium 
line.  In  the  short  wavelength  wing  we  observe  asymmetries  and  a  satellite  at  about 
377- 1  nm.  In  the  long  wavelength  quasi-static  wing  we  clearly  observe  the  transition  from 
the  Lorentzian  to  the  van  der  Waals  type  wing. 


In  recent  studies  of  self-broadening  of  the  quasi-static  wings  of  the  resonance  lines 
in  the  alkali  elements,  satisfactory  agreement  has  been  obtained  between  experiment 
(Niemax  and  Pichler  1975)  and  theory  (Movre  and  Pichler  1977).  In  this  letter  we 
describe  an  extension  of  these  studies  to  thallium,  a  group  three  element  which  is 
essentially  also  a  one-electron  system.  Some  early  studies  of  the  spectra  of  thallium 
vapour  in  emission  were  made  by  Hamada  (1931).  Since  the  ground  state  of  Tl 
is  a  p  level  and  the  first  excited  state  is  an  s  level  (see  figure  1),  we  expected  similarities 
with  the  alkalis  in  the  self-broadening  of  the  resonance  lines.  In  particular,  in  studying 
the  62P, ;2-72S1/2  transition  at  377-6  nm  in  thallium  we  expected  that  the  main  results 
from  the  study  of  the  alkali  resonance  lines  could  be  applied.  The  similarity  arises 
since  the  interaction  potential  curves  emerging  from  the  62P,/2  +  72S,;2  level  are 
similar  to  the  potential  curves  for  the  ns  2S1/2  +  np  2P1/2  resonance  levels  of  the 
alkalis.  However,  one  essential  difference  is  very  apparent — the  very  large  spin-orbit 
interaction  in  the  ground  levels  of  thallium.  Initially  we  expected  that  the  inner  wing 
satellites,  observed  near  the  heavy  alkali  resonance  line  doublets,  would  occur  very 
far  into  the  line  wings,  based  on  figure  3  in  the  paper  of  Movre  and  Pichler  (1977). 
However,  no  satellite  in  the  far  blue  wing  spectrum  of  the  377-6  nm  line  of  thallium 
was  observed.  Analysis  has  shown  that  this  is  due  to  the  fact  that  the  oscillator 
strength  of  the  377-6  nm  line  is  very  small  (/  =  0T4),  which  together  with  the  short 
wavelength,  results  in  a  smaller  interaction  constant  for  the  resonance  interaction 
expected  by  analogy  to  the  alkalis  (C3  oc  //.).  Instead  we  observed  the  blue  satellite 
structure  much  closer  to  the  centre  of  the  line.  This  letter  reports  our  preliminary 
findings. 

In  our  experiment  we  used  a  cw  tungsten-halogen  lamp  with  current  stabilisation 
to  study  the  absorption  profile  of  the  3776  A  line  of  Tl.  The  light  beam  was  collimated 
and  sent  through  a  heat-pipe  oven  having  an  approximate  10  cm  column  of  pure 
thallium  vapour.  Argon  was  used  as  a  buffer  gas  to  protect  the  quartz  windows 
from  the  hot  thallium  vapour.  The  working  temperature  range  was  in  the  interval 
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Figure  t.  Partial  term  diagram  for 
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between  900  and  1100CC,  in  which  the  thallium  vapour  pressure  was  changed  by 
means  of  changing  the  argon  pressure  in  the  interval  of  1  and  20  Torr.  In  order 
to  determine  the  centre  of  the  strongly  absorbed  377-6  nm  line,  we  used  a  low  current 
T1  emission  lamp.  During  a  spectral  scan  of  the  377-6  nm  line,  the  absorption  in 
the  heat-pipe  oven  was  saturated  within  the  region  of  about  0-05  nm  around  the 
line  centre.  In  this  spectral  region  a  folding  mirror  was  turned  to  allow  the  light 
from  the  T1  emission  lamp  to  enter  the  monochromator  and  mark  the  centre  of 
the  resonance  line.  The  light  coming  out  from  the  absorption  vessel  and  the  T1 
emission  lamp  was  focused  on  the  entrance  slit  of  a  scanning  monochromator 
(McPherson,  Model  EU-700)  with  resolution  of  0-05  nm.  A  photomultiplier  with  an 
S-20  photocathode  was  used  for  the  light  detection.  The  signal  was  amplified  by 
means  of  the  lock-in  amplifier  and  recorded  on  an  XT  recorder. 

In  the  preliminary  results  reported  here  we  have  not  measured  the  length  of 
the  thallium  vapour  column,  which  varied  when  we  changed  the  electrical  power 
into  the  oven  or  the  argon  pressure.  Hence  the  data  have  been  analysed  in  terms 
of  optical  depth  r  =  kL  and  therefore  have  no  absolute  values  of  the  absorption 
coefficientt.  A  typical  absorption  profile  is  shown  in  figure  2  for  6  3  Torr  of  thallium 
pressure.  About  0-5  nm  from  the  line  centre,  toward  shorter  wavelengths,  a  broad 
shoulder  appears,  beyond  which  there  is  a  sharp  decrease  of  absorption.  We  believe 
this  satellite  is  formed  by  the  0,*  state  as  discussed  below.  Even  closer  to  the  line 
centre  the  red  wing  is  more  absorbed  than  the  blue  wing.  However,  our  instrumental 
width  precluded  any  detailed  measurements  within  the  region  of  the  line  centre.  In 
order  to  verify  that  the  oven  was  running  in  a  heat-pipe  mode  and  that  the  broadening 
was  due  to  T1  T1  collisions,  and  not  Tl-Ar  collisions  we  increased  the  argon  pressure 
to  the  point  where  the  effect  of  the  Tl-Ar  collisions  could  be  seen.  We  found  that 
the  absorption  lineshape  was  then  consistent  with  a  previous  study  by  Cheron  ci 
ul  (1976.  1977). 

We  have  analysed  our  data  by  plotting  the  logarithm  of  optical  depth  t  against 
the  logarithm  of  the  wavelength  separation  A/,  from  the  line  centre.  Different  sets 

♦  In  order  to  obtain  k  V 2  (ctn’l  values  (absorption  coefficient  divided  by  the  square  of  the  particle  density! 
it  is  necessary  to  multiply  our  relative  r  values  given  in  figure  '  by  a  factor  of  (I -ft  +  0-51  »  10 
which  was  obtained  by  a  rough  estimate  of  /.  =  10  cm. 


Figure  2.  Typical  absorption  profile  of  the  T1  377-6  nm  line.  Arrow  shows  the  line  centre 
determined  by  the  emission  line.  The  temperature  was  1000  C  and  the  pressure  of  the 
thallium  vapour  was  6-3  Tort. 

of  data  were  put  together  by  sliding  along  the  log  t  axis.  Figure  3  shows  the  resulting 
diagram.  The  satellite  in  the  blue  wing  at  0-5  nm  is  clearly  seen.  Since  an  exponential 
■  fall-off  is  expected  after  a  satellite  (Sando  and  Wormhoudt  1973,  Szudy  and  Baylis 

I  1975,  Carrington  and  Gallagher  1974)  we  also  analysed  the  optical  depth  beyond 

I 
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the  satellite  at  0-5  nm  in  terms  of  a  semi-logarithmic  plot  (log  t,  A/).  Exponential 
behaviour  was  verified  up  to  0-8  nm  at  a  thallium  pressure  of  20  Torr. 

Careful  study  of  figure  3  reveals  that  the  red  wing  in  the  region  from  01  to 
0-5  nm  from  the  line  centre  is  approximately  Lorentzian  [r  ~  (A/)-2]  and  then  it 
continuously  converges  to  a  van  der  Waals  type  wing  (r  -  (A2)  1  *,  see  e.g.  Mihalas 
1978)  up  to  381  nm  where  a  red  satellite  appears.  This  change  from  Lorentzian  to 
van  der  Waals  dependence  is  indicative  of  the  importance  of  the  van  der  Waals 
type  interaction  in  the  shape  of  the  resonance  potential  curves.  As  will  be  discussed 
below  this  van  der  Waals  interaction  is  also  of  great  importance  in  the  formation 
of  the  satellite  in  the  blue  wing. 

In  figure  4  we  have  reproduced  the  resonance  interaction  potential  curves  (broken 
curves)  emerging  from  the  72S,  ,  +  62P,  2  level  in  a  thallium  collision-induced  mol¬ 
ecule  (Fontana  1962,  Movre  and  Pichler  1977)  with  the  following  C3  interaction 
constants : 


C3(0g+)  =  -C3(0U+)  =  — 

C3(1J=  -C3(lg)  =  ~- 

C3(0g-)  =  C3(0U-)  =  0 
where  C  =  e2|<0|r|  1  >|2  and 


K  0|  r  1 1  >  | 2 


9hfX 
8  n2mc 


(1) 

(2) 

(3) 

(4) 


Figure  4.  Resonance  interaction  curves  emerging  from  the  72S,  2  +  62P,  2  molecular  level 
(broken  curves)  and  the  bending  of  these  curves  caused  by  the  van  der  Waals  interaction 
(full  curves). 
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Using  an  /-value  of  0-14  (Gallagher  and  Lurio  1964,  Penkin  et  al  1973)  for  the  Tl 
377-6  nm  line  we  obtain  interaction  constants  (in  rad  s  “ 1  cm  J)  for  the  repulsive 
states  as  follows 


C3(0,+  )  =  14-2  x  10" 9 

(5) 

C3(l„)  =  7-1  x  10'9. 

(6) 

The  full  curves  in  figure  4  show  how  the  resonance  potential  curves  bend  down 
under  the  influence  of  the  attractive  van  der  Waals  interaction.  In  this  way  the  two 
maxima  H( 0,+  )  and  H{ lu)  are  formed.  We  took  the  wavelength  separation  of 
0-5  ±  0-2  nm  as  a  position  where  the  satellite  has  approximately  65%  of  its  peak 
intensity  (Carrington  and  Gallagher  1974,  and  references  therein),  which  corresponds 
to  the  maximum  H  in  the  0*  potential  curve.  We  assume  that  this  maximum  is 
formed  by  the  competition  of  the  R~3  and  Rb  terms  in  the  potential  expansion 
(we  also  assume  implicitly  that  the  ground  state  is  flat  in  this  long-range  region): 

^  ~  ^  (7) 


By  means  of  the  condition  for  an  extremum 

d(Am) 

dR 


(8) 


we  can  obtain  both  the  C6  constant  and  the  internuclear  separation,  where  the  extre¬ 
mum  occurs,  by  means  of  the  following  relations: 

=  H  -  Cj/4C6  (9) 

and 

-  (2C6/C3),/j.  (10) 

From  (9),  using  the  C3  constant  from  (5)  and  Ax,,,  =  0-5  +  0-02  nm  for  the  maximum, 
we  obtain  (in  rad  s  " 1  cm6) 

C6  =  (7-6  ±  0-3)  x  10' 30 

and  from  (10), 

Kn,.«(0,+  )=  10-2  A. 

We  have  used  the  same  C6  constant  in  order  to  estimate  the  position  of  the  maximum 
//( lu).  The  analysis  indicates  that  it  should  appear  at  A/  =  0-12  nm  in  the  blue  wing 
(with  Rma,  =  12-8  K). 

Unfortunately,  the  resolution  of  our  monochromator  did  not  allow  us  to  resolve 
the  spectral  region  where  this  satellite  should  appear,  but  we  have  observed  the 
red  asymmetry  in  the  wings  outside  this  region,  which  would  occur  due  to  the 
exponential  fall-off  beyond  the  satellite. 

We  have  estimated  the  limit  between  the  quasi-static  and  impact  approximations 
using  both  the  C3  and  C6  constants.  We  found  that  the  region  around  AA  =  0-1  nm 
is  still  close  to  the  edge  of  the  impact  region,  a  fact  which  will  broaden  the  satellite 
related  to  the  maximum  ff(lu).  An  additional  broadening  of  the  satellite  could  come 
from  the  mixing  among  near-by  slates  in  the  72S,/2  +  62P1/2  manifold. 
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Using  the  same  continuum  background  source  and  quartz  optics  we  have  also 
made  several  scans  across  the  T1  276-8  nm  line  (62D3/2-62P1/2  transition)  and  we 
have  observed  one  blue  satellite  even  closer  to  the  line  centre  at  about  AA  =  0-18  nm. 
Beyond  this  satellite  the  steep  decrease  of  the  absorption  was  observed  and  a  very 
extended  red  wing.  We  have  not  studied  this  line  in  detail  since  there  is  no  analogous 
line  in  the  resonance  lines  of  the  alkali  atoms.  We  intend  to  extend  calculations 
of  the  potential  curves  for  this  line  in  the  future. 

We  would  like  to  point  out  one  major  difficulty  encountered  here  in  the  thallium 
self-broadened  resonance  lines,  which  is  not  present  in  alkali  cases.  In  figure  1  we 
have  reproduced  a  partial  term  diagram  of  atomic  thallium  and  a  diatomic  thallium 
collisional  molecule.  All  strong  resonance  lines  are  marked  in  both  diagrams.  How¬ 
ever  we  immediately  see  that  all  atomic  lines  can  be  generated  in  the  diatomic  mol¬ 
ecule  from  two  distinctively  different  molecular  levels.  If  both  transitions  could  occur 
at  the  same  time  then  one  spectral  line  would  bear  information  about  the  splitting 
structure  of  all  four  levels  involved.  In  absorption  measurements  this  ambiguity  is 
easily  avoidable  since  thermal  population  of  the  62P3/2  level  at  temperatures  of  about 
1000°C  is  still  relatively  low.  Nevertheless  we  have  observed  all  five  resonance  lines 
in  absorption,  but  the  535-0  nm,  351-9  nm,  and  352-9  nm  lines  have  not  been  absorbed 
enough  to  allow  broadening  studies  in  the  quasi-static  wings. 

The  authors  acknowledge  the  helpful  comments  of  A  Gallagher.  This  work  was  par¬ 
tially  supported  by  National  Science  Foundation  Grant  No  PHY76-04761  and  Office 
of  Naval  Research  Contract  No  N00014-76-C-0611,  both  through  the  University  of 
Colorado. 
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Abstract.  We  have  studied  collisional  redistribution  of  near-resonance  radiation  in  thal¬ 
lium  vapour  caused  by  collisions  with  argon  buffer  gas.  We  spectrally  resolved  electronic 
Raman  scattering  and  collision-induced  fluorescence,  both  to  the  weakly  populated 
metastable  61P3/2  level  in  thallium.  Thus  we  partially  avoided  the  problem  with  radiative 
trapping  encountered  in  previous  experiments.  By  measuring  the  total  collision-induced 
fluorescence  as  a  function  of  the  incident  laser  frequency,  we  obtained  the  absolute 
collisional  redistribution  line  shape.  We  have  compared  this  line  shape  with  the  emission 
tine  shape  measured  by  Cheron,  Scheps  and  Gallagher. 


The  combined  effect  of  collisional  interaction  and  radiative  interaction  upon  an  atom  in 
a  vapour  continues  to  be  a  subject  of  interest  (Nayfeh  etal  1977,  Ballagh  and  Cooper 

1977,  Raymer  and  Carlsten  1977,  Voslamber  and  Yelnik  1978).  A  question  of  basic 
importance  is:  what  is  the  spectrum  of  the  scattered  radiation  when  monochromatic 
light  from  a  laser  is  incident  upon  a  vapour  in  which  collisions  are  present?  Although  it 
is  intuitively  expected  that  this  spectrum,  known  as  the  redistribution  function,  is 
related  in  some  manner  to  the  absorption  and  emission  profiles  obtained  from  con¬ 
ventional  studies,  it  is  not  known  what  the  explicit  relationship  is  in  the  general  case.  At 
present  there  are  several  theoretical  treatments  which  address  this  problem  (Cooper 

1978,  Voslamber  and  Yelnik  1978,  Nienhuis  and  Schuller  1977,  Huber  1969).  In  this 
letter  we  present  the  results  of  a  redistribution  experiment  and  compare  them  to  those 
of  an  emission  experiment  by  Cheron  et  al  (1977). 

To  describe  the  redistribution  process  we  refer  to  figure  1,  which  shows  the  three 
lowest  energy  levels  of  atomic  thallium.  A  low-intensity  laser  is  incident  at  frequency 
«L,  near  the  zPi/2-JSj/2  resonance  line  and  emission  occurs  at  four  frequencies: 
Rayleigh  scattering  at  wray.  electronic  Raman  scattering  at  coRAM  and  collision- 
induced  fluorescence  at  both  0121  and  o>23-  The  collisions  can  be  thought  of  as  providing 
the  energy  needed  to  transfer  population  from  the  laser-induced  virtual  level  (broken 
line)  to  the  real  level  2.  The  process  can  also  be  viewed  in  a  molecular  picture  as  in 
figure  2,  which  shows  the  potential  curves  for  a  Tl-Ar  collision  (Cheron  etal  1977).  The 
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Figure  1.  The  three  lowest  levels  of  atomic  thallium. 
A  low-intensity  laser  of  frequency  aiL  incident  on  a 
vapour  of  three-level  atoms  generally  produces  four 
components  in  the  scattered  spectrum:  Rayleigh 
scattering  at  u>ray( =<•'!.).  electronic  Raman  scatter¬ 
ing  at  uram(=«l_"3i)  and  collision-induced 
fluorescence  at  both  uj,  and  *)u. 


Figure  2.  Potential  curves  for  Tl-Ar  (Cheron  er  al 
1977),  showing  collisional  redistribution,  viewed  at 
absorption  by  the  quasi-molecule  formed  during  the 
collision .  The  excited  11  atom  then  radiates  at  w2i or 
w 2 3  after  the  collision  is  completed. 


interaction  with  an  Ar  perturber  causes  the  T12Si/2  level  to  come  into  resonance  with 
the  laser  frequency  (for  the  red  detuning  shown)  at  some  internuclear  separation, 
allowing  absorption  of  a  photon  and  creating  an  excited  quasi-molecule.  Usually  the 
quasi-molecule  has  enough  internal  kinetic  energy  to  fly  apart  into  a  ground-state  Ar 
and  an  excited  T1  atom,  which  can  then  emit  a  photon  of  frequency  ton  or  on.  This  type 
of  collision  has  been  called  a  dephasing  collision  or  an  elastic  (or  quasi-elastic)  collision 
(Omont  et  al  1972,  Mollow  1977). 

The  theory  of  light  scattering,  in  the  presence  of  collisions,  is  well  understood  at  low 
intensities  when  the  incident  and  scattered  light  have  frequencies  very  close  to  the  line 
centre  (in  the  impact  region).  In  this  region  the  line  shape  is  Lorentzian,  regardless  of 
the  form  of  the  interaction  potentials.  Referring  to  figure  1 ,  the  normalised  absorption 
or  emission  profile  for  the  1-2  transition  is  given  by  Omont  et  al  (1972)  as 


/(A,) 


(rE+yN)/ir 
Ai  +  tyE+yN)^ 


(l) 


where  2>n  =  2-yN  +  2y”  is  the  sum  of  the  spontaneous  decay  rates  from  level  2  to  levels 
1  and  3,  2yE  is  the  rate  of  elastic  collisions,  and  Ai  =  on  -ton  is  the  detuning  of  the 
incident  light  from  resonance. 

In  a  redistribution  experiment,  the  differential  cross  section  (in  units  cm2  rad"1  s)for 
absorbing  a  photon  with  frequency  om  and  emitting  one  with  frequency  a>:  (integrated 
over  outgoing  direction  and  polarisation)  can  be  written  as  (Omont  et  al  1972,  Batlagh 
and  Cooper  1977): 


g  2  9 

- — <t(u»i,  <uj)  =  (2n  e  //mc)f(u>i  -  u>ji)P(u>i,  wj) 
do»2 


(2) 
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where 


P(t ill,  «j)  =  ( — —  S(<02  -  0>i)  + ^ /(«I*2  -  <i»2»)) 

ys'ye+yN  yE+yN  • 

yE  \ 

2  —  0>l  +  Ct)Jl)  + - - - ton)  I 

ye+yN  / 


y”/ 

'  yN 

ys' 

► 

+ 

ID 

(3) 


and  /(&>)  is  always  the  function  in  equation  (1).  Here  P(o>i,  on)  is  the  redistribution 
function.  Given  that  a  photon  of  frequency  <oi  has  been  absorbed,  P(u> i,  w 2)  gives  the 
probability  that  one  will  be  emitted  at  w2-  Since  we  have  ignored  quenching  collisions, 
we  have  J  P(«m,  «j2)  dw2  =1. 

The  terms  in  equation  (3)  proportional  to  yl!  describe  the  emission  near  a>2i:  the 
Rayleigh  scattering  is  a  delta  function  at  wi(=a»L)  while  the  collision-induced  fluores¬ 
cence  is  a  Lorentzian  centred  at  «2i.  Similarly,  the  terms  proportional  to  yl?  describe 
the  Raman  scattering  at  o>ram  and  fluorescence  centred  at  <023-  If  equation  (2)  is 
integrated  over  both  the  incident  frequency  a>t  and  the  scattered  frequency  a »2  one 
obtains  the  usual  result  for  the  total  cross  section  <xtot  =  / 1  me,  where  /  is  the 

absorption  oscillator  strength  for  the  1-2  transition  (Mihalas  1970).  In  order  for 
equation  (3)  to  be  valid,  the  broadening  of  levels  1  and  3  must  be  negligible  compared  to 
that  of  level  2.  When  the  incident  light  is  far  off  resonance  (Ai  »yE  +  yN),  the  four 
components  can  be  spectrally  resolved  and  integrated  separately  over  <o2  to  give  a  cross 
section  (in  units  of  cm2)  for  each  process; 


o'Ray  =  OTOTfyN/rrA2) 

Rayleigh 

(4a) 

o-f1  =  o’TOTfyN /yN)(yE/irA?) 

2-1  fluorescence 

(46) 

oram  -  Oror(yN3  /  tr  A 1) 

Raman 

(4c) 

0”  =oTor(yN /yN)(yE/oA?) 

2-3  fluorescence. 

(4cf) 

These  cross  sections  can  be  used  to  calculate  the  scattered  intensity;  for  example,  the 
2-3  fluorescence  intensity  /”  (in  units  of  photons/cm3/s)  is  equal  to  N^aVli.,  where 
N-n  is  the  atomic  thallium  density  (cm-3)  and  IL  is  the  incident  laser  intensity  (in  units  of 
photons/cm2/s).  Thus  the  profile  for  collisionally  redistributed  fluorescence  is  Lorent¬ 
zian  in  the  impact  region. 

For  larger  detunings,  outside  the  impact  region,  the  situation  is  not  completely 
understood.  It  is  known  that  the  absorption  (emission)  profile  becomes  non-Lorent- 
zian,  and  Cooper  (1978)  has  shown  it  can  be  obtained  by  replacing  yE  in  equation  (1)  by 
TeIAi) 


/(A.)  = 


(yE(Ai)  +  yN)/tr 

A?  +  (yE(Ai)  +  yN)2 


(5) 


where  yE(A0  is  a  frequency-dependent  collision  rate.  It  has  been  observed  (Carlsten  et 
al  1977)  that  the  redistribution  profile  also  becomes  non-Lorentzian.  Cooper  (1978) 
has  predicted  that,  when  collisions  dominate  the  spontaneous  emission  (yE(Ai) »  yN), 
the  non-impact  redistribution  profile  may  also  be  obtained  by  replacing  yE  (in  equation 
(4d))  by  yE(Ai): 
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The  condition  yE(A|) »  yn  ensures  that  the  redistributed  light  comes  from  uncorrelated 
absorption  and  emission  events  and  thus  will  be  isotropic  and  unpolarised.  When  this 
condition  is  not  upheld,  the  angle-  and  polarisation-resolved  redistribution  cross 
section  should,  in  principle,  contain  more  information  about  the  mixing  of  the  atomic 
states  during  the  collision  than  the  absorption  cross  section  (Cooper  1978).  When 
Ye(Ai)  »  yN,  the  absorption  (emission)  profile,  equation  (5),  can  be  approximated  in  the 
line  wing  as 


/(A,) 


ye(A,) 

TaT 


(7) 


Thus  when  collisions  dominate  spontaneous  emission,  the  redistribution  profile  (6)  and 
the  absorption  (emission)  profile  (7)  are  predicted  to  have  the  same  shape. 

To  study  the  relation  between  emission  and  redistribution  we  chose  the  Tl-Ar 
system  for  two  reasons:  it  has  already  been  studied  by  Cheron  et  al  (1977)  using  the 
emission  technique,  and  it  avoids  some  of  the  experimental  problems  associated  with 
radiative  trapping.  In  the  experiment  by  Cheron  etal,  the  T1  resonance  line  was  excited 
at  its  line  centre  by  a  spectrally  narrow  light  source  and  the  spectrum  of  fluorescence  in 
the  line  wings  was  measured.  Since  in  that  experiment  collisions  dominated  spon¬ 
taneous  emission,  the  emission  profile  they  observed  was  essentially  that  in  equation 
(7).  In  a  previous  study  (Carlsten  et  al  1977)  we  observed  redistribution  in  a  two-level 
atom  Sr,  but  the  fluorescence  was  heavily  trapped,  making  an  absolute  measurement  of 
Ye(Ai)  difficult.  By  detecting  fluorescence  at  o»23  in  T1  we  avoid  part  of  this  problem  (see 
later  in  the  discussion),  as  the  2P3/2  level  is  essentially  unpopulated  at  our  operating 
temperature. 

The  apparatus  used  in  this  experiment  is  shown  in  figure  3  and  is  similar  to  that 
described  by  Carlsten  etal  (1977).  A  dye  laser,  pumped  by  an  N2  laser,  was  tuned  in  the 
region  of  the  62Pi/2-72Si/2  TI  resonance  line  at  377-6  nm.  Emission  was  observed  near 
the  72St/2-62P3/2  line  at  535-0  nm.  The  dye  laser  had  a  pulse  duration  of  10  ns,  a 
spectral  width  of  0-03  nm  and  an  energy  of  about  50  The  beam  was  directed  into  an 
oven  at  1000  K  containing  0-1  Torr  of  Tl  vapour  and  20  Torr  of  Ar  buffer  gas.  The 
column  length  of  Tl  was  approximately  10  cm.  As  our  present  oven  did  not  allow 
observation  at  right  angles  to  the  laser  beam,  we  chose  to  detect  back-scattered  light  by 
use  of  a  dichroic  mirror,  which  transmitted  the  uv  laser  light,  but  reflected  the  green 
fluorescence  signal.  A  cut-off  filter  was  also  inserted  to  further  reject  laser  light.  The 
fluorescence  was  then  focused  onto  the  slits  of  a  monochromator  of  0-06  nm  resolution 
with  an  f/ 10  optical  system.  The  signal  from  a  photomultiplier  was  boxcar  averaged 
and  the  spectral  profiles  recorded  on  a  chart  recorder. 
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The  ratio  of  fluorescence  intensity,  7f\  to  Raman  intensity,  /ram,  was  measured  at 
each  incident  laser  frequency  and  the  quantity  (/“//ramKyn/wA^/V  ‘2wc  is  plotted 
as  points  in  figure  4,  where  N  is  the  number  density  of  argon  and  yn  is  the  half -width  (in 
units  of  rad  s~l)  due  to  natural  broadening  (Gallagher  and  Lurio  1964).  From  equa¬ 
tions  (4c)  and  (6)  one  can  see  that  this  quantity  should  have  the  same  shape  as  /(A,)  in 
equation  (7).  The  constant  factors  guarantee  normalisation  in  units  of  cm4.  Because  we 
measured  only  the  ratio  /”//ram.  and  not  /”  separately,  one  can  say  that  we  used  the 
Raman  to  monitor  our  dye  laser  intensity  as  we  moved  it  through  its  tuning  range.  This 
procedure  relies  on  the  well  founded  assumption  (Cooper  1978,  Carlsten  et  al  1977) 
that  the  Raman  is  Lorentzian  over  the  entire  spectral  range  (i.e.  equation  (4c)  holds 
even  outside  the  impact  region).  The  linearity  of  the  signal  with  laser  power  was 
checked  to  verify  that  there  were  no  saturation  or  stimulated  effects  present.  This  test  is 
important  because,  at  higher  laser  power,  stimulated  emission  at  both  wram  and  cj2j  has 
been  observed  to  occur  (Raymer  and  Carlsten  1977). 

at  icm'1) 


Figure  4.  Line  shape  of  the  JPl/2_IS|/2  transition  in  Ti,  broadened  by  Ar,  obtained  by  two 

techniques:  •,  collisional  redistribution  (present  work);  - ,  emission  spectroscopy 

(Cheron  et  al  1977).  The  broken  curve  is  the  absorption  profile,  obtained  from  the  emission 
profile  by  multiplying  by  a  Boltzmann  factor  (see  text). 


Also  plotted  in  figure  4  as  a  full  curve  is  the  normalised  emission  profile  of  Cheron  et 
al  (1977).  The  agreement  between  this  emission  profile  and  our  redistribution  profile  is 
striking.  The  difference  at  large  detunings  on  the  red  side  (+AA)  is  due  to  a  tempera¬ 
ture-dependent  population  effect,  discussed  by  Hedges  et  al  ( 1972).  According  to  their 
discussion,  we  should  really  compare  our  data  with  an  absorption  profile,  obtained  by 
multiplying  the  emission  profile  by  the  factor  exp(A k/kT),  where  A k  is  in  cm-1  and 
kT  =  5 16  cm-1  at  743  K,  the  temperature  at  which  the  emission  experiment  was  done. 
This  derived  absorption  profile  is  shown  as  a  broken  curve  in  figure  4.  However,  since 
our  redistribution  experiment  was  done  at  a  higher  temperature  (1000  K),  it  is  not 
surprising  that  our  points  lie  above  the  derived  absorption  profile,  because  our 


LI  24 


Letter  to  the  Editor 


collisions  at  higher  temperature  should  probe  farther  up  the  repulsive  wall  of  the 
ground  state  (figure  2)  than  those  at  the  lower  temperature  used  in  the  other  experi¬ 
ment.  The  effect  on  the  blue  side  (-AA )  is  less  pronounced  (and  not  plotted)  due  to  the 
steep  fall-off  in  the  line  shape.  The  deviation  beyond  30  A  on  the  blue  side  is  believed 
to  be  caused  by  impurities  in  the  oven,  as  these  points  do  not  scale  linearly  in  argon 
pressure  as  the  rest  of  the  profile  does.  The  present  data  were  all  taken  outside  the 
impact  region  (due  to  instrumental  resolution)  and  so  we  never  reached  the  limit,  near 
line  centre,  where  the  red  and  blue  sides  of  the  line  become  symmetric.  A  complete 
discussion  of  the  line  shape  in  figure  4  is  gr'en  by  Cheron  et  al  (1977)  and  will  not  be 
repeated  here.  Clearly,  it  is  related  to  the  detailed  shapes  of  the  potential  curves  shown 
in  figure  2. 

Although  the  line  shapes  agree  very  well,  there  is  uncertainty  in  our  experiment 
about  the  absolute  scale  in  figure  4.  We  plotted  our  data  as  if  the  emission  were  coming 
from  a  single  T1  atom  which  behaves  in  the  way  described  by  equation  (3).  That  is,  the 
branching  ratios  for  the  2-3  and  2-1  transitions  are  ynVyn  and  YnVyn-  However, 
when  we  took  for  emission  at  <u2i  we  find  that  it  is  trapped  by  the  surrounding 
ground-state  atoms.  If  all  of  this  trapped  radiation  is  subsequently  emitted  at  «23  then 
we  should  see  a  fluorescence  signal  at  wj3  enhanced  by  the  factor  (yn3  +yn  )/yn3  = 
2-01.  This  interpretation  would  require  us  to  lower  all  of  our  points  by  the  factor 
(2-01)-1.  This  uncertainty  must  be  resolved  before  the  absolute  cross  sections  from  the 
two  experiments  can  be  compared  more  accurately. 

In  summary,  we  have  demonstrated  the  relationship  between  collisional  redistri¬ 
bution  and  emission  line  profiles.  We  found  that  in  the  TI-Ar  system  these  two  profiles 
are  essentially  the  same  when  collisions  dominate  spontaneous  emission.  When 
collisional  and  radiative  rates  are  comparable,  systems  may  exist  in  which  a  difference  in 
the  line  shapes  could  be  found,  due  to  angular  and  polarisation  effects  in  the  redistribu¬ 
tion  process. 

We  would  like  to  acknowledge  helpful  discussions  with  R  J  Ballagh,  J  Cooper,  A 
Gallagher  and  A  Szoke.  This  work  was  supported  by  Office  of  Naval  Research  Contract 
No  N00014-76-C-0611  and  National  Science  Foundation  Grant  No  PHY76-04761, 
both  through  the  University  of  Colorado. 
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Theory  of  stimulated  Raman  scattering  with  broadband  lasers 
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The  authors  have  extended  the  theory  of  stimulated  Ramar  scattering  to  include  the  effects  of  laser 
bandwidths,  in  both  the  transient  and  steady-state  regimes.  The  rase  o."  two  interacting  laser  beams,  a  pump 
laser  and  a  probe  (Stokes)  laser,  is  treated.  Using  the  phase-diffusion  model  for  laser  bandwidth,  the  authors 
demonstrate  that  in  the  absence  of  dispersion,  the  forward  Raman  gain  is  essentially  independent  of  the  laser 
bandwidths  in  the  high-gain  limit,  while  in  the  low-gain  limit  the  gain  coefficient  is  inversely  proportional  to 
the  sum  of  the  bandwidths.  It  is  further  shown  that  when  the  pump-laser  bandwidth  is  much  larger  than  the 
linewidth  of  the  Raman  medium,  the  stimulated  Stokes  output  assumes  the  same  spectrum  as  the  pump  laser 
in  the  high-gain  limit.  A  possible  interpretation  of  these  results  is  discussed  assuming  a  “phase-locking”  of 
the  Stokes  phase  to  the  fluctuations  in  the  pump-laser  phase,  due  to  the  nonlinear  interaction  of  the  two 
beams  through  the  Raman  medium. 


I.  INTRODUCTION 

The  effects  of  finite  l.iser  bandwidth  are  being 
recognized  as  important  in  the  study  of  nonlinear 
optical  processes.  Resonance  fluorescence,1** 
two-photon  absorption,  7,3  second  harmonic  genera¬ 
tion,*  multiphoton  ionization,10*13  and  stimulated 
Raman  scattering  are  all  areas  in  which  key  ele¬ 
ments  of  understanding  depend  on  the  ability  to 
model  the  laser,  not  as  a  monochromatic  wave  of 
definite  phase  and  amplitude,  but  as  a  multimode 
broad-band  wave  with  fluctuating  phase  and  ampli¬ 
tude.  The  problem  of  stimulated  Raman  scattering 
(SRS)  is  especially  timely  in  light  of  ongoing  ef¬ 
forts  to  use  it  as  method  for  developing  new  co¬ 
herent  light  sources1*  as  well  as  compressing 
high-energy  laser  pulses  to  achieve  higher  peak 
powers  for  use  in  laser  fusion.13  In  applications 
of  these  types  a  detailed  understanding  of  all  the 
factors  influencing  the  efficiencies  of  the  process¬ 
es  is  obviously  desirable.  However,  one  import¬ 
ant  factor,  laser  bandwidth,  has  not  yet  been  fully 
explored. 

Recently,  two  groups13* l*  have  observed  a  large 
forward-backward  asymmetry  of  the  Raman  gain, 
which  they  attribute  to  the  broad-band  nature  of 
the  pump  laser  used.  These  were  in  the  absence 
of  other  effects,  such  as  self-focusing  or  extran¬ 
eous  feedback,  which  are  known  to  produce  anom¬ 
alous  gains.17  These  asymmetries  are  consistent 
with  several  theoretical  predictions13* 19  that  in  the 
backward  direction  (counterpropagating  pump  and 
Stokes  waves)  the  gain  coefficient  is  proportional 
to  (r>r£)*',  where  r  and  r£  are  the  spectral 
widths  of  the  Raman  medium  and  the  pump  laser, 
respectively;  while,  in  the  forward  direction,  in 
the  absence  of  dispersion  of  the  Stokes  wave  re¬ 
lative  to  the  pump  wave,  the  gain  coefficient  is 


proportional  to  r*1  alone.  Carman  et  al.1*  refer 
to  this  as  “the  rather  startling  conclusion. . .  that 
the  Stokes  gain  is  independent  of  the  frequency 
spectrum  of  the  (pump)  laser. . .  even  if  this 
spectrum  is  much  broader  than  T.”  Thus  when 
Tt  is  much  larger  than  r  the  forward  gain  is  much 
larger  than  the  backward  gain.  These  results  go 
against  intuition  based  on  the  idea  that  gain  should 
depend  on  the  number  of  photons  per  unit  frequency 
in  the  pump  beam.  Apparently,  the  concept  of 
photons  as  independent  incoherent  bundles  is  in¬ 
adequate  to  describe  the  subtleties  in  the  SRS 
problem. 

The  purpose  of  this  paper  is  to  further  develop 
the  theory  of  SRS,  including  the  effects  of  finite 
laser  bandwidths,  in  a  way  that  allows  explicit 
calculation  of  the  gain  and  spectrum  of  the  Stokes 
wave.  We  consider,  as  in  Fig.  1,  a  medium  of 
three-level  atomB  interacting  with  two  classical 


1 

T 


FIQ.  1.  Three-level  atom  Interacting  with  a  pump 
laser  with  frequency  uL  and  a  probe  (Stokes)  laser  with 
frequency  ut .  The  cumulative  detunings  are  At  and  Ax. 
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electromagnetic  waves;  one  a  pump  (laser)  wave 
and  the  other  a  Stokes  wave,  differing  in  frequency 
by  the  Raman  shift  of  the  medium.  Both  waves 
have  constant  amplitudes,  but  phases  which  fluc¬ 
tuate  randomly,  giving  rise  to  bandwidths  rL  and 
rs.  This  is  the  “phase-diffusion  model,”  which 
has  been  used  recently  as  a  means  to  include 
bandwidth  effects  into  the  calculation  of  light¬ 
scattering  spectra,  2‘*  as  well  as  multiphoton 
ionization.10’"  The  calculations  are  performed 
as  statistical  averages  over  the  random  phase 
variables  of  the  two  waves.  Previous  treat¬ 
ments18’19  of  SRS  attempted  to  accommodate  both 
fluctuating  phases  and  amplitudes.  We  will  discuss 
later  why  the  inclusion  of  amplitude  fluctuations 
in  this  problem  is  such  a  difficult  task.  Within 
the  stated  model,  we  verify  the  independence  of 
the  forward  gain  from  the  pump  width  rt  and  the 
Stokes  width  Ts ,  in  the  high-gain  limit.  Further 
we  will  show  that  in  the  case  that  the  width  of  the 
pump  laser  rt  is  broader  than  the  Raman  line- 
width  T,  the  amplified  Stokes  signal  assumes  the 
spectral  width  of  the  pump  laser,  regardless  of  its 
initial  width. 

Akhmanov,  D’yakov  and  Pavlov1*  have  separated 
the  problem  into  four  regimes  of  interest  (i)  rt 
«T,  with  no  dispersion;  (ii)  rL«r,  with  dis¬ 
persion;  (iii)  rt»r,  with  no  dispersion;  (iv)  Tt 
»T,  with  dispersion.  Case  (i)  was  considered 
by  Bloembergen  and  Shen, 20  who  predicted  an 
enhancement  of  the  forward  gain  for  a  multimode 
laser.  In  this  paper  we  treat  mainly  case  (iii), 
where  the  laser  linewidth  is  broader  than  the 
atomic  linewidth.  Here  there  is  no  enhancement, 
but  neither  is  there  a  significant  suppression  of 
the  gain,  compared  to  that  calculated  in  case  (i) 
in  the  limit  rt  =  0.  Carman  et  al .‘*  have  also 
treated  case  (iii),  and  they  reached  essentially 
the  same  conclusions  by  calculating  numerical 
solutions  to  the  problem.  Cases  (ii)  and  (iv)  treat 
the  effects  of  dispersion.  There  is  a  consensus18’ 
19,21  that  broadening  of  the  laser  in  the  presence 
of  dispersion  does  result  in  a  lowering  of  the  gain, 
due  to  the  inability  of  the  Stokes  wave  to  stay 
correlated  with  the  pump  fluctuations  as  they  pro¬ 
pagate.  However,  Akhmanov  et  al ,19  have  shown, 
further,  that  there  is  a  critical  pump  intensity, 
above  which  the  effects  of  dispersion  are  overcome 
and  the  gain  coefficient  increases  again  to  nearly 
the  narrow-band  value. 

Dzhotyan  et  al.**  have  treated  the  problem  by 
assuming  the  pump  and  Stokes  waves  to  be  com¬ 
posed  of  many  monochromatic  modes,  with  uniform 
frequency  spacings  large  compared  to  the  Raman 
linewidth  T.  This  results  in  significant  interaction 
only  between  certain  resonant  pairs  of  modes  (one 
pump  and  one  Stokes).  This  approach  can  be 


thought  of  as  complementary  to  the  present  ap¬ 
proach,  in  which  the  energy  in  the  waves  is  taken 
to  be  spread  continuously  over  a  small  frequency 
interval.  The  multimode  approach  of  Dzhotyan 
et  al.  is  a  generalization  of  an  idea  developed  by 
Giordmaine  and  Kaiser22  (and  discussed  by  Byer 
and  Herbst22),  in  which  the  pump  and  Stokes  waves 
each  consist  of  two  modes.  This  treatment  illus¬ 
trates  the  relationship  of  SRS  with  four-wave 
parametric  interactions.  Another  related  dis¬ 
cussion  is  that  of  Harris, 22  in  which  the  threshold 
for  parametric  oscillation  with  multimode  lasers 
is  shown  to  depend  only  on  the  total  power  in  the 
pump  laser. 

In  Sec.  II  we  derive  the  equations  of  motion  for 
the  Raman  problem  in  a  novel  way  by  using  the 
“two- photon  vector  model”  of  Takatsuji24  and 
Grischkowsky  et  al.Zi  In  Sec.  III  we  review  the 
general  solutions  of  the  equations,  following 
Carman  et  al. 18  and  evaluate  the  gain  with  mono¬ 
chromatic  pumping  for  both  the  transient  and 
steady-state  limits.  Then  we  apply  the  phase- 
diffusion  model  to  evaluate  the  gain  under  ar¬ 
bitrary  broad-band  pumping  conditions,  again  in 
both  the  transient  and  steady-state  limits.  In 
Sec.  IV  we  calculate  the  spectrum  of  the  amplified 
Stokes  wave  by  considering  its  autocorrelation 
function  in  the  steady-state  limit.  In  Sec.  V  we 
discuss  a  possible  interpretation  of  the  results 
obtained,  and  in  Sec.  VI  we  summarize  the  main 
results  of  the  paper. 


II.  EQUATIONS  OF  MOTION 


Here  we  give  a  novel  derivation  of  the  usual 
equations  of  motion  for  the  Raman  problem  and  a 
discussion  of  the  physical  model  leading  to  them. 
For  simplicity,  we  treat  the  case  of  near- reson¬ 
ance  Raman  Stokes  scattering,  in  which  only  three 
atomic  levels  need  to  be  considered.  Thus  we  con¬ 
sider  a  vapor  of  atoms  with  energy  levels  shown 
in  Fig.  1.  A  pump  laser  is  tuned  near  (but  not  on) 
the  1-2  transition  and  a  probe  laser  is  tuned  near 
(but  not  on)  the  2-3  transition.  It  is  sufficient  to 
treat  the  pump  laser  as  a  prescribed  field  as  long 
as  it  is  not  depleted.  The  probe  laser  will  ex¬ 
perience  gain  in  a  manner  dependent  on  both  the 
amplitude  and  phase  structure  of  the  pump  laser. 

Consider  the  fields  Et  (pump  laser)  and  E( 
(Stokes,  or  probe  laser)  acting  on  the  three-level 
atom  of  Fig.  1.  Let 

Et  =  S  LStL  cos(u >Lt  -  ktt  *<PL)^S  cos0t ,  (la) 
Es  =SsSs  cos(u>sfTfes*  +  y>,)« <Ssi*  cos$t ,  (lb) 
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where  6 1  and  Ss  are  the  (real)  amplitudes  of  the 
waves,  with  linear  polarization  vectors  $L  andj?s, 
carrier  frequencies  and  ws,  propagation 
vectors  kt  =  kLi,  lts  = ksi ,  and  slowly  varying 
phases  <fit  and  <ps.  The  state  of  the  atom  can  be 
written 

l  +  a,*"*"  i  ***1*1/15 

+«J<r«-t***i-s*^>  (2) 

where  av  a„  and  a,  are  the  slowly  varying  co¬ 
efficients  in  a  “doubly  rotating  frame,"”  and  ipt, 
4>2,  and  <p,  are  the  stationary  eigenstates  of  the 
atomic  Hamiltonian,  with  energies  Aw,,  Rw2,  and 
Aw2.  Atf  =  0  the  atom  is  in  the  ground  state  (a, 

=  1,  a,  =  a3  =  0)  and  afterwards  the  coefficients 
evolve  according  to  Schrodinger’s  equation 


c? 

1 

II 

••a 

(3a) 

*aa  =  (\  -  4>L)a2  -  £os<*j . 

(3b) 

ia j  =  (As  -  vL  +  i>$  )«3  -  . 

(3c) 

where  the  detunings  are  At  =  w2,  - 

and  As  =  Wj, 

+  ws  -  wL.  The  Rabi  frequencies  for  the  two 
transitions  are  given  by  (lL  -  dl2S  t/A  and  fis 
=  dnSs/K,  where<i12  =  dI2  -Eiandd21  =  d23-/rs  aredi- 
pole  matrix  elements.  The  rotating- wave  approxi¬ 
mation”  (RWA)  has  been  invoked  in  writing  Eq,  (3). 
This  is  valid  when  the  detunings  are  small  enough 
(*1**21'  *s  ***,)■ 

Several  authors”'”  have  discussed  the  case  in 
which  Level  2  may  be  eliminated  from  the  Eqs. 
(3a)-(3c).  When  At  is  much  larger  than  As  and 
the  fields  have  no  appreciable  Fourier  components 
at  the  atomic  frequencies,  we  may  set  a2»  0  in 
Eq.  (3b)  and  get  (neglecting  <Pt) 

da**  +  sfija3)/^t  ■  (4) 


This  approximation  is  the  basis  of  the  “two- photon 
vector  model”  of  Takatsuji”  and  Grischkowsky 
ef  of.,25  and  is  discussed  more  fully  in  Appendix 
A.  Using  this  approximation  in  the  Schrodinger 
equation  [Eq.  (3)],  one  may  obtain  two  equations 
for  o,  and  a3,  which  are  identical  in  form  to  those 
of  a  one-photon  transition  with  effective  Rabi  fre¬ 
quency  ft,  =  |ntI2s/At  and  effective  detuning  \ 

=  A,  +J(fit2-ns2)/At,  which  shows  the  effect  of 
ac  Stark  shifting.  We  write  the  resulting  equations 
in  the  convenient  Bloch  form, 23  using  (J  +  iV  =  20,83*, 
and  W=a3aJ-ala*, 

f/  =  _(A#-tf>t  +  £s)F-ru,  (5a) 

V=(*.-i>L*i>a)U  +  n,w-rv ,  (5b) 

w=-n,v.  (5c) 

Here  we  have  included  the  phenomenological  col- 


lisional  dephasing  rate  T ,  which  is  the  halfwidth 
at  half  maximum  (HWHM)  in  rad  sec'1  for  the 
(Raman)  transition  between  levels  1  and  3. 

To  describe  Raman  amplification,  one  must 
solve  the  wave  equation  for  the  Stokes  wave 


8%  1  »2ES  4 it  82P 

8z*  “v2  at1  V  at1’ 


(6) 


where  v  is  the  velocity  in  the  medium,  and  P  is 
the  polarization  of  the  medium.  Considering  only 
Fourier  components  close  to  the  carrier  frequency 
ws ,  and  only  the  linear  polarization  xs ,  leads  to 
a  polarization 

P=N  (4i|d-Js \if>)  =  2Nd2} Re(o2a*) cos0s 


+  2Ndu  Im(a2oJ)  sin0s  , 


(7) 


where  N  is  the  atomic  number  density.  Making 
use  of  Eq.  (7)  for  the  polarization,  Eq.  (4)  to  again 
eliminate  o2,  and  the  slowly  varying  envelope 
approximation,  ”  one  can  write  Eq.  (6)  in  the  form 


3<Ss  18 6s 

*  8*  +v  at 


-k2SlV, 


(8a) 


*■  '■‘•‘'-sS**1 w*  i) . 


(8b) 


where  k2  =  nN  wsvdI2d22/ czK&L  and  the  plus  and 
minus  signs  are  for  copropagating  and  counter- 
propagating  pump  and  probe  beams,  respectively. 
Equations  (5)  and  (8),  along  with  a  similar  one  for 
the  pump  laser  St,  completely  describe  the  pro¬ 
pagation  and  material  response  for  the  Raman 
problem,  including  the  effects  of  phase  modulation 
(<PL,  <ps ).  They  have  been  derived  by  Takatsuji2* 
and  by  Courtens28  in  essentially  the  same  form  and 
used  for  considerations  of  optical  transient 
phenomena. 

In  contrast,  we  are  interested  in  the  special 
case  that  the  atoms  are  weakly  excited  and  the 
pump  laser  is  a  prescribed,  external  field.  Thus 
we  take  for  the  inversion  W=-l,  W=0  in  Eq.  (5). 
W.e  also  assume  exact  resonance  (A,  =  0).  The 
remaining  four  equations,  (5a),  (5b),  (8a),  and 
(8b),  can  be  combined  into  two  complex  equations. 
Defining  the  quantities 


Et  =Ste**i;  Es  =  Sse**s; 

9  =  (U  +  fV>M*i***\  (9) 

we  obtain  directly 

(i0a) 


8  Q* 

at 


-r<J* +  *«,£££ 


5  • 


(10b) 


where  «,  =d,2d21/ 2R2A1.  These  two  coupled  equa¬ 
tions  are  the  starting  point  for  many  theories  of 
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stimulated  Raman  scattering.  They  have  usually 
been  derived  in  the  coupled  wave  approach  of  non¬ 
linear  optics”  for  the  case  of  molecular  Raman 
scattering,  where  Q  is  the  normal-mode  coordinate 
of  a  molecular  vibration  and  is  often  called  an 
optical- phonon  wave.  In  those  treatments,  per¬ 
turbation  theory  was  used  from  the  beginning  and 
the  coupling  constants  *,  and  k2  were  given  in 
terms  of  molecular  polarizabilities.  Here  we  have 
provided  a  connection  between  the  “two- photon 
vector  model"  and  the  standard  theories  of  Raman 
propagation.  We  have  given  the  explicit  relations, 
Eq.  (9),  between  the  variables  used  in  the  earlier 
nonlinear-optics  theories  and  the  more  modern 
optical  resonance  or  Bloch  vector  picture,  which 
has  been  used  here,  and  continues  to  give  insight 
into  many  laser-related  problems. 

III.  EVALUATION  OF  RAMAN  GAIN 

General  solutions  of  Eq.  (10)  have  been  obtained 
in  the  case  of  copropagating  waves  by  Carman  . 
et  al .‘8  In  this  case  the  prescribed  undepleted 
pump-laser  field  EL  depends  only  on  the  local 
time  variable  r  =  t  -  z/v.  It  is  assumed  that  the 
waves  travel  with  equal  velocity  (kt  =  ks),i.e., 
there  is  no  dispersion.  Denoting  by  £s{ 0,r)  the 
Stokes  field  at  the  input  of  the  cell  (z  =  0),  the 
solution  for  the  Stokes  output  field  is1* 


approximations  and  numerical  evaluations  of  the 
Stokes  output,  given  by  Eq.  (12t.  for  two  different 
limits. 


I.  Transient  limit 


The  transient  limit  occurs  for  times  much  short¬ 
er  than  the  reciprocal  of  the  Raman  linewidth 
(I't«  1).  After  replacing  e~T'  by  1  in  Eq.  (12),  the 
integral  can  be  done  to  give£s  (* ,  t)  =  <5S04[(a*T)1/’]. 
which  for  large  enough  o izr  leads  to  the  asymptotic 
form  for  the  output  Stokes  intensity 


eiu.-n 


£2  -2(««T)l/2 

2 v  (<w)‘/j 


(high  gain,  rT«l).  (13) 

We  have  used  the  property  f,(x)  — e,/(2nx),/*,  for 
x-«,  for  any  i.2B  Equation  (13)  is  the  usual  re¬ 
sult  for  the  transient  Raman  effect.”  It  is  inter¬ 
esting  to  note  that,  in  the  transient  limit,  the 
Raman  gain  given  by  Eq.  (13)  does  not  depend  on 
the  Raman  linewidth  T. 


2.  Steady-state  limit 


The  steady-state  limit  occurs  for  times  much 
larger  than  the  reciprocal  of  the  Raman  linewidth 
(rr  »  1).  Extending  the  upper  limit  to  infinity, 
the  integral  in  Eq.  (12)  can  be  done  exactly”  to 
give  for  the  Stokes  intensity 


Jrr  --r*r-r  J 

o  [/>Cr)-/>(r')],7i 

x/.tfox.Xjzt^r)-^')]}*'*) 

xEL(T)El(T')Es<fl,T')dr' ,  (11) 

where  lt(x)  is  the  Bessel  function  of  imaginary 
argument,”  and 

P(r)=  f  |Ut(T')|*<fr* 

•'o 

is  the  integrated  power  in  the  pump  laser  up  to 
time  t. 

A.  Stolen  gain  for  monochromatic  pump  and  input  warn 

It  is  instructive  to  review  the  properties  of  the 
solution  (11)  for  the  case  that  the  t>ump  wave  and 
the  input  Stokes  wave  are  constant  and  monochro¬ 
matic.  In  this  case  we  have  £t(T)  =  E£(T')  =  A’t  and 
£*(0,r')  =  $to.  This  leads  to 

(a*)1'* 

xJf^2i,((o«),,‘)«k,  (12) 


Es(z,T)  =  &loe“  (arbitrary  gain,  rr»l),(14a) 
where 


ifl 


(Mb) 


The  steady- state  gain  coefficient/:  is  the  usual 
one  derived  for  stimulated  electronic  Raman 
scattering.”  It  does  depend  on  the  Raman  line- 
width  T.  in  contrast  to  the  transient  case. 

Equations  (13)  and  (14a)  for  the  output  Stokes 
intensity,  along  with  numerical  evaluation  of 
Eq.  (12),  are  plotted  in  Fig.  2,  as  a  function  of 
gz,  or  equivalently,  pump  laser  intensity,  for 
both  a  transient  case  (rr  =  10"*)  and  a  steady-state 
case  (rr  =  10J).  Here  we  interpret  r  as  the  pulse 
duration  of  the  pump  laser.  Equation  (13)  for  the 
transient  gain  (Ft  =  10*2)  is  seen  to  agree  well  with 
the  exact  numerical  results  when  log,„(£j/Aj0) 

>1,  while  Fq.  (14 a)  for  the  steady-state  gain  (rr 
*  102)  agrees  everywhere.  Note  that  the  values 
below  log1()(£j/AJo)  =  0  are  unphysical.  The  other 
four  curves  in  Fig.  2  show  the  effects  of  laser 
bandwidth  on  the  gain,  as  described  in  Sec.  WB. 


R.  Raman  gain  for  broad-band  pump  and/or  input  warns 


where  we  have  used  x  *  r  -  t',  p(t)  =  S\t  and  a.  We  now  evaluate  the  Raman  gain  in  the  case  that 

Following  Wang”  we  present  analytic  the  spectral  width  of  either  the  pump  laser  or  the 
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Stokes  input  wave  (or  both)  is  greater  than  the 
Raman  linewidth  (rt  *-rs  »  T).  This  can  be  ac¬ 
complished  by  performing  an  average  o(  the 
general  solution,  Eq.  (11),  over  a  statistical  en¬ 
semble  chosen  to  model  .he  bandwidths.  An 
especially  useful  model  is  that  of  phase  diffu- 
slon,  ’**  in  which  the  field  amplitude  6  L  is  constant, 
but  the  phase  suffers  abrupt  changes  at  an  average 
rate  2Et  (see  Appendix  B).  The  field  autocorrela¬ 
tion  functions  are  then 


«£t(T)£2(r,)»^ie-ri'f-r' 

U5a) 

«£*  (0,r)£j(0,  r')»  =  , 

(15b) 

which  lead  directly  to  Lorentzian  line  shapes  with 
halfwidths  (HWHM)  equal  to  rt  for  the  pump  laser 
and  Tj  for  the  input  Stokes  wave.  The  brackets  ((  )) 
indicate  an  ensemble  average  over  the  statistical 
fluctuations  of  the  field.11  This  model  describes 
a  stabilized  laser  operating  far  above  threshold, 
but  it  also  proves  to  be  very  convenient  mathe¬ 
matically. 

The  intensity  of  the  output  Stokes  wave  is  given 
by  ((|£s Cz,t)|2)).  To  evaluate  the  intensity  of  the 
Stokes  wave  we  first  introduce  some  notation: 


/W- (#*•//*  )/,[(«*)*'*],  (16a) 

6,0  (18b) 

Then  from  Eq.  (11),  using  p(r)  =  S\r  in  accordance 
with  the  phase  diffusion  model  for  the  pump  laser, 
we  have 

Es  (*,  r)  =  Es  (0,  r)  +  5S0{i[(az)1/2]}E(r) ,  (17a) 

«|«sU,t)|  »»  =  510  + 5, „(«)*/* 
x<(£J(0,t)F(t)» 
>«!ofi(«*)K<|F(T)|I».  (17b) 

The  second  term  in  Eq.  (17b)  can  be  easily 
evaluated  using  Eq.  (IS), 


«EJ(0,t)F(t)»=$Jo  f/(T- r') 

•'o 

x  »e-ri «*•*•> dr' ,  (18) 

where  we  have  assumed  statistical  independent  of 
the  waves  at  z  =  0;  that  is 

«Et(T)E£(T')£,(0,T')E|(0,T)» 

*  «£t(T)££(r')»«£,  (0,  r')£|(  0,  t)» 

(see  Ref.  32).  This  integral  is  identical  to  the 
Integral  in  Eq.  (12),  but  with  r  replaced  by  r  +rt 

♦  Tj,  Thus  when  Tj  +  Tj  »T,  this  term  grows 
with  a  very  small  gain  coefficient.  This  is  in 


FIG.  2.  Normalized  Stokes  output  intensity  as  a  func¬ 
tion  of  gain  coefficient  gz  (or  equivalently,  pump  laser 
intensity)  under  various  physical  conditions  and  differ¬ 
ing  levels  of  approximation.  The  curve  labeled  “NB.  SS” 
is  the  narrow-band  steady-state  result  obtained  from 
Eq.  (14a),  or  Eq.  (12)  with  Tt-102,  where  T  is  the 
Raman  linewidth  and  t  is  the  laser  pulse  length.  The 
“NB.TB”  curves  are  the  narrow-band  transient  results, 
obtained  exactly  from  Eq.  (12)  with  Tt  =  1(T 2  (solid 
curve),  or  approximately  from  Eq.  (13)  (dashed  curve). 
The  curves  labeled  "BB,  SS”are  the  broad-band  steady- 
utate  results  obtained  exactly  from  Eqs.  (17b),  (18), 
and  (20b)  with  Tt«102  (solid  curve),  or  approximately 
from  Eq,  (25)  (dashed  curve).  The  “BB.TR”  curves 
are  the  h road-band  transient  results  obtained  exactly 
from  Eqs.  (17b),  (18),  and  (20b)  with  Tr»  10* 5  (solid 
curve),  or  approximately  from  Eq.  (22)  (dashed  curve). 
The  broad-band  curves  are  for  a  bandwidth  ratio  (rt 
♦  r#)/r-108,  where  rt  and  r5  are  the  bandwidths  of 
the  pump  and  probe  (Stokes)  lasers. 


contrast  to  the  third  term  in  Eq.  (17b),  which,  as 
we  will  see,  grows  with  essentially  the  narrowband 
gain  given  in  Eq.  (14a).  Thus  we  expect  the  third 
term  of  Eq.  (17b)  to  be  dominant  in  the  high-gain 
limit  Using  the  correlation  functions  given  by 
Eq.  (15),  one  can  evaluate  thiB  third  term  as 
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« |F(t, |>»  =  f  rAT  .  ry(r  -  r>)£|(Q,  r-)))^, 

o  */0  &L  O$o 

*  f  f  f(r  - 1  dr'dr*  f 

•'o  Jb 


where  we  used  the  fact  that  £t(T)£2(T)=£*  is  in¬ 
dependent  of  the  statistical  averaging.  To  evaluate* 
the  double  integral  in  Eq.  (19b),  we  note  that  the 
exponential  factor  is  much  different  from  zero 
only  near  the  tine  r'  =  r*.  In  the  limit  that  /(x) 
changes  very  slowly  in  a  time  (rL  +  rs)*‘  (i.e., 
®2<<I'z  +  rs),  we  can  replace  the  exponential 
factor  by  the  properly  normalized  6  function 
2(r x,  +  rs  r*  6(t'  -  t').  We  then  get 

«|^(T)|2»=~-j;-^/*(T-r')/ir'  (20a) 

"r~!rs  I0  (20b) 

In  similar  fashion  to  the  integral  in  Eq.  (12),  this 
integral  can  be  evaluated  analytically  as  well  as 
numerically,  in  the  two  limits: 

1.  Trvaknt  limit 

As  before,  for  Tr  « 1,  the  exponential  can  be 
replaced  by  1  and  the  integral  done  (in  this  case 
asymptotically,  using  the  asymptotic  form  for 
I,)  to  give 


«|Jr(T)|a»  = 


®(rz+r4) 


Thus,  in  the  high-gain  limit  where  Eq.  (17b)  iB 
dominated  by  the  last  term,  we  find  that  the  Stokes 
output  intensity  in  the  transient  limit  is 

g2  02<a^t),,2 

((|£s(z,T)|*»=^a  2(rt+rf) 

(high  gain,  rT«l).  (22) 

Because  of  the  form  of  the  exponential,  this  result 
for  the  broad-band  transient  will  be  nearly  indis¬ 
tinguishable  from  the  result,  Eq.  (13),  for  the 
narrowband  transient 

2.  Strmty-tUUt  limit 

To  evaluate  the  steady-state  limit  of  Eq.  (20b), 
we  extend  the  upper  limit  of  the  integration  to 
infinity,  and  do  the  integral  to  give11 

S»> 

where  ,F,  is  the  generalized  hypergeometrlc 
function,  which  can  be  evaluated  asymptotically 
for  large  argument  as** 

»F,(1,  l;2,3,x)-(4/V"i)  e*/x*/2.  (24) 


This  leads  to  an  asymptotic  form  for  the  Stokes 
intensity  under  broad-band  steady-state  high-gain 
conditions: 

(high  gain,  Tt  » 1) .  (25) 

It  can  be  seen  by  comparing  the  Stokes  intensities 
given  by  Eqs.  (25)  and  (14a)  that  under  the  condi¬ 
tions  assumed  (copropagating  waves,  no  disper¬ 
sion),  the  growth  of  a  Stokes  wave,  in  the  steady- 
state  high-gain  limit,  is  virtually  unaffected  by 
either  its  input  bandwidth  or  that  of  the  pump 
laser.  We  can  demonstrate  this  result  by  writing 
the  output  Stokes  intensity  as 

«|£sU,,<>)|2S»  =  ilioe6.  (26) 

G=GHB_ln{I(ri+rs)/r](^I1B)*">,  (27) 

where  GKB=gz  is  the  narrow-band  gain  coefficient 
from  Eq.  (14b).  Thus  for  large  gain  the  difference 
between  G  and  GNB  becomes  relatively  insignifi¬ 
cant.  We  will  present  a  possible  interpretation  for 
this  result  in  Sec.  V.  Equation  (27)  is  similar  to 
the  result,  conjectured  by  Carman  et  al., 11  that 
G=GMB-ln{(ri/r)GNBl.  The  difference  between 
our  result  and  theirs  (when  Ts  =  0)  may  be  due  to 
the  fact  that  they  allowed  also  for  amplitude 
fluctuations  of  the  pump  laser,  whereas  we  have 
restricted  ourselves  to  phase  modulation  alone, 
in  order  to  make  an  explicit  calculation  tractible. 
The  calculation  becomes  intractible  when  ampli¬ 
tude  fluctuations  are  present  because,  then,  p{r) 
in  Eq.  (11)  is  a  random  variable  which  makes  the 
statistical  average  difficult  to  perform. 

As  in  the  narrow-band  case,  we  plot  the  broad¬ 
band  Stokes  intensities,  Eqs.  (22)  and  (25),  in 
Fig.  2,  along  with  numerical  evaluation  of  the 
Stokes  intensity  from  ((|£s(*,t)|3»,  defined  by 
Eqs.  (17b),  (18),  and  (20b).  Since  Eq.  (20b)  is 
valid  only  for  rt  ♦  rs  »  r,  az,  we  have  plotted  the 
extreme  case  that  (r  t  +  r4)/T  =  105,  in  order  to 
demonstrate  the  validity  of  the  asymptotic  forms 
Eqs.  (22)  and  (25).  Again  we  have  plotted  the 
transient  and  steady-state  cases:  rr  =  lCT*  and 
103.  Again  we  see  agreement  of  the  asymptotic 
forms  with  the  exact  numerical  results  when 
log«!£s(z,T)|*»/$J#>|.. 

We  may  now  compare  the  narrow-band  and 
broad-band  results.  For  this  extreme  steady-state 
easel  (rt+r,)/r  =  10i,  we  see  a  Significant  sup- 
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pression  of  the  gain  in  the  turn-on  region  for  both 
transient  and  steady- state  limits.  However, 
according  to  Eq.  (27),  at  very  high  gains  the  dif¬ 
ference  between  the  narrow-  and  broad-band 
stimulated  outputs  becomes  less  and  less,  relative 
to  their  absolute  magnitudes.  It  is  interesting  that 
the  broad-band  output  exhibits  a  threshold- type 
behavior,  in  contrast  to  the  exponential  behavior 
of  the  narrow-band  output. 

Although  the  principal  interest  here  is  in  the 
high-gain  limit,  some  comments  can  also  be  made 
about  the  low-gain  steady-state  limit,  important 
to  such  experimental  techniques  as  CARS  (coher¬ 
ent  anti-Stokes  Raman  scattering).  Thus  we  will 
find  the  terms  in  Eq.  (17b)  which  are  of  lowest 
order  (linear)  in  az.  It  can  be  shown34  that  the 
third  term,  containing  ^|£(t)|j^,  is  quadratic 
in  at  as  at  —  0  and  can  thus  be  neglected.  The 
first  and  second  terms  in  Eq.  (17b)  are  evaluated, 
using  Eq.  (18)  with  the  upper  limit  r  taken  to  in¬ 
finity,  to  give29 

«|£4(z,0O)|2»  =51o(2«'»s*/,-1),  (28) 

*„  =  a/2(r+rt+rs).  (29) 

In  the  limit  at  —  0,  this  reduces  to 

(low  gain,  Tr  » 1) .  (30) 

We  see  that  in  the  low-gain  steady-state  limit, 
the  SRS  grows  linearly  with  the  (“broad-band”) 
gain  coefficient  gBB.  This  is  the  result  that  would 
naively  be  predicted  on  the  basis  of  photons  per 
mode,  as  discussed  in  Sec.  L 

C.  Raman  gain  for  arbitrary  bandwidths 

Here  we  analyze,  numerically,  the  properties 
of  the  stimulated  output  when  the  condition  r  L  *  Ts 
»  r,  at  is  not  necessarily  upheld,  as  was  assumed 
in  Sec.  Ill  B.  First  note  that  if  we  take  (rt  +  rs)/r 
=  102,  rather  than  105  as  used  in  Fig.  2,  the  analy¬ 
sis  of  Sec.  HI  Bis  valid  only  for  gz<  102,  makingpre- 
diction  of  the  broad-band  transient  above  gt  =  102 
impossible  by  those  methods.  However,  also  note 
that  the  solution  in  Eq.  (19b),  before  approximation 
to  obtain  Eq.  (20),  contains  the  information  we  are 
seeking  in  the  general  case.  Thus  we  evaluated 
Eq.  (19b),  by  a  numerical  method  discussed  in 
Appendix  C,  and  obtained  the  output  Stokes  in¬ 
tensity  «|£,(z,t)|2»  defined  by  Eqs.  (17b),  (18), 
and  (19b).  These  results  are  shown  in  Fig.  3, 
where  we  have  covered  a  large  region  of  the  in¬ 
teresting  parameters:  Tr  and  (rt  +rs)/r  both 
vary  between  10*2  and  102. 

Beginning  with  Fig.  3(a),  we  see  that  the  laser 
bandwidth  has  little  effect  on  the  gain  in  the 
transient  limit  (fr  =  10*2).  This  is  not  surprising, 


as  a  short  laser  pulse  of  duration  t  -  10*2/r  has  a 
spectral  width  of  KPr  without  phase  diffusion 
(r£  =  0).  Thus  we  see  no  effect  of  additional 
broadening  by  phase  diffusion  until  rt  ♦  rs  >  lO2!", 
at  which  point  the  gain  becomes  slightly  depressed. 
Progressing  in  Figs.  3(b)  and  3(c)  to  longer  pulse 
duration  r,  we  see  the  general  result  that  no  effect 
of  phase  diffusion  broadening  is  apparent  until 
rz  +  rs  *  l/T.  The  steady-state  limit  occurs  in 
Fig.  3(d),  where  no  difference  is  seen  between 
Tt  =  10  and  Tt  =  102. 

IV.  SPECTRUM  OF  THE  STOKES  OUTPUT  IN  STEADY  STATE 

In  Secs.  I— III  the  input  Stokes  wave  has  been 
taken  to  have  a  width  rs .  But  because  all  of  the 
power  in  the  broad-band  pump  laser  is  effective 
for  amplifying  the  Stokes  wave,  it  is  interesting 
to  ask  what  becomes  of  the  spectral  distribution 
of  the  Stokes  wave  after  it  has  been  amplified.  In 
this  section  we  calculate  the  spectrum  of  the 
Stokes  wave  in  the  steady-state  high-gain  limit, 
in  two  different  cases:  Tt  =  0  and  »  F. 


'•f  <<!> 


FIG.  3.  Normalized  Stokes  output  Intensity  aa  a  func¬ 
tion  of  gt  (or  pump-laser  Intensity),  evaluated  numerl- 
tcally  using  Eqs.  (17b),  (18),  and  (19b),  for  various 
values  of  the  bandwidth  ratio  ( rt+  r t]/T,  where  r  Is 
the  Raman  linewtdth,  and  TL  and  r$  are  the  bandwidths, 
due  to  phase  diffusion  broadening,  of  the  pump  and  probe 
(Stokes)  lasers.  Four  different  pulse  lengths  r  are 
shown:  (a)  rr»  10*2,  transient  limit;  <b)  I't-IOT1;  (c) 
rr-1;  (d)  rr-10  and  103,  steady-state  limit.  In  all 
cases  the  effect  of  the  laser  bandwidths  becomes  rela¬ 
tively  unimportant  at  high  gains. 
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For  a  stationary  wave  £(r),  the  definition  of  the 
power  spectrum  P(u>)  is 

P(w)=i£**,"K(s)‘£*-  (3U 

where  u  is  the  frequency  as  measured  from  the 
frequency  of  the  carrier  wave  and  P{ o>)  is  the 
Fourier  transform  of  the  electric  field  autocor¬ 
relation  function  )C(s), 

£(s)  =  «£(t)£*(t+*)».  (32) 

It  is  easy  to  show  from  Eq.  (31)  that  P(w)  is  nor¬ 
malized  as  follows: 

f "  P(<i>)d<i>  =X(0)  =  « |  £(r)  |  *» .  (33) 

_ I 


As  an  example,  when  Eq.  (31)  is  used  to  calculate 
the  laser  spectrum  P^w),  from  the  correlation 
function  Eq.  (IS),  one  finds 

PL(*)=££rSl,  (34) 

a  Lorentzian,  as  stated  in  Sec.  III.  We  can  use 
the  general  solution,  Eq.  (11),  to  determine  K (s) 
for  the  Stokes  wave 

£(s)“«!oti(«z)K(P(r)F*(T  +  s)» ,  (35) 

where  we  have  again  kept  only  the  term  which 
dominates  in  the  high-gain  limit.  In  steady  state 
(r  -.«)  we  expect  that  K(s)  will  depend  only  on  s 
and  not  r.  K(s)  can  be  evaluated  as 


*  •'o  J0 

r dy/(xr(y)G(x,y,s)e■^*,-^,, 


where 

G(*,y,s)  =  exp[rt(|s  +  *|  +  [s-y|-  |s-y+*{ 


i(j)  -f  (2f  r  +  )  tu)+rt 

A *  (2rt+rs  )4*u*+rf  * 


(36a) 

(36b) 

(38a) 


-|*|-|y|-|s|)].  (36c) 

Here  G(*,y,s)  is  the  four-time  correlation  func¬ 
tion  of  the  pump-laser  field,  assuming  the  phase 
diffusion  model,  and  is  evaluated  in  Appendix  B. 

In  deriving  Eq.  (36b)  we  have  used  x  =  T  -  T'  and 
y  =  r  +  s  -  t'.  In  steady  state  the  upper  integration 
limits  are  extended  to  infinity  and  K (s )  becomes 
independent  of  r.  In  order  to  simplify  the  absolute 
values,  the  integral  is  transformed  to  the  triangu¬ 
lar  region  above  the  y  =x  line  by  use  of  the  pro¬ 
perty  C(y,*,s)  =  G(*,y,-s).  Then  for  the  Stokes 
spectrum  we  have 


P . 


*  2£ 
s„4- 


we  can  write 

£1  =  2ReA_e-‘ri*r*,,«”» 

+  2  coswfy  ~x) 

-2ImA.e-ri*»-*>sinw(y-x),  (38b) 

£2  =  2(-ReA.  costox  +  ImA_  sinwx) 

-2(ReA.  cosojy  +  ImA.  sinwy ) 
xe-rttr»w>-rt»  (38c) 

Equations  (37)  and  (38)  are  now  used  to  evaluate 
the  Stokes  output  spectrum  in  two  different  cases. 


xjT  dy  <fx/(x)f(y)£(x,y,w),  (37a) 

where 

£(x,y,  <*>)  =  2  Re<rrz'»**»  f“  ds  «-*“• 
x  exp(rt(|s  **|  +  |*  -y  | 

-  I* -y >x| -  |  si)  -r,|s-y+*|]. 

(37b) 

The  transform  £  can  be  calculated  under  the  condi¬ 
tion  x«y.  We  write  £  as  the  sum  of  two  parts 
£=£,+£,.  Defining 


A.  Stokes  output  spectrum  for  s  monochromatic  pump  laser 

Here  we  treat  the  case  that  the  pump  laser  is 
monochromatic  (Tt  =  0)  and  the  spectral  width  of 
the  input  Stokes  is  allowed  to  assume  two  different 
limits:  T,  =  0  or  Ts  »  T.  The  steady-state  gain 
for  these  two  limits  has  already  been  given  in 
Eqs.  (14)  and  (25). 

The  spectrum  is  easily  obtained  by  setting  Tt 
=  0  in  Eq.  (38).  Then,  because  A.  =  0,  we  have 
£,  =  0  and 

£,  =  4[r,/(<^  ♦  r* )]  cosw(y  -  *) .  (39) 

Transforming  back  to  the  full  x,y  quadrant  gives 
for  the  spectrum  of  the  Stokes  output 
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„  ,  ,  i  n  4rs 

I 

f  dy  ]  dxf(x)f{y)  cosu>(y  -x)  (40a) 

*  ■'o  ^0 


l/l'-  <** 

where” 

/=  /  e,u'f(x)dx 

JQ 

’(^WfrT^)"1]-  (40c 

In  the  high-gain  limit  the  Stokes  output  spectrum 

is 

PsM-SIo-^yt^p  <41 


.  B.  Si  oka  output  ipectrum  for  a  broadband  pump  laser 

Here  we  treat  the  case  that  the  pump  laser 
width  is  larger  than  the  atomic  width  (Tt »  r) 
and  the  input  Stokes  width  rs  is  arbitrary.  The 
spectrum  is  obtained  by  applying  several  approxi¬ 
mations  to  Eq.  (38).  First  note  that  in  the  high- 
gain  limit,  only  X,  will  give  a  significant  con¬ 
tribution  to  the  spectrum  because  its  exponentials 
damp  as  ( y -x),  rather  than  (>*  +x).  Thus  along 
the  line  y  =  x,  X,  is  large,  while  X2  becomes  neg¬ 
ligible.  Second,  note  that  because  I\  +  rs  is 
assumed  large  (rt  +  Ts  »  T,  az),  we  may  replace 
the  exponential  factors  by  properly  normalized 
delta  functions,  as  in  the  argument  leading  to  Eq. 
(20); 

e‘  Tz*r  s  >  -  (rt  +  rs )- « 6(  y  _  x) ,  (44a) 

e'rt‘»*,,COsw(y  -x)-[rt/(«a  +  r|))6(y  -x), 


To  put  this  result  into  proper  form  we  must  as¬ 
sume  one  of  the  two  above-mentioned  limits. 

The  first  limit  is  that  of  a  monochromatic  input 
Stokes  wave  (rs  =  0).  Here  we  can  use  Eq.  (14a) 
for  the  output  Stokes  intensity:  (( |£s  (z, 00 )  |a» 
=5§0e".  Then  Eq.  (41)  can  be  rewritten,  in  the 
limit  r4-0, 

Ps(a,)  =  6M«|Es(z,«)|2».  (42) 

Thus  we  see  that  the  output  Stokes  wave  is  mono¬ 
chromatic  when  the  input  Stokes  wave  and  the  pump 
laser  are  monochromatic,  as  expected. 

The  second  limit  is  rs  »  r.  Here  we  again 
consider  a  monochromatic  pump  laser  (rt  =  0) 
and  use  Eq.  (25)  for  the  Stokes  output  <(|£s(z,  ®)|*> 
to  rewrite  Eq.  (41)  for  the  Stokes  output  spectrum 
as 


x«|£,(e,-)|*>,  (43a) 

*«(">*  (tf)  exp(~  r*"’ w’)  *,)fa * » 

(43b) 

where  the  last  step  is  valid  because  the  halfwidth 
[Cln of  the  exponential  factor  is  much 
less  than  T  and  rs.  Note  that  Eq.  (43b)  is  nor¬ 
malized  as  in  Eq.  (33).  Equation  (43)  describes 
a  Lorentzian-shaped  atomic  line  of  width  T  that 
has  been  gain  narrowed.  The  ideal  that  the  center 
of  the  line  will  experience  more  gain  than  the 
wings  is  a  familiar  idea  in  laser  theory.  A  com¬ 
parison  of  the  atomic  line  and  the  gain-narrowed 
Stokes  line  is  shown  in  Fig.  . 4  for  the  case  gz  « 15. 


(44b) 

e-rt  if* » sln  ^  _  x)  _  [u/(uj7  +  )]  6(y  -  x) . 

(44c) 


-2r  -r  o  r  2r 


FREQUENCY,  w 

FIG.  4.  Comparison  of  the  Lorentzlan  Raman  line 
shape  (curve  (a))  with  atomic  halfwtdth  r,  and  the  gain- 
narrowed  output  Stokes  spectrum  (curve  (b))  with 
halfwtdth  (ln2r*/g*),/1  for  gt « 15,  plotted  from  Eq. 
(43b).  This  Stokes  spectrum  narrowing  results  when  the 
pump  laser  is  monochromatic  and  the  input  Stokes  laser 
is  broad  band. 
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Combining  these  results  gives 

£,-[4/(ri  +  rs )](  rt/(«» + ri)]e(y  -  *) .  (45) 

Now  the  spectrum  of  the  output  Stokes  intensity 
Ps  (<*>)  is  easily  evaluated  from  Eq.  (37a): 


/*,(»> 


PM 


at 

T 


* 

rt+rs  STTrJ 


x/’/Wdy, 

J  0 

-)!*»• 


(46a) 

(46b) 


where  we  used  Eqs.  (17b)  and  (20a)  to  define  the 
output  Stokes  intensity  i  j£s(z,«°)|JJ.  Note  that 
Ps(v)  is  normalized  as  required  by  Eq.  (33).  This 
indicates  'hat  our  neglect  of  £2  is  justified. 

Comparing  Eqs.  (46b)  and  (34)  shows  that  the 
Stokes  wave  assumes  exactly  the  same  spectrum 
as  the  pump  laser  during  the  amplification  process, 
regardless  of  the  spectral  width  of  the  input  Stokes 
wave. 


V.  DISCUSSION 

In  the  case  just  treated,  that  the  pump  laser 
width  is  greater  than  the  atomic  width  (rt»  r), 
our  interpretation  is  that  the  fluctuations  in  the 
phase  of  the  pump  laser  dominate  the  time  be¬ 
havior  of  the  amplification  process.  Indeed, 

Carman  et  al.1*  found  numerically,  in  the  case  Of 
a  quadratic  phase  sweep  in  the  pump  laser,  that 
the  Stokes  phase  closely  followed  this  sweep  after 
a  brief  initial  period.  Our  result  for  the  Stokes 
spectrum  (Eq.  (46)]  is  consistent  with  the  conjecture 
of  Carman  et  al.  that,  when  rt »  r,  the  Stokes 
phase  always  follows  the  pump  phase  in  the  high- 
gain  limit,  regardless  of  the  phase  structure  of 
the  input  Stokes  wave.  If  correct,  this  effect, 
which  we  will  call  "phase  locking,”  also  explains 
the  fact  that  the  gain  is  unaffected  by  the  phase 
flucutations  which  lead  to  the  bandwidth.  As  Car¬ 
man  et  al.  pointed  out,  if  the  phases  <ps  and  <pt 
differ  at  all  points  by  a  constant  [<ps(t,  t )  =  <pL(t,  0 
♦  <p0],  the  phases  drop  out  entirely  from  Eq.  (11), 
leading  to  the  narrow-band  gain  result,  Eq.  (14). 
Thus,  the  idea  of  "phase  locking”  leads  to  results 
consistent  with  our  results  for  VL  »  T.  When  both 
rt  and  rs  are  larger  than  r,  we  can  say  that  the 
amplified  Stokes  wave  builds  up  from  the  broad¬ 
band  input  noise  in  a  way  which  automatically 
satisfies  <ps  *<pt  ♦  <p0-  That  is,  only  that  part  of 
the  noise  which  satisfies  this  relation  will  ex¬ 
perience  large  gain. 

To  illustrate  the  idea  of  phase  locking  we  have 
compared,  in  Fig.  5,  several  steady-state  gain 
curves.  We  have  reproduced  curves  from  Fig. 


3(d),  calculated  from  the  exact  equations  (labeled 
‘phase  locked”).  We  have  also  plotted  curves  using 
Eq.  (14a)  (labeled  “narrow  band”),  and  also  using 
Eq.  (14a)  with#  replaced  by  #BB  of  Eq.  (29)  (labeled 
“unlocked").  We  see  that  at  low  gains  the  exact 
curve  follows  the  “unlocked"  curve,  consistent 
with  the  idea  that  there  is  no  correlation  between 
the  output  Stokes  and  pump  laser  waves.  This 
low-gain  behavior  was  predicted  at  the  end  of 
Sec.  Ill  B.  However,  at  high  gains  the  exact  curve 
approaches  the  “narrow-band”  curves,  consistent 
with  the  idea  that  it  has  become  “phase  locked," 
resulting  in  an  enhanced  gain.  We  thus  see  that 
phase  locking  appears  to  occur  only  above  a  cer¬ 
tain  (threshold)  gain.  In  contrast  to  the  behavior 
found  in  the  present  treatment,  Dzhotyan  et  al.M 
found  in  the  multimode  approach  (see  Sec.  I)  that 
the  “narrow-band”  gain  was  appropriate  even  at 
low  gains.  This  is  a  major  difference  between  the 
two  models. 

Finally,  we  point  out  that  we  have  treated  only 
the  case  of  Raman  amplification,  and  not  SRS 
which  grows  from  the  initial  Stokes  photons  spon¬ 
taneously  emitted  with  frequencies  near  ios,  in 
the  absence  of  an  external  input  Stokes  wave  at 
that  frequency.  Here  we  wish  to  make  some  con¬ 
jectures  on  the  outcome  in  the  latter  case.  We 
may  consider  the  spontaneous  photons  as  making 
up  the  source  term  £s(0,  r).  Although  here  we 
certainly  cannot  make  the  decorrelation  of  the 
pump  wave  E^r)  and  source  term  £s( 0,  r)  that 
we  made  in  connection  with  Eq.  (18),  we  still  ex¬ 
pect  that,  at  high  gains,  the  major  results  we  have 
obtained  do  apply  to  spontaneously  generated  SRS. 
That  is,  we  expect  the  gain  to  be  essentially  in- 


FIG.  S.  Normalized  Stokes  output  intensity,  in 
steady-state,  as  a  function  of  gt  (or  pump-laser  Inten¬ 
sity),  for  two  different  laser  bandwidth  ratios:  (a) 
(rt+rs)/r-l  and  (b)  (rt*rs)/r-10.  In  both  cases, 
the  exact  results  (labeled  “phase  locked")  are  seen  to 
agree,  at  low  gains,  with  the  results  one  would  expect 
In  the  absence  of  phase  locking  (labeled  'Unlocked”), 
which  were  obtained  by  replacing  T  by  T«  I"*  ♦  I*  In  the 
expression  |Eq.  (14b))  for  the  gain  coefficient#.  How¬ 
ever,  at  higher  gains,  the  exact  results  approach  the 
narrow-band  curve  (Eqs.  (14a)  and  (14b)]. 
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a2(f ) ■  f  c’** >£(/') dt' , 
J  0 

*(/).»[  |  nta,(/M  «,«,&)]. 


(Al) 

(A2) 


dependent  of  pump-laser  bandwidth,  and  when  the 
pump  bandwidth  is  greater  than  the  atomic  width 
we  expect  the  output  SRS  to  assume  the  spectrum 
of  the  pump. 

VI.  SUMMARY 

Using  the  phase  diffusion  model,  we  have  ex¬ 
tended  the  theory  of  stimulated  Raman  scattering, 
in  the  case  of  two  interacting  classical  waves 
(pump  and  input  Stokes),  to  allow  for  arbitrary 
bandwidth  of  either  wave.  In  the  forward  direction 
If  there  is  no  dispersion,  we  showed  that,  in  the 
high-gain  limit,  the  gain  of  the  Stokes  wave  is 
essentially  independent  of  the  input  bandwidth  of 
either  wave.  In  the  low-gain  limit  the  gain  co¬ 
efficient  was  found  to  be  inversely  proportional  to 
the  sum  of  the  bandwidths.  We  also  calculated  the 
spectrum  of  the  output  Stokes  wave,  in  the  high- 
gain  limit,  under  various  conditions.  We  found 
that  when  the  pump  bandwidth  r L  is  greater  than 
the  atomic  width  r,  the  Stokes  wave  assumes  ex¬ 
actly  the  spectrum  of  the  pump  laser,  regardless 
of  the  spectral  width  of  the  input  Stokes  wave. 
When  both  input  waves  are  monochromatic 
(rt,  rs  =  0),  we  found  that  the  Stokes  spectrum  is 
unchanged  by  the  amplification  process.  Finally, 
when  =  0  and  Ts  »  T,  we  found  that  the  output 
Stokes  wave  has  a  spectrum  which  is  a  gain- 
narrowed  atomic  profile;  that  is,  the  Stokes  width 
becomes  much  narrower  than  the  atomic  width. 

Note  added  in  proof.  A  recent  preprint  by  W.  R. 
Trutna,  Y.  K.  Park,  and  R.  L.  Byer  [to  appear  in 
IEEE  J.  Quant.  Electron.  (July  1979)]  has  come  to 
our  attention.  Broad-band  SRS  was  treated  using 
the  coupled- wave  approach  (similar  to  that  in  Ref. 
22)  and  qualitative  agreement  was  found  with  our 
work  in  the  high-gain  limit.  At  low  gains,  how¬ 
ever,  their  treatment  indicates  no  suppression  of 
the  gain,  in  contrast  to  our  results  [Eq.  (30)]. 
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APPENDIX  A 

Here  we  discuss  more  carefully  the  elimination 
of  a,  from  Eq.  (3).  We  first  neglect  <pt  In  Eq. 

(3b),  as  At  is  assumed  to  be  much  larger  than 
the  pump  laser  bandwidth.  The  formal  solution 
of  Eq.  (3b)  can  then  be  written 


Repeated  integration  by  parts  gives 

a’(,)  - TTt - {ia 7? - 

*  — - •  i*3) 

Note  that  since  n,(0)  =  1  and  n,(0)  =  0,  c(0)  =  i{\)nL. 
When  one  assumes  »  As,  ^>L,  t>s,  ,  R s,  it  can 
be  shown  from  Eq.  (3)  that£(/)<^  Thus 

when  At  is  large,  one  is  left  with 

«a(f)  =  \g(t)  -  .  (A4) 

However,  because  is  large,  the  exponential 
term  oscillates  rapidly  compared  tojf(f).  Hence, 
in  the  spirit  of  the  RWA,  we  neglect  the  rapidly 
oscillating  part  and  retain  only  the  slowly  varying 
part; 

«2(f)s^(f)/*'Al  =  (50t«I(f)insa,(/))/At .  (AS) 

It  is  interesting  that  the  same  result  is  obtained 
by  merely  setting  a,  *  0  in  Eq.  (3b). 

APPENDIX  B:  FHASE-DIFFUSION  MODEL 

The  phase-diffusion  model  for  laser  bandwidth 
describes,  to  good  approximation,  a  cw  laser 
operating  well  above  threshold,  where  the  in¬ 
tensity,  J(f)=/0+/'(f),  is  nearly  constant,  with 
average  value  J„  and  small  fluctuations  I'Q  ).** 
However,  well  above  threshold  the  phase  <p(t) 
fluctuates  randomly,  in  a  way  reminiscent  of  a 
diffusing  Brownian  particle.  Simple  laser  theory 
gives  the  equations  for  the  intensity  and  phase 
as” 

/»(/)  —  X/'(f)+F/f),  (Bl) 

<J(f)"F,(f),  (b2) 

where  F,(t)  and  F,(f)  are  random  Langevin  forces 
with  correlation  functions 

CF,(f,)F/(ft)»=  2m,  -/,), 

«F,(ft)F.(/,)j  .irMf,  -  ft) , 
and 

«WF.(f,)».  0. 

Here  1 A  is  the  correlation  time  of  the  intensity 
fluctuations,  with  mean  value  D/\,  and  r  is  the 
bandwidth  of  the  lighL  These  6  correlations  sim¬ 
ply  imply  that  the  forces  fluctuate  on  a  time  scale 
shorter  than  any  other  interesting  time  scale. 

The  phase-dlffusioa  model  is  based  on  the  as¬ 
sumption  that  the  intensity  exhibits  no  fluctuations, 
/'(f)  >0,  and  that  the  phase  fluctuates  according 
to  Eq.  (B2).  The  correlation  function  for  the  phase 
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can  be  derived  from  Eq.  (B2)  as 

«*«>«,)»=  fUdt 
J  0  •'0 

=  r(fI+fJ-|iI-fl|),  (B3) 

where  we  have  taken  #>(0)  =  0,  since  the  results 
calculated  later  cannot  depend  on  <p(0)  for  a 
stationary  process.  In  the  present  context,  the 
aim  of  the  model  is  to  calculate  correlation  func¬ 
tions  for  the  field  E(t)*&  ei0(,\  where  we  are 
using  the  notation  of  Eq.  (9).  Here  we  have  as¬ 
sumed  that  the  field  amplitude  &  (and  thus  the 
intensity)  is  a  constant.  So  the  correlation  func¬ 
tions  can  be  written 

«E(f ,)  •  •  • E(tm)E*(tml )  •  .•£*(/„)» 

=  S’i  exp[iV(f1)  + . . .  +  i<p(tm) 

^,) - MOD  •  <B4) 

In  order  to  calculate  these  correlations  it  is 
expedient  to  further  assume  that  the  phase  ip.(t) 
is  a  Gaussian  stochastic  quantity,  that  is,  cor¬ 
relation  functions  of  any  order  can  be  expressed 
in  terms  of  the  two-time  correlation  function  of  the 
phase  Specifically,”  we  have 

Mtt)  •••«’(*».,)>  =0, 

«*(*,)  -  •  •  <P(tJ»  -  T  ««>(«, >(*,. to  •  •  • 

9»r» 

(B5> 

where  the  summation  is  taken  over  all  unique 
permutations  of  t .,...,  fw.  A  useful  relation  can 
be  derived  from  Eq.  (B5),  which  makes  it  easy 
to  calculate  the  correlation  functions  in  Eq.  (B4). 
This  is 

((exP('/'  dt’Atwtyj 

=  exp ^  f  dt’J  dt’d(t')J{tm) 

x  «o<f'Ml'to) ,  (B6) 

where  J{t’)  is  an  arbitrary  function."  This  re¬ 
lation  can  be  proven,  term  by  term,  after  ex¬ 
panding  the  exponentials  and  using  the  property 
Eq.  (B5). 

We  can  now  calculate  the  desired  correlation 
functions.  By  letting  -tt)  in  Eq.  (B6) 

we  get 

((Ea.to-^explMO]) 

■  5  exp(-  17,)- 0,  (BY) 

where  we  have  taken  the  stationary  limit  n,  -  «, 
where  the  initial  transients  have  died  out.  Thus 
the  average  field  is  zero,  as  expected  for  a  fluc¬ 


tuating  field.  By  letting  J(t ')  *  6(/'  -  f,)  -  6(1'  -  f,) 
in  Eq.  (B6)  we  get 

«E(f  ,)£*(!,)»  =  S\{  expi  «>(f ,)  -  *>«,)]> 

=  5*exp(-r|f1-/1|).  (B8) 

The  power  spectrum  of  the  field,  given  by  the 
Fourier  transform  of  the  two-time  correlation 
function  in  Eq.  (B8),  is  thus  a  Lorentzian  with 
halfwidth  T.  The  four-time  correlation  function 
used  in  Sec.  IV  can  be  calculated  by  letting 

«/(*')  + 

-6(f'_/,)-6(f'-.#4), 
which  gives 
«£(f,)E(fj)£*(fj)E*(f4to 

=  $‘exp{r(|fI-f,|  +  |f3-/4|-|f,-/J| 

(B9) 

This  result  can  be  used  to  illustrate  one  of  the 
basic  assumptions  of  the  phase-diffusion  model. 

By  letting  /3  =  <,  and  /,  =  /,,  and  defining  the  in¬ 
tensity  as  /(f)  =  |  E(t)  | 2,  we  can  see  from  Eq.  (B9) 
that  the  intensity  correlation  function  is  given  by 

«W,)/(*2to-S4«/J,  (BIO) 

i.e.,  the  intensity  is  always  perfectly  correlated 
with  Itself  in  the  phase-diffusion  model,  because 
it  does  not  fluctuate. 

APPENDIX  C 

Here  we  describe  the  numerical  technique  used 
to  evaluate  the  double  integral  in  Eq.  (19b).  Let 
r  =  (T-T')*/!,  s-( r-r")1*,  and  a  =  (a*)l/*.  Then, 
we  have 

Ti/»  ,i t* 

i  |F(r)  |  to  *  4  /  dr  I  ds  e-^l^ar) 

J  o  J  o 

(Cl) 

*8  f  dr  f  d 
Jo  Jo 

x e  r,I/1(as )e'<rt*r*> ‘,*'**‘,  (C2) 

where  we  used  ihe  symmetry  of  the  integrand  in 
Eq.  (Cl)  with  respect  to  interchange  of  r  and  s. 
Now  defining 

«(*)-K|F(*‘)|to  (c3) 

and 

t>(x)  .e*«r'rt*rtM  f’ds  e-<r-rt-r,>.*/i(M)f 
"'0 

(C4) 
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we  have 


«(*)= /"  dr  I  t(ar)v(r) , 

J  a 

from  which  we  can  obtain 
u'(x)  =/,(<!*)«’(*), 


(C5) 

(C6) 


v'(x)m-2(r+rt  +  rs)xr(x)+e-’r’2/,M.  (C7) 

Thus  we  iiave  transformed  the  double  integral  into 
a  set  of  two  coupled  ordinary  differential  equations, 
Eqs.  (C6)  and  (C7),  which  can  be  solved  readily 
by  standard  numerical  techniques. 
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The  author*  show  that  the  laser  bandwidth  effects  in  resonance  fluorescence  at  low  intensities  can  be 
treated  most  readily  in  terms  of  simple  convolutions,  and  further  that  this  approach  allows  one  to  consider  an 
incident  field  with  an  arbitrary  spectral  distribution,  with  or  without  the  presence  of  atomic  (line-btoadenttn) 
collisions. 


1.  INTRODUCTION 

m  a  recent  paper  by  Knight,  Molander,  and 
Stroud1  (KMS)  an  interesting  pedagogical  discus¬ 
sion  was  given  of  the  effects  of  laser  bandwidth  on 
atomic  resonance  fluorescence  spectra  at  low  in¬ 
tensities.  Assuming  a  Lorentzian-shaped  laser 
spectrum  (phase-diffusion  model2*2)  they  showed 
that  the  finite  bandwidth  gives  rise  to  an  asymme¬ 
try  in  the  scattered  spectrum  when  the  incident 
field  is  detuned  from  the  atomic  resonance.  They 
offered  a  physical  explanation  for  the  asymmetry, 
first  found  by  Kimble  and  Mandel,4  in  terms  of  a 
continuing  reinitiation  of  the  transient  response  of 
the  atom  owing  to  the  fluctuations  in  the  laser  field. 
While  this  explanation  is  most  certainly  correct, 
one  point  remained  unclear,  that  concerning  the 
role  of  the  convolution  in  this  problem,  which  is  a 
particular  case  of  linear  response  theory. 

The  purpose  of  this  comment  is  twofold:  to  show 
that  the  laser  bandwidth  effects  in  resonance  fluo¬ 
rescence  at  low  intensities  can  be  treated  most 
readily  in  terms  of  a  simple  convolution  and  to 
_ _ I 


show  further  that  this  aporoach  allows  one  to  con¬ 
sider  an  incident  field  with  an  arbitrary  spectral 
distribution,  with  or  without  the  presence  of  atom¬ 
ic  (line-broadening)  collisions.  . 

II.  EFFECT  OF  LASER  UNE  SHAPE 

For  ease  of  comparison  we  will  follow  KMS  and 
use,  where  possible,  their  notation.  We  start  with 
the  low -intensity  operator  Bloch  equation  for  the 
positive-frequency  dipole  operator  d., 

«*(*)  =  (*w„  -  T„)o.(f)  -  ( id/K)t(t ) ,  (1) 

where  io0  is  the  atomic  resonance  frequency,  y„  is 
one-half  the  natural  decajr  width,  d  is  the  transi¬ 
tion  dipole  moment,  and  E'(t )  is  the  negative-fre¬ 
quency  electric  field  operator  for  the  driving  field, 
which  oscillates  with  frequency  uL.  The  scattered 
spectrum  g(u>,u L)  is  given  in  steady  state  by  the 
Fourier  transform  of  the  dipole  autocorrelation 
function  (GJ.(f)o-(f  +  t)K  which  is  obtained  by  for¬ 
mally  solving  Eq.  (1).  This  yields 


g(a>,vt)=j  drexp(»'u>T)|j  j  dt’  j*  dt*  expl(tco0  -  »)(!  -f')]expU-f«i>0 -y#)(f  +r -t”)]iE~(t')E*(t*)} ,  (2) 


where  o>  is  the  frequency  of  the  scattered  light. 

The  angular  brackets  &  >  indicate  both  a  quantum 
expectation  value  and  an  ensemble  average  over  the 
statistics  of  the  driving  field. 

For  a  monochromatic  laser  the  field-correlation 

i 


I - 

function  is 

«E-(f')E’(t')»  =  |  E0|2  explfwt(f' -<")] ,  (3) 

giving  for  the  narrow-band  spectrum  gna(<r>,vL) 


gj'rfrex p(iwT)rf*ffi'*  f  dt'  f"dt"  exp[(«w0 -  »)(f  -f')]expl(-  »u>0  -y„)(f  +  r  -r)]expl«u>t(f'  -I')]. 


(4) 


This  result  can  be  evaluated  to  give 

•<*-**>»  <5) 


l - 

KMS  treat  the  broad-band  laser  by  assuming  the 
phase -diffusion  model2*1  for  the  field 

«E'(f')E'(/')» 


which  shows  that  only  elastic  scattering  occurs 
for  an  isolated  atom  suffering  no  collisions. 


=  j  E„  |2 expl»wt(f '  -  r )Jexp(-  yL \t'  -  f  f ) , 


(«) 


j 

i 


20 


223A 


i«!  1979  The  American  Phyucel  Society 


RESONANCE  FLUORESCENCE  IN  A  WEAK  RADIATION  FIELD... 


2239 


20 

where  yL  is  one-half  the  laser  bandwidth.  Using 
this  in  Eq.  (2)  leads  directly  to  the  main  result  of 
KMS*: 

/  \  n*y,/2 

T**+6.->tn*Hv,+rtr\ 

x  /Aa+2A(w-h»,)+(r3-r?) 

\  (»-<*LY+yl 

(u,-w0)>+yf  )*  W 

where  the  detuning  is  a  =  w„ -ul  and  the  Rabi  fre¬ 
quency  is  fi  =  2d  I E0 1 H.  This  spectrum  clearly 
shows  two  maxima:  elastic  scattering  centered  at 
wL  and  fluorescence  centered  at  ui0.  However,  it 
also  seems  to  contain  interesting  interference 
structure,  e.g.,  (i*>  -  u>L)  in  the  numerator.  It  is 
stated  by  KMS  that  “the  spectrum  lf(w,(uL)]  is  not 
the  narrow-band  response  convolved  with  the  spec¬ 
tral  line  shape  of  the  incident  field”  (Lorentzlan  in 
this  case).  This  statement  is  certainly  correct  for 
intense,  saturating  incident  fields.  However,  we 
will  show  here  that  j(ki, <jl)  in  Eq.  (7),  valid  for 
low  fields,  is  obtainable  by  a  simple  convolution. 
This  is  as  expected  in  linear  response  theory. 

We  will  solve  the  problem  in  a  different,  more 
revealing,  way.  The  power  spectrum  PL(w'L )  of 
the  incident  field  is  generally  given  by* 

pi^D  dT  expfcXr)  iE'(t)b(t  +t)»  ,  (8)  . 

and  the  inverse  can  be  written  as 

)£*(/*)»  =/[<Kexp{i<(f'  -f*)]/*x(wj )  •  « 

Inserting  this  form  for  the  correlation  function 
into  Eq.  (2)  for  the  spectrum and  com¬ 
paring  the  resulting  equation  with  the  narrow -band 
spectrum  £wfl(w,wt)  in  Eq.  (4),  it  becomes  appar¬ 
ent  that  the  spectrum  g(w,wl)  can  be  written  as 

g(w,wL)  =  \£0\->  f^du’LgNB(»,u'L)PL(u'L).  (10) 

This  is  the  main  result  of  this  paper  and  shows 
that,  for  an  arbitrary  laser  spectrum  PL(u'L),  the 
scattered  spectrum  at  low  intensity,  is 

given  by  the  narrow-band  scattered  spectrum 
gHB{u,wL)  convolved  with  the  laser  spectrum.  Us¬ 
ing  the  8-function  form  of  g„a(w,wt)  from  Eq.  (5) 
in  Eq.  (10)  results  in 

g(«,u,t)..£0pl— (ID 

which  shows  that  the  scattered  light  retains  a  dis¬ 
torted  version  of  the  incident  light  spectrum 

For  the  special  case  of  the  phase-diffusion  mod¬ 


el,  used  by  KMS  and  described  in  Eq.  (6),  the  la¬ 
ser  spectrum  becomes  [by  inserting  Eq.  (6)  into 
Eq.  (8)J 


/ywi)=i^oi2 


yjr 

(»£ -»*)■+>*• 


(12) 


a  Lorentzian,  as  stated  earlier.  This  then  gives, 
from  Eq.  (11),  the  scattered  spectrum  as 

it  O2  yj* 

^w>wi)=2l(u^i2+y2j[(w-Wl)J+ylJ’  (13) 

a  result  which  is  simple  in  form,  but  in  fact  (after 
some  algebra)  is  identical  to  the  more  complex¬ 
looking  expression  given  by  KMS  in  Eq.  (7).  This 
shows  that  the  scattered  spectrum  is  the  product 
of  two  Lorentzians,  one  peaked  at  the  laser  fre¬ 
quency  wt  and  the  other  at  the  atomic  resonance 
frequency  u„. 


III.  COMMENTS 

We  have  shown  [Eq.  (10)]  that,  for  low-intensity 
resonance  fluorescence  in  the  absence  of  colli¬ 
sions,  the  scattered  spectrum  is  given  by  that  de¬ 
rived  by  assuming  a  monochromatic  driving  laser 
convolved  with  the  spectrum  of  the  actual  driving 
laser.  This  result  holds  for  arbitrary  laser  line 
shapes  and,  in  particular,  for  the  Lorentzian  line 
shape  derived  from  the  phase-diffusion  model1** 
and  used  by  KMS.  The  convolution  is  no  longer 
valid  when  an  intense,  saturating  driving  laser  is 
considered,  as  the  response  of  the  atom  is  then 
nonlinear.*  The  conclusion  of  KMS,  that  the  in¬ 
elastic  component  of  the  scattered  light  (centered 
on  u>0)  is  due  to  the  continuing  reinitiation  of  the 
transient  response  of  the  atom  by  fluctuations  in 
the  driving  field,  is  strengthened  by  the  present 
result.  The  conclusion  holds  regardless  of  the 
form  of  the  power  spectrum  of  the  fluctuations. 

Asa  further  demonstration  of  the  generality  of 
the  convolution  result,  we  consider  the  effects  of 
atomic  line-broadening  collisions  on  the  scattered 
spectrum.  These  include  dephasing  (or  elastic) 
collisions  and  quenching  collisions,  as  well  as 
alignment-  and  orientation-changing  collisions. 

The  effect  of  collisions  nns  been  studied  in  detail 
for  the  case  of  a  low- intensity  monochromatic 
driving  field.’’**  In  particular,  it  is  known  that  col¬ 
lisions  redistribute  the  frequency  of  the  scattered 
light  and  thus  Eq.  (5),  for  the  scattered  spectrum, 
is  not  appropriate  in  this  case.  However,  denoting 
the  appropriate  narrow-band  scattered  spectrum10 
by  Cyfl(w,wI),  the  scattered  spectrum  C(u),wt)  for 
brond-band  excitation  can  be  shown  still  to  be  a 
simple  convolution 

G(u.,«l)*IE0l**£d<..;c„(w,u»i)/*l(wi),  (14) 
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where  again  l’L{u'L)  is  the  (arbitrary)  spectrum  of 
the  low-intensity  driving  field.  This  result  can  be 
verified  by  inspecting  Eqs.  (4)  and  (5)  of  Hef.  7, 
which  give  the  scattered  spectrum  G,yB(u),  wL)  as 
an  integral  over  several  terms  including  the  factor 

f  ,(?t •  ?)(«  i  •  iOexpl— -  O), 

where  /,  is  the  intensity  of  the  incident  laser  with 
polarization  vector  «,  and  «T  is  the  atomic  dipole 
operator.  For  an  arbitrary  driving  field  spectrum 
this  factor  should  be  replaced  by 


where  £,(/)  is  the  electric  field  of  the  laser,  and 
the  field  average  has  been  performed.  This  field- 
correlation  function  is  analogous  to  the  one  ap¬ 
pearing  in  the  integrand  in  Eq.  (2)  for 
where  collisions  were  ignored.  Thus  the  deriva¬ 
tion  of  in  Eq.  (14)  follows  identically  to 

that  of /,'(u>,wt)  in  Eq.  (10).  The  convolution  re¬ 
sult  Eq.  (14)  for  the  scattered  spectrum  C(ut,u>L) 
holds  fcr  arbitrary  collisional  line  shapes,  both 
inside  and  outside  the  impact  region. 
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We  report  observation  of  the  scattered  spectrum  from  a  Zeeman  degenerate  atom  at  high  laser  intensities.  I-or  a  7  ■  0  to 
J  =  1  transition  in  the  presence  of  collisions  we  observe  an  asymmetric  triplet  in  the  polarization  parallel  to  the  incident 
laser  and  an  asymmetric  doublet  in  the  polarization  perpendicular  to  the  incident  laser. 


I.  Introduction.  In  recent  years  there  has  been 
widespread  interest  in  high-intei.sity  near-resonant 
light  scattering  in  the  presence  of  collisions  [1-6]. 
Mollow  [1,2]  has  developed  the  theory  for  a  simple 
two-state  atom,  including  the  effects  of  collisions  by 
introducing  phenomenological  relaxation  terms  to  the 
equations  of  a  strongly  driven  two-state  system.  This 
approximation,  which  is  equivalent  to  the  Markoff  ap¬ 
proximation,  leads  to  equations  which  are  valid  in  the 
impact  region  of  the  spectrum.  The  theory  predicts 
the  well-known  asymmetric  triplet  structure  in  the 
scattered  spectrum.  These  effects  have  been  observed 
by  Carlsten,  Szoke  and  Raymer  (CSR)  [4]  in  nonlin¬ 
ear  light  scattering  from  strontium  perturbed  by  argon. 

Recently,  Copper,  Ballagh  and  Burnett  (CBB)  [5], 
and  Fiutak  and  Van  Kranendonk  [6]  have  calculated 
the  expected  polarized  scattering  spectrum  at  high  in¬ 
tensities  when  Zeeman  degeneracy  effects  are  included. 
CBB  predict,  for  the  case  of  a  /  =  0  to  J  =  1  transition 
and  a  linearly  polarized  incident  laser,  an  asymmetric 
triplet  scattered  in  it  polarization,  and  an  asymmetric 
doublet  in  o  polarization.  The  predictions  for  the  it- 
polarized  scattering  spectrum  agree  with  ref.  [6],  how¬ 
ever,  Fiutak  and  Van  Kranendonk  do  not  quote  results 
for  the  o-polarized  spectrum.  In  this  letter  we  report 
the  observation  of  these  spectral  features  in  the  case 
of  strontium  perturbed  by  argon.  The  qualitative  agree¬ 
ment  with  the  results  of  CBB  |5]  is  i.e  good. 
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2.  Theoretical  considerations.  The  theory  for  in¬ 
tense  field  scattering  from  a  Zeeman  degenerate  atom 
in  a  collisional  environment  has  beer,  given  by  CBB 
[5]  in  the  impact  approximation  and  steady-state 
limit.  We  use  the  notation  of  CBB  throughout.  The  im¬ 
pact  approximation  restricts  the  validity  of  their  re¬ 
sults  to  pressures  low  enough  such  that  the  duration 
of  a  collision  (rc)  is  small  compared  to  the  time  be¬ 
tween  collisions  (  Tc).  In  addition  the  results  are  applic¬ 
able  for  detunings  from  line  center  (o>L  -  coy  =  A) 
small  compared  to  the  inverse  of  a  collision  duration, 

A  <  1  lrc.  For  the  case  of  van  der  Waals  broadening 
1/tc  *=  1012  s  1 ,  thus  their  results  should  be  valid  for 
detunings  A  <,5  cm-1 .  The  theory,  furthermore,  is 
valid  in  the  regime  where  the  collision  rates  are  inde¬ 
pendent  of  laser  intensity,  i.e.,  when  S2tc  <  1 ,  where 
ft  is  the  Rabi  frequency. 

Similarly  the  steady-state  assumption  restricts  their 
results  to  laser  pulse  lengths  (rp)  large  compared  to  the 
time  scales  associated  with  the  system  relaxation  rates 
(rp  >  Tc,  rN  where  tn  is  the  radiative  lifetime).  For 
buffer  gas  pressures  of  ^10  Tort  and  again  assuming  a 
van  der  Waals  interaction  Tc  «  tn  5  ns. 

For  the  case  of  interest  here,  a  linearly  polarized 
laser  tuned  near  the  Sr  resonance  line  (aJ  =  0  ioJ  =  1 
transition)  and  perturbed  by  Ar  (van  der  Waals  broad¬ 
ening)  CBB  [5]  predict  a  triplet  in  the  scattered  spec¬ 
trum  polarized  parallel  to  the  incident  laser  and  a 
doublet  in  the  scattered  spectrum  polarized  perpendic- 
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l'ig.  1 .  Processes  involved  in  the  formation  of  the  spectrum 
(the  straight  arrows  represent  absorption  of  a  laser  photon, 
the  wavy  arrows  represent  an  emitted  photon,  and  the  double 
arrows  represent  a  collisional  process). 

ular  to  the  incident  laser.  In  the  parallel  polarization 
(fig.  la)  the  three  peaks  correspond  to  (i)  Rayleigh 
scattering  off  the  ground  or  my  =  0  sublevel  of  the  ex¬ 
cited  state,  (ii)  “three-photon”  scattering  followed  by 
fluorescence,  and  (iii)  collision-induced  fluorescence 
from  the  excited  mj  =  0  state.  These  features  are  exact¬ 
ly  those  associated  with  nonlinear  scattering  from  an 
idealized  two-state  atom  in  the  presence  of  collisions. 

In  the  perpendicular  polarization  the  two  peaks  cor¬ 
respond  to  collisional  excitation  of  the  ttij-t]  sub- 
levels  of  the  excited  state  (fig.  lb)  followed  by  (i)  fluo¬ 
rescence  to  the  ground  state  or  (ii)  an  inverse  Raman 
excitation  of  the  my  =  0  excited  state. 

In  the  experimental  results  reported  here,  neither 
the  impact  approximation  nor  the  steady-state  approx¬ 
imation  is  strictly  valid.  Thus  detailed  quantitative 
comparisons  with  the  theory  are  not  possible.  How¬ 
ever.  the  experimental  results  are  in  agreement  with 
the  qualitative  conclusions  of  CBB  (51. 

3.  Experimental  apparatus.  The  experimental  appa¬ 
ratus  is  basically  the  same  as  that  used  by  CSR  [4] . 

The  apparatus  consists  of  an  N2-pumped  Hansch-type 
dye  laser  tuned  near  the  strontium  resonance  line  (6s2 
•Sfl-bsOp  1  Pj  at  460.7  nm)  and  tightly  focused  into 
a  stainless-steel  crossed  oven  containing  strontium  va¬ 
por  at  =»5008C-600°C  and  an  argon  buffer  gas  at  pres¬ 
sures  of  *2  to  500  Ton.  The  laser  power  at  the  center 
of  the  oven  was  estimated  to  be  *=1  GW/cm2.  The  scat¬ 


tered  light  was  detected  at  90°  through  the  side  arm 
of  the  oven,  and  focused  onto  the  entrance  slit  of  the 
monochromator.  The  signal  was  detected  by  a  photo¬ 
multiplier,  fed  into  an  rf  shielded,  gated,  analog  inte¬ 
grate  and  hold  circuit  and  recorded  on  a  chart  recorder 
with  a  2.2  s  time  constant.  Further  details  of  the  appa¬ 
ratus  can  be  found  in  CSR  [4).  The  major  difference 
in  our  experiment  was  the  addition  of  a  glan-prism  po¬ 
larizer  on  the  dye-laser  output,  an  analyzing  polarizer 
at  the  side  window  output  of  the  oven,  and  a  reference 
polarizer  oriented  at  45°  and  positioned  before  the  en¬ 
trance  slit  of  the  monochromator  to  remove  the  polar¬ 
izing  influence  of  the  grating. 

4.  Experimental  results;  spectrum.  Typical  spectra 
are  shown  in  fig.  2.  These  spectra  were  taken  at  a 
temperature  of  =»600°C  and  argon  pressures  of  5,  50 
and  500  Torr,  respectively.  The  Sr  density,  as  mea- 
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Fig.  2.  Typical  scattered  light  spectra  (red  detuning),  at  (a! 

5  Ton,  (b)  50  Ton  and  (c)  500  Ton.  In  each  case  the  upper 
curve  is  the  parallel  polarization  and  the  lower  curve  the  per¬ 
pendicular  polarization.  The  anows  indicate  the  unshifted  po¬ 
sitions  of  the  peaks. 
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